Changeset 12854 for NEMO/branches/2020/r4.0-HEAD_r12713_clem_dan_fixcpl/src/ICE/icethd_zdf_bl99.F90

Timestamp:

2020-05-01T19:06:23+02:00 (4 years ago)

Author:

clem

Message:

1st implementation of snow fraction (impact on albedo). Light transmission is still not ok since we need a non-zero penetration of solar flux when snow thickness is > 0

File:

• NEMO/branches/2020/r4.0-HEAD_r12713_clem_dan_fixcpl/src/ICE/icethd_zdf_bl99.F90 (modified) (5 diffs)

NEMO/branches/2020/r4.0-HEAD_r12713_clem_dan_fixcpl/src/ICE/icethd_zdf_bl99.F90

@@ -130 +130 @@
       REAL(wp), DIMENSION(jpij)            ::   zq_ini      ! diag errors on heat
       REAL(wp), DIMENSION(jpij)            ::   zghe        ! G(he), th. conduct enhancement factor, mono-cat
+      REAL(wp), DIMENSION(jpij)            ::   za_s_fra    ! ice fraction covered by snow 
       !
       ! Mono-category
@@ -143 +144 @@
       END DO
 
+      CALL ice_var_snwfra( h_s_1d(1:npti), za_s_fra(1:npti) )   ! calculate ice fraction covered by snow
+      
       !------------------
       ! 1) Initialization
@@ -149 +152 @@
 
          ! If the snow thickness drops below zhs_min then reduce the snow fraction instead
-         IF( h_s_1d(ji) < zhs_min ) THEN
-           isnow(ji) = h_s_1d(ji) / zhs_min
-           zh_s(ji) = zhs_min * r1_nlay_s
-         ELSE
-           isnow(ji) = 1.0_wp
-           zh_s(ji) = h_s_1d(ji) * r1_nlay_s
-         END IF
-
+         !!IF( h_s_1d(ji) < zhs_min ) THEN
+         !!  isnow(ji) = h_s_1d(ji) / zhs_min
+         !!  zh_s(ji) = zhs_min * r1_nlay_s
+         !!ELSE
+         !!  isnow(ji) = 1.0_wp
+         !!  zh_s(ji) = h_s_1d(ji) * r1_nlay_s
+         !!END IF
+         isnow(ji) = za_s_fra(ji) !!clem: 2 variables for the same thing
+         zh_s(ji) = h_s_1d(ji) * r1_nlay_s
+         
+         
          ! layer thickness
          zh_i(ji) = h_i_1d(ji) * r1_nlay_i
@@ -166 +172 @@
       END WHERE
       !
-      WHERE( zh_s(1:npti) > 0._wp   )   ;   z1_h_s(1:npti) = 1._wp / zh_s(1:npti)
+      WHERE( zh_s(1:npti) >= epsi10 )   ;   z1_h_s(1:npti) = 1._wp / zh_s(1:npti)
       ELSEWHERE                         ;   z1_h_s(1:npti) = 0._wp
       END WHERE
@@ -881 +887 @@
          !--- Snow-ice interfacial temperature (diagnostic SIMIP)
          zfac = rn_cnd_s * zh_i(ji) + ztcond_i(ji,1) * zh_s(ji)
-         IF( h_s_1d(ji) >= zhs_min ) THEN
+         IF( h_s_1d(ji) >= zhs_min ) THEN !!clem change that
             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) ) / MAX( epsi10, zfac )