Changeset 1849 for LMDZ5/trunk/libf/phylmd/physiq.F

Timestamp:

Aug 27, 2013, 12:55:18 PM (11 years ago)

Author:

Ehouarn Millour

Message:

Including the thermodynamics of ice in the convection scheme (iactive only if iflag_ice_thermo==1).
CR+JYG

File:

• LMDZ5/trunk/libf/phylmd/physiq.F (modified) (3 diffs)

LMDZ5/trunk/libf/phylmd/physiq.F

@@ -1878 +1878 @@
       DO i = 1, klon
          zlvdcp=RLVTT/RCPD/(1.0+RVTMP2*q_seri(i,k))
-c        zlsdcp=RLSTT/RCPD/(1.0+RVTMP2*q_seri(i,k))
-         zlsdcp=RLVTT/RCPD/(1.0+RVTMP2*q_seri(i,k))
+cjyg<
+c      Attention : Arnaud a propose des formules completement differentes
+c                  A verifier !!!
+         zlsdcp=RLSTT/RCPD/(1.0+RVTMP2*q_seri(i,k))
+         IF (iflag_ice_thermo .EQ. 0) THEN
+           zlsdcp=zlvdcp
+         ENDIF 
+c>jyg
+
          zdelta = MAX(0.,SIGN(1.,RTT-t_seri(i,k)))
          zb = MAX(0.0,ql_seri(i,k))
@@ -2263 +2270 @@
                 nbtr_tmp=nbtr
              END IF
-          CALL concvl (iflag_con,iflag_clos,
+cjyg   iflag_con est dans clesphys
+cc          CALL concvl (iflag_con,iflag_clos,
+          CALL concvl (iflag_clos,
      .        dtime,paprs,pplay,t_undi,q_undi,
      .        t_wake,q_wake,wake_s,
@@ -2765 +2774 @@
      .           frac_impa, frac_nucl, beta_prec_fisrt,
      .           prfl, psfl, rhcl, 
-     .           zqasc, fraca,ztv,zpspsk,ztla,zthl,iflag_cldcon )
+     .           zqasc, fraca,ztv,zpspsk,ztla,zthl,iflag_cldcon,
+     .           iflag_ice_thermo)
 
       WHERE (rain_lsc < 0) rain_lsc = 0.