Changeset 2448


Ignore:
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:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ5/trunk/libf/phylmd/physiq_mod.F90

    r2443 r2448  
    21032103           zx_qs  = zx_qs*zcor
    21042104        ELSE
    2105            IF (zx_t.LT.t_coup) THEN
     2105!!           IF (zx_t.LT.t_coup) THEN             !jyg
     2106           IF (zx_t.LT.rtt) THEN                  !jyg
    21062107              zx_qs = qsats(zx_t)/pplay(i,k)
    21072108           ELSE
     
    33713372        zx_t = t_seri(i,k)
    33723373        IF (thermcep) THEN
    3373            if (iflag_ice_thermo.eq.0) then
     3374!!           if (iflag_ice_thermo.eq.0) then                 !jyg
    33743375           zdelta = MAX(0.,SIGN(1.,rtt-zx_t))
    3375            else
    3376            zdelta = MAX(0.,SIGN(1.,t_glace_min-zx_t))
    3377            endif
     3376!!           else                                            !jyg
     3377!!           zdelta = MAX(0.,SIGN(1.,t_glace_min-zx_t))      !jyg
     3378!!           endif                                           !jyg
    33783379           zx_qs  = r2es * FOEEW(zx_t,zdelta)/pplay(i,k)
    33793380           zx_qs  = MIN(0.5,zx_qs)
     
    33813382           zx_qs  = zx_qs*zcor
    33823383        ELSE
    3383            IF (zx_t.LT.t_coup) THEN
     3384!!           IF (zx_t.LT.t_coup) THEN             !jyg
     3385           IF (zx_t.LT.rtt) THEN                  !jyg
    33843386              zx_qs = qsats(zx_t)/pplay(i,k)
    33853387           ELSE
Note: See TracChangeset for help on using the changeset viewer.