Changeset 8920 for branches/2017/dev_CNRS_2017/NEMOGCM/NEMO/OPA_SRC

Timestamp:

2017-12-06T15:33:57+01:00 (6 years ago)

Author:

clem

Message:

debug ice option in which conductivity is calculated from bulk formulae and used in ice thermodynamics (np_jules_active in forced mode)

File:

• branches/2017/dev_CNRS_2017/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk.F90 (modified) (3 diffs)

branches/2017/dev_CNRS_2017/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk.F90

@@ -68 +68 @@
    PUBLIC   sbc_blk       ! called in sbcmod
 #if defined key_lim3
-   PUBLIC   blk_ice_tau   ! routine called in icestp module
-   PUBLIC   blk_ice_flx   ! routine called in icestp module
-   PUBLIC   blk_ice_qcn   ! routine called in icestp module
+   PUBLIC   blk_ice_tau   ! routine called in iceforcing
+   PUBLIC   blk_ice_flx   ! routine called in iceforcing
+   PUBLIC   blk_ice_qcn   ! routine called in iceforcing
 #endif 
 
@@ -1041 +1041 @@
       CASE ( 1 , 3 )
          !
-         zfac    =   1._wp /  ( rn_cnd_s + rcdic )
-         zfac2   =   EXP(1._wp) * 0.5_wp * zepsilon
-         zfac3   =   2._wp / zepsilon
+         zfac  = 1._wp /  ( rn_cnd_s + rcdic )
+         zfac2 = EXP(1._wp) * 0.5_wp * zepsilon
+         zfac3 = 2._wp / zepsilon
          !   
          DO jl = 1, jpl                
             DO jj = 1 , jpj
                DO ji = 1, jpi
-                  !                                ! Effective thickness
-                  zhe               =   ( rn_cnd_s * phi(ji,jj,jl) + rcdic * phs(ji,jj,jl) ) * zfac
-                  !                                ! Enhanced conduction factor
-                  IF( zhe >=  zfac2 )  &
-                     zgfac(ji,jj,jl)   =   MIN( 2._wp, ( 0.5_wp + 0.5 * LOG( zhe * zfac3 ) ) )
+                  zhe = ( rn_cnd_s * phi(ji,jj,jl) + rcdic * phs(ji,jj,jl) ) * zfac                             ! Effective thickness
+                  IF( zhe >=  zfac2 )   zgfac(ji,jj,jl) = MIN( 2._wp, 0.5_wp * ( 1._wp + LOG( zhe * zfac3 ) ) ) ! Enhanced conduction factor
                END DO
             END DO
@@ -1063 +1060 @@
       ! -------------------------------------------------------------!
       !
-      zfac   =   rcdic * rn_cnd_s 
-      !                             ! ========================== !
-      DO jl = 1, jpl                !  Loop over ice categories  !
-         !                          ! ========================== !
+      zfac = rcdic * rn_cnd_s
+      !
+      DO jl = 1, jpl
          DO jj = 1 , jpj
             DO ji = 1, jpi
-               !                    ! Effective conductivity of the snow-ice system divided by thickness
-               zkeff_h =   zfac * zgfac(ji,jj,jl) / ( rcdic * phs(ji,jj,jl) + rn_cnd_s * phi(ji,jj,jl) )
-               !                    ! Store initial surface temperature
-               ztsu    =   ptsu(ji,jj,jl)
-               !                    ! Net initial atmospheric heat flux
-               zqa0    =   qsr_ice(ji,jj,jl) - qsr_ice_tr(ji,jj,jl) + qns_ice(ji,jj,jl)
+               !                    
+               zkeff_h = zfac * zgfac(ji,jj,jl) / &                                    ! Effective conductivity of the snow-ice system divided by thickness
+                  &      ( rcdic * phs(ji,jj,jl) + rn_cnd_s * MAX( 0.01, phi(ji,jj,jl) ) )
+               ztsu    = ptsu(ji,jj,jl)                                                ! Store initial surface temperature
+               zqa0    = qsr_ice(ji,jj,jl) - qsr_ice_tr(ji,jj,jl) + qns_ice(ji,jj,jl)  ! Net initial atmospheric heat flux
                !
-               DO iter = 1, nit     ! --- Iteration loop
-                   !                      ! Conduction heat flux through snow-ice system (>0 downwards)
-                   zqc   =   zkeff_h * ( ztsu - ptb(ji,jj) )
-                   !                      ! Surface energy budget
-                   zqnet =   zqa0 + dqns_ice(ji,jj,jl) * ( ztsu - ptsu(ji,jj,jl) ) - zqc
-                   !                      ! Temperature update
-                   ztsu  =   ztsu - zqnet / ( dqns_ice(ji,jj,jl) - zkeff_h )
+               DO iter = 1, nit     ! --- Iterative loop
+                  zqc   = zkeff_h * ( ztsu - ptb(ji,jj) )                              ! Conduction heat flux through snow-ice system (>0 downwards)
+                  zqnet = zqa0 + dqns_ice(ji,jj,jl) * ( ztsu - ptsu(ji,jj,jl) ) - zqc  ! Surface energy budget
+                  ztsu  = ztsu - zqnet / ( dqns_ice(ji,jj,jl) - zkeff_h )              ! Temperature update
                END DO
                !