Changeset 3624


Ignore:
Timestamp:
Jan 27, 2020, 10:50:54 AM (4 years ago)
Author:
jyg
Message:

Correction to commit 3614:
A new formula is used for "phinu2p" (hopefully
correct) in subroutine cv3_undilute2 in
cv3_routines.F90.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/trunk/libf/phylmd/cv3_routines.F90

    r3614 r3624  
    11901190  REAL tg, qg, dqgdT, ahg, alv, alf, s, tc, es, esi, denom, rg, tca, elacrit
    11911191  REAL                                               :: phinu2p
    1192   REAL als
     1192  REAL                                               :: qhthreshold
     1193  REAL                                               :: als
    11931194  REAL                                               :: qsat_new, snew
    11941195  REAL, DIMENSION (nloc,nd)                          :: qi
     
    14941495            elacrit = elcrit*min(max(1.-(tp(i,k)-T0)/Tlcrit, 0.), 1.)                   !!jygprl
    14951496!!!!            qcld(i,k) = min(clw(i,k), elacrit)                                          !!jygprl
    1496             qcld(i,k) = min(clw(i,k), elacrit*(1.-qta(i,k-1))/(1.-elacrit))             !!jygprl
    1497             phinu2p = max(qhsat(i,k-1) + qcld(i,k-1) - (qhsat(i,k) + qcld(i,k)),0.)   !!jygprl
     1497            qhthreshold = elacrit*(1.-qta(i,k-1))/(1.-elacrit)
     1498            qcld(i,k) = min(clw(i,k), qhthreshold)             !!jygprl
     1499!!!!            phinu2p = max(qhsat(i,k-1) + qcld(i,k-1) - (qhsat(i,k) + qcld(i,k)),0.)   !!jygprl
     1500            phinu2p = max(clw(i,k) - max(qta(i,k-1) - qhsat(i,k-1), qhthreshold), 0.)
    14981501            qpl(i,k) = qpl(i,k-1) + (1.-frac(i,k))*phinu2p                            !!jygprl
    14991502            qps(i,k) = qps(i,k-1) + frac(i,k)     *phinu2p                            !!jygprl
Note: See TracChangeset for help on using the changeset viewer.