Changeset 13169
- Timestamp:
- 2020-06-27T19:10:19+02:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/r4.0-HEAD_r12713_clem_dan_fixcpl/src/ICE/icethd_zdf_bl99.F90
r13148 r13169 133 133 REAL(wp), DIMENSION(jpij) :: za_s_fra ! ice fraction covered by snow 134 134 REAL(wp), DIMENSION(jpij) :: isnow ! snow presence (1) or not (0) 135 REAL(wp), DIMENSION(jpij) :: isnow_comb ! snow presence for met-office 135 136 ! 136 137 ! Mono-category … … 182 183 isnow (ji) = 0._wp 183 184 ENDIF 185 ! for Met-Office 186 IF( h_s_1d(ji) < zh_min ) THEN 187 isnow_comb(ji) = h_s_1d(ji) / zh_min 188 ELSE 189 isnow_comb(ji) = 1._wp 190 ENDIF 184 191 END DO 185 192 ! clem: we should apply correction on snow thickness to take into account snow fraction … … 331 338 DO ji = 1, npti ! Snow-ice interface 332 339 IF ( .NOT. l_T_converged(ji) ) THEN 333 ! Calculate combined surface snow and ice conductivity to pass through the coupler 334 zkappa_comb(ji) = isnow (ji) * zkappa_s(ji,0) + ( 1._wp - isnow(ji) ) * zkappa_i(ji,0)340 ! Calculate combined surface snow and ice conductivity to pass through the coupler (met-office) 341 zkappa_comb(ji) = isnow_comb(ji) * zkappa_s(ji,0) + ( 1._wp - isnow_comb(ji) ) * zkappa_i(ji,0) 335 342 ! If there is snow then use the same snow-ice interface conductivity for the top layer of ice 336 343 IF( h_s_1d(ji) > 0._wp ) zkappa_i(ji,0) = zkappa_s(ji,nlay_s) … … 896 903 DO ji = 1, npti 897 904 IF( h_i_1d(ji) >= zhi_ssl ) THEN 898 !!clem metOcnd_ice_1d(ji) = 2._wp * zkappa_comb(ji)899 cnd_ice_1d(ji) = 2._wp * zkappa_i(ji,0)905 cnd_ice_1d(ji) = 2._wp * zkappa_comb(ji) 906 !!cnd_ice_1d(ji) = 2._wp * zkappa_i(ji,0) 900 907 ELSE 901 908 cnd_ice_1d(ji) = 2._wp * ztcond_i(ji,0) / zhi_ssl ! cnd_ice is capped by: cond_i/zhi_ssl
Note: See TracChangeset
for help on using the changeset viewer.