 Timestamp:
 20200501T19:06:23+02:00 (4 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

NEMO/branches/2020/r4.0HEAD_r12713_clem_dan_fixcpl/src/ICE/icethd_zdf_bl99.F90
r12827 r12854 130 130 REAL(wp), DIMENSION(jpij) :: zq_ini ! diag errors on heat 131 131 REAL(wp), DIMENSION(jpij) :: zghe ! G(he), th. conduct enhancement factor, monocat 132 REAL(wp), DIMENSION(jpij) :: za_s_fra ! ice fraction covered by snow 132 133 ! 133 134 ! Monocategory … … 143 144 END DO 144 145 146 CALL ice_var_snwfra( h_s_1d(1:npti), za_s_fra(1:npti) ) ! calculate ice fraction covered by snow 147 145 148 ! 146 149 ! 1) Initialization … … 149 152 150 153 ! If the snow thickness drops below zhs_min then reduce the snow fraction instead 151 IF( h_s_1d(ji) < zhs_min ) THEN 152 isnow(ji) = h_s_1d(ji) / zhs_min 153 zh_s(ji) = zhs_min * r1_nlay_s 154 ELSE 155 isnow(ji) = 1.0_wp 156 zh_s(ji) = h_s_1d(ji) * r1_nlay_s 157 END IF 158 154 !!IF( h_s_1d(ji) < zhs_min ) THEN 155 !! isnow(ji) = h_s_1d(ji) / zhs_min 156 !! zh_s(ji) = zhs_min * r1_nlay_s 157 !!ELSE 158 !! isnow(ji) = 1.0_wp 159 !! zh_s(ji) = h_s_1d(ji) * r1_nlay_s 160 !!END IF 161 isnow(ji) = za_s_fra(ji) !!clem: 2 variables for the same thing 162 zh_s(ji) = h_s_1d(ji) * r1_nlay_s 163 164 159 165 ! layer thickness 160 166 zh_i(ji) = h_i_1d(ji) * r1_nlay_i … … 166 172 END WHERE 167 173 ! 168 WHERE( zh_s(1:npti) > 0._wp) ; z1_h_s(1:npti) = 1._wp / zh_s(1:npti)174 WHERE( zh_s(1:npti) >= epsi10 ) ; z1_h_s(1:npti) = 1._wp / zh_s(1:npti) 169 175 ELSEWHERE ; z1_h_s(1:npti) = 0._wp 170 176 END WHERE … … 881 887 ! Snowice interfacial temperature (diagnostic SIMIP) 882 888 zfac = rn_cnd_s * zh_i(ji) + ztcond_i(ji,1) * zh_s(ji) 883 IF( h_s_1d(ji) >= zhs_min ) THEN 889 IF( h_s_1d(ji) >= zhs_min ) THEN !!clem change that 884 890 t_si_1d(ji) = ( rn_cnd_s * zh_i(ji) * t_s_1d(ji,1) + & 885 891 & ztcond_i(ji,1) * zh_s(ji) * t_i_1d(ji,1) ) / MAX( epsi10, zfac )
Note: See TracChangeset
for help on using the changeset viewer.