 Timestamp:
 20200530T08:17:37+02:00 (4 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

NEMO/branches/UKMO/NEMO_4.0.1_remove_0.1m_snow_test/src/ICE/icethd_zdf_bl99.F90
r12834 r12995 124 124 REAL(wp), DIMENSION(jpij,0:nlay_s) :: zkappa_s ! Kappa factor in the snow 125 125 REAL(wp), DIMENSION(jpij,0:nlay_s) :: zeta_s ! Eta factor in the snow 126 REAL(wp), DIMENSION(jpij) :: zkappa_combi ! Combined snow and ice surface conductivity 126 127 REAL(wp), DIMENSION(jpij,nlay_i+3) :: zindterm ! 'Ind'ependent term 127 128 REAL(wp), DIMENSION(jpij,nlay_i+3) :: zindtbis ! Temporary 'ind'ependent term … … 317 318 END DO 318 319 DO ji = 1, npti ! Snowice interface 319 IF ( .NOT. l_T_converged(ji) ) & 320 zkappa_i(ji,0) = zkappa_s(ji,nlay_s) * isnow(ji) + zkappa_i(ji,0) * ( 1._wp  isnow(ji) ) 320 IF ( .NOT. l_T_converged(ji) ) THEN 321 ! Calculate combined surface snow and ice conductivity to pass through the coupler 322 zkappa_combi(ji) = zkappa_s(ji,0) * isnow(ji) + zkappa_i(ji,0) * ( 1._wp  isnow(ji) ) 323 324 ! If there is snow then use the same snowice interface conductivity for the top layer of ice 325 IF ( h_s_1d(ji) > 0._wp ) THEN 326 zkappa_i(ji,0) = zkappa_s(ji,nlay_s) 327 END IF 328 321 329 END DO 322 330 ! … … 847 855 DO ji = 1, npti 848 856 IF( h_i_1d(ji) > 0.1_wp ) THEN 849 cnd_ice_1d(ji) = 2._wp * zkappa_ i(ji,0)857 cnd_ice_1d(ji) = 2._wp * zkappa_combi(ji) 850 858 ELSE 851 859 cnd_ice_1d(ji) = 2._wp * ztcond_i(ji,0) * 10._wp
Note: See TracChangeset
for help on using the changeset viewer.