Changeset 3807 for branches/2012/dev_v3_4_STABLE_2012/NEMOGCM
- Timestamp:
- 2013-02-13T07:29:43+01:00 (11 years ago)
- Location:
- branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/NEMO/LIM_SRC_3
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/NEMO/LIM_SRC_3/limthd_dh.F90
r3801 r3807 38 38 39 39 !!---------------------------------------------------------------------- 40 !! NEMO/LIM3 4.0, UCL - NEMO Consortium (2010)40 !! NEMO/LIM3 3.4 , UCL - NEMO Consortium (2010) 41 41 !! $Id$ 42 42 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) … … 317 317 318 318 DO ji = kideb, kiut 319 ! if qla is positive (upwards), heat goes to the atmosphere, therefore 320 ! snow sublimates, if qla is negative (downwards), snow condensates 319 ! qla_ice is always >=0 (upwards), heat goes to the atmosphere, therefore snow sublimates 320 #if defined key_coupled 321 zdh_s_sub(ji) = 0._wp ! coupled mode: sublimation already included in emp_ice (to do in limsbc_ice) 322 #else 323 ! ! forced mode: snow thickness change due to sublimation 321 324 zdh_s_sub(ji) = - parsub * qla_ice_1d(ji) / ( rhosn * lsub ) * rdt_ice 325 #endif 322 326 dh_s_tot (ji) = dh_s_tot(ji) + zdh_s_sub(ji) 323 327 zdhcf = ht_s_b(ji) + zdh_s_sub(ji) -
branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/NEMO/LIM_SRC_3/limthd_dif.F90
r3558 r3807 34 34 35 35 !!---------------------------------------------------------------------- 36 !! NEMO/LIM3 4.0, UCL - NEMO Consortium (2011)36 !! NEMO/LIM3 3.4 , UCL - NEMO Consortium (2011) 37 37 !! $Id$ 38 38 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) … … 717 717 ! 11) Fluxes at the interfaces ! 718 718 !-------------------------------------------------------------------------! 719 ! duplicate the loop for performances issues 719 720 DO ji = kideb, kiut 720 ! ! update of latent heat fluxes 721 qla_ice_1d (ji) = qla_ice_1d (ji) + dqla_ice_1d(ji) * ( t_su_b(ji) - ztsuold(ji) ) 721 #if ! defined key_coupled 722 ! forced mode only : update of latent heat fluxes (sublimation) (always >=0, upward flux) 723 qla_ice_1d (ji) = MAX( 0._wp, qla_ice_1d (ji) + dqla_ice_1d(ji) * ( t_su_b(ji) - ztsuold(ji) ) ) 724 #endif 722 725 ! ! surface ice conduction flux 723 726 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.