Ignore:
Timestamp:
Dec 14, 2019, 10:15:37 AM (5 years ago)
Author:
jyg
Message:

Bug fix:
Missing a max( . . . ,0.) in the computation of
the conversion of cloud water in precipitation
(subroutine cv3_undilute2 in cv3_routines.F90).

File:
1 edited

Legend:

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

    r3502 r3614  
    14951495!!!!            qcld(i,k) = min(clw(i,k), elacrit)                                          !!jygprl
    14961496            qcld(i,k) = min(clw(i,k), elacrit*(1.-qta(i,k-1))/(1.-elacrit))             !!jygprl
    1497             phinu2p = qhsat(i,k-1) + qcld(i,k-1) - (qhsat(i,k) + qcld(i,k))                     !!jygprl
     1497            phinu2p = max(qhsat(i,k-1) + qcld(i,k-1) - (qhsat(i,k) + qcld(i,k)),0.)   !!jygprl
    14981498            qpl(i,k) = qpl(i,k-1) + (1.-frac(i,k))*phinu2p                            !!jygprl
    14991499            qps(i,k) = qps(i,k-1) + frac(i,k)     *phinu2p                            !!jygprl
Note: See TracChangeset for help on using the changeset viewer.