Changeset 9423 for branches/2017/dev_merge_2017/NEMOGCM

Timestamp:

2018-03-22T20:38:54+01:00 (6 years ago)

Author:

clem

Message:

small correction for the temperature diffusion when snow thickness is small

File:

• branches/2017/dev_merge_2017/NEMOGCM/NEMO/LIM_SRC_3/icethd_zdf_bl99.F90 (modified) (4 diffs)

branches/2017/dev_merge_2017/NEMOGCM/NEMO/LIM_SRC_3/icethd_zdf_bl99.F90

@@ -95 +95 @@
       REAL(wp) ::   ztsu_err  =  1.e-5_wp     ! range around which t_su is considered at 0C 
       REAL(wp) ::   zdti_bnd  =  1.e-4_wp     ! maximal authorized error on temperature 
+      REAL(wp) ::   zhs_min   =  0.01_wp      ! minimum snow thickness for conductivity calculation 
       REAL(wp) ::   ztmelt_i                  ! ice melting temperature
       REAL(wp) ::   zdti_max                  ! current maximal error on temperature 
@@ -156 +157 @@
       END WHERE
       !
-      WHERE( zh_s(1:npti) >= epsi10 )   ;   z1_h_s(1:npti) = 1._wp / zh_s(1:npti)
+      WHERE( zh_s(1:npti) > 0._wp   )       zh_s(1:npti) = MAX( zhs_min * r1_nlay_s, zh_s(1:npti) )
+      !
+      WHERE( zh_s(1:npti) > 0._wp   )   ;   z1_h_s(1:npti) = 1._wp / zh_s(1:npti)
       ELSEWHERE                         ;   z1_h_s(1:npti) = 0._wp
       END WHERE
@@ -182 +185 @@
          DO ji = 1, npti
             !                             ! radiation transmitted below the layer-th snow layer
-            zradtr_s(ji,jk) = zradtr_s(ji,0) * EXP( - zraext_s * zh_s(ji) * REAL(jk) )
+            zradtr_s(ji,jk) = zradtr_s(ji,0) * EXP( - zraext_s * h_s_1d(ji) * r1_nlay_s * REAL(jk) )
             !                             ! radiation absorbed by the layer-th snow layer
             zradab_s(ji,jk) = zradtr_s(ji,jk-1) - zradtr_s(ji,jk)
@@ -798 +801 @@
             !--- Snow-ice interfacial temperature (diagnostic SIMIP)
             zfac = rn_cnd_s * zh_i(ji) + ztcond_i(ji,1) * zh_s(ji)
-            IF( zh_s(ji) >= 1.e-3 .AND. zfac > epsi10 ) THEN
+            IF( h_s_1d(ji) >= zhs_min ) THEN
                t_si_1d(ji) = ( rn_cnd_s       * zh_i(ji) * t_s_1d(ji,1) +   &
-                  &            ztcond_i(ji,1) * zh_s(ji) * t_i_1d(ji,1) ) / zfac
+                  &            ztcond_i(ji,1) * zh_s(ji) * t_i_1d(ji,1) ) / MAX( epsi10, zfac )
             ELSE
                t_si_1d(ji) = t_su_1d(ji)