Changeset 2448

Timestamp:

Feb 16, 2016, 6:37:21 PM (9 years ago)

Author:

jyg

Message:

Correction in the computation of qsat in
physiq_mod.F90: the transition of humidity from
saturation over ice to saturation over liquid
water occurs now at rtt (=273.15 K) instead of
t_glace_min or t_coup.

File:

• LMDZ5/trunk/libf/phylmd/physiq_mod.F90 (modified) (3 diffs)

LMDZ5/trunk/libf/phylmd/physiq_mod.F90

@@ -2103 +2103 @@
            zx_qs  = zx_qs*zcor
         ELSE
-           IF (zx_t.LT.t_coup) THEN
+!!           IF (zx_t.LT.t_coup) THEN             !jyg
+           IF (zx_t.LT.rtt) THEN                  !jyg
               zx_qs = qsats(zx_t)/pplay(i,k)
            ELSE
@@ -3371 +3372 @@
         zx_t = t_seri(i,k)
         IF (thermcep) THEN
-           if (iflag_ice_thermo.eq.0) then 
+!!           if (iflag_ice_thermo.eq.0) then                 !jyg
            zdelta = MAX(0.,SIGN(1.,rtt-zx_t))
-           else
-           zdelta = MAX(0.,SIGN(1.,t_glace_min-zx_t))
-           endif
+!!           else                                            !jyg
+!!           zdelta = MAX(0.,SIGN(1.,t_glace_min-zx_t))      !jyg
+!!           endif                                           !jyg
            zx_qs  = r2es * FOEEW(zx_t,zdelta)/pplay(i,k)
            zx_qs  = MIN(0.5,zx_qs)
@@ -3381 +3382 @@
            zx_qs  = zx_qs*zcor
         ELSE
-           IF (zx_t.LT.t_coup) THEN
+!!           IF (zx_t.LT.t_coup) THEN             !jyg
+           IF (zx_t.LT.rtt) THEN                  !jyg
               zx_qs = qsats(zx_t)/pplay(i,k)
            ELSE