Changeset 991 for branches/dev_003_CPL/NEMO/LIM_SRC_3/limthd_dh.F90

Timestamp:

2008-05-23T17:55:55+02:00 (16 years ago)

Author:

smasson

Message:

dev_003_CPL: preliminary draft (not working), see ticket #155

File:

• branches/dev_003_CPL/NEMO/LIM_SRC_3/limthd_dh.F90 (modified) (1 diff)

branches/dev_003_CPL/NEMO/LIM_SRC_3/limthd_dh.F90

@@ -384 +384 @@
          ! if qla is positive (upwards), heat goes to the atmosphere, therefore
          ! snow sublimates, if qla is negative (downwards), snow condensates
-         zdh_s_sub(ji)   =  - parsub * qla_ice_1d(ji) / ( rhosn * lsub ) * rdt_ice
+         
+         IF ( .NOT. lk_cpl ) THEN 
+            zdh_s_sub(ji) = - parsub * qla_ice_1d(ji) / ( rhosn * lsub ) * rdt_ice
+         ELSE IF (parsub == 1) THEN 
+            CALL ctl_stop( 'In coupled mode, use parsub = 0 or send dqla' ) 
+         ELSE
+            zdh_s_sub(ji)  = 0.0
+         ENDIF
          dh_s_tot(ji)    =  dh_s_tot(ji) + zdh_s_sub(ji)
          zdhcf           =  ht_s_b(ji) + zdh_s_sub(ji)