Changeset 3807 for branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/NEMO/LIM_SRC_3/limthd_dh.F90

Timestamp:

2013-02-13T07:29:43+01:00 (11 years ago)

Author:

gm

Message:

dev_v3_4_STABLE_2012: #997 and #898: in LIM3, no use of qla_ice in coupled mode and its capping in forced mode

File:

• branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/NEMO/LIM_SRC_3/limthd_dh.F90 (modified) (2 diffs)

branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/NEMO/LIM_SRC_3/limthd_dh.F90

@@ -38 +38 @@
 
    !!----------------------------------------------------------------------
-   !! NEMO/LIM3 4.0 , UCL - NEMO Consortium (2010)
+   !! NEMO/LIM3 3.4 , UCL - NEMO Consortium (2010)
    !! $Id$
    !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt)
@@ -317 +317 @@
 
       DO ji = kideb, kiut
-         ! if qla is positive (upwards), heat goes to the atmosphere, therefore
-         ! snow sublimates, if qla is negative (downwards), snow condensates
+         ! qla_ice is always >=0 (upwards), heat goes to the atmosphere, therefore snow sublimates
+#if defined key_coupled
+         zdh_s_sub(ji)    =  0._wp      ! coupled mode: sublimation already included in emp_ice (to do in limsbc_ice)
+#else
+         !                              ! forced  mode: snow thickness change due to sublimation
          zdh_s_sub(ji)    =  - parsub * qla_ice_1d(ji) / ( rhosn * lsub ) * rdt_ice
+#endif
          dh_s_tot (ji)    =  dh_s_tot(ji) + zdh_s_sub(ji)
          zdhcf            =  ht_s_b(ji) + zdh_s_sub(ji)