Changeset 3807 for branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/NEMO/LIM_SRC_3/limthd_dif.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_dif.F90 (modified) (2 diffs)

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

@@ -34 +34 @@
 
    !!----------------------------------------------------------------------
-   !! NEMO/LIM3 4.0 , UCL - NEMO Consortium (2011)
+   !! NEMO/LIM3 3.4 , UCL - NEMO Consortium (2011)
    !! $Id$
    !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt)
@@ -717 +717 @@
       !   11) Fluxes at the interfaces                                          !
       !-------------------------------------------------------------------------!
+         ! duplicate the loop for performances issues
       DO ji = kideb, kiut
-         !                                ! update of latent heat fluxes
-         qla_ice_1d (ji) = qla_ice_1d (ji) + dqla_ice_1d(ji) * ( t_su_b(ji) - ztsuold(ji) )
+#if ! defined key_coupled
+         ! forced mode only : update of latent heat fluxes (sublimation) (always >=0, upward flux) 
+         qla_ice_1d (ji) = MAX( 0._wp, qla_ice_1d (ji) + dqla_ice_1d(ji) * ( t_su_b(ji) - ztsuold(ji) ) )
+#endif
          !                                ! surface ice conduction flux
          isnow(ji)       = INT(  1._wp - MAX( 0._wp, SIGN( 1._wp, -ht_s_b(ji) ) )  )