Changeset 663 for trunk

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

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/phymars/improvedclouds.F (modified) (2 diffs)

trunk/LMDZ.MARS/README ¶

@@ -1679 +1679 @@
   local array).
 
-== 21/05/12 == JYC
+== 21/05/12 == JYC+FL
 - corrected bug in computation of HCO2+ in paramfoto_compact.F
+- avoid infinite H2O saturation ratio at very low temperatures in improvedclouds.F

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)
 
 !=============================================================