Changeset 663 for trunk/LMDZ.MARS/libf

Timestamp:

May 21, 2012, 1:35:29 PM (13 years ago)

Author:

flefevre

Message:

Improvement in the stability of the water cycle scheme at very low
temperatures: avoid infinite values of the H2O saturation ratio
in improvedclouds.F

File:

• trunk/LMDZ.MARS/libf/phymars/improvedclouds.F (modified) (2 diffs)

trunk/LMDZ.MARS/libf/phymars/improvedclouds.F

@@ -36 +36 @@
 #include "dimradmars.h"
 #include "microphys.h"
+#include "conc.h"
 c------------------------------------------------------------------
 c     Inputs:
@@ -258 +259 @@
 
 c       Get the partial pressure of water vapor and its saturation ratio
-        ph2o = zq(ig,l,igcm_h2o_vap) * (44./18.) * pplay(ig,l)
-        satu = zq(ig,l,igcm_h2o_vap) / zqsat(ig,l)
+        ph2o = zq(ig,l,igcm_h2o_vap) * (mmean(ig,l)/18.) * pplay(ig,l)
+        satu = zq(ig,l,igcm_h2o_vap) / max(zqsat(ig,l), 1.e-30)
 
 !=============================================================