Changeset 689


Ignore:
Timestamp:
Jun 4, 2012, 10:12:14 AM (13 years ago)
Author:
flefevre
Message:

Minor change. Better handling of very small water mixing ratio at saturation
(for very low temperatures).

Location:
trunk/LMDZ.MARS/libf/phymars
Files:
2 edited

Legend:

Unmodified
Added
Removed

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

    r663 r689  
    260260c       Get the partial pressure of water vapor and its saturation ratio
    261261        ph2o = zq(ig,l,igcm_h2o_vap) * (mmean(ig,l)/18.) * pplay(ig,l)
    262         satu = zq(ig,l,igcm_h2o_vap) / max(zqsat(ig,l), 1.e-30)
     262        satu = zq(ig,l,igcm_h2o_vap) / zqsat(ig,l)
    263263
    264264!=============================================================

  • trunk/LMDZ.MARS/libf/phymars/watersat.F

    r38 r689  
    3636        es=10.**es
    3737        qsat(i)=(es*100.)/p(i)*0.41
     38        qsat(i) = max(qsat(i), 1.e-30)
    3839
    3940c        qss = specific humidity
Note: See TracChangeset for help on using the changeset viewer.