Changeset 3808
- Timestamp:
- 2013-02-13T07:31:26+01:00 (11 years ago)
- Location:
- branches/2012/dev_MERGE_2012/NEMOGCM/NEMO/LIM_SRC_3
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2012/dev_MERGE_2012/NEMOGCM/NEMO/LIM_SRC_3/limthd_dh.F90
r3625 r3808 315 315 316 316 DO ji = kideb, kiut 317 ! if qla is positive (upwards), heat goes to the atmosphere, therefore 318 ! snow sublimates, if qla is negative (downwards), snow condensates 317 ! qla_ice is always >=0 (upwards), heat goes to the atmosphere, therefore snow sublimates 318 #if defined key_coupled 319 zdh_s_sub(ji) = 0._wp ! coupled mode: sublimation already included in emp_ice (to do in limsbc_ice) 320 #else 321 ! ! forced mode: snow thickness change due to sublimation 319 322 zdh_s_sub(ji) = - parsub * qla_ice_1d(ji) / ( rhosn * lsub ) * rdt_ice 323 #endif 320 324 dh_s_tot (ji) = dh_s_tot(ji) + zdh_s_sub(ji) 321 325 zdhcf = ht_s_b(ji) + zdh_s_sub(ji) -
branches/2012/dev_MERGE_2012/NEMOGCM/NEMO/LIM_SRC_3/limthd_dif.F90
r3625 r3808 716 716 !-------------------------------------------------------------------------! 717 717 DO ji = kideb, kiut 718 ! ! update of latent heat fluxes 719 qla_ice_1d (ji) = qla_ice_1d (ji) + dqla_ice_1d(ji) * ( t_su_b(ji) - ztsuold(ji) ) 718 #if ! defined key_coupled 719 ! forced mode only : update of latent heat fluxes (sublimation) (always >=0, upward flux) 720 qla_ice_1d (ji) = MAX( 0._wp, qla_ice_1d (ji) + dqla_ice_1d(ji) * ( t_su_b(ji) - ztsuold(ji) ) ) 721 #endif 720 722 ! ! surface ice conduction flux 721 723 isnow(ji) = INT( 1._wp - MAX( 0._wp, SIGN( 1._wp, -ht_s_b(ji) ) ) )
Note: See TracChangeset
for help on using the changeset viewer.