49 |
|
|
50 |
! Local: |
! Local: |
51 |
REAL cevapcu(klev) |
REAL cevapcu(klev) |
52 |
REAL ztmsmlt, zdelta, zqsat |
REAL ztmsmlt, zqsat |
53 |
|
|
54 |
!jq The variable maxpdmfdp(klon) has been introduced by Olivier Boucher |
!jq The variable maxpdmfdp(klon) has been introduced by Olivier Boucher |
55 |
!jq 14/11/00 to fix the problem with the negative precipitation. |
!jq 14/11/00 to fix the problem with the negative precipitation. |
160 |
zsnmlt=MIN(pmflxs(i, k), zfac*(pten(i, k)-ztmelp2)) |
zsnmlt=MIN(pmflxs(i, k), zfac*(pten(i, k)-ztmelp2)) |
161 |
pdpmel(i, k)=zsnmlt |
pdpmel(i, k)=zsnmlt |
162 |
ztmsmlt=pten(i, k)-zsnmlt/zfac |
ztmsmlt=pten(i, k)-zsnmlt/zfac |
163 |
zdelta=MAX(0., SIGN(1., RTT-ztmsmlt)) |
zqsat = R2ES * FOEEW(ztmsmlt, RTT >= ztmsmlt) / pap(i, k) |
|
zqsat = R2ES * FOEEW(ztmsmlt, zdelta) / pap(i, k) |
|
164 |
zqsat = MIN(0.5, zqsat) |
zqsat = MIN(0.5, zqsat) |
165 |
zqsat = zqsat / (1. - RETV * zqsat) |
zqsat = zqsat / (1. - RETV * zqsat) |
166 |
pqsen(i, k) = zqsat |
pqsen(i, k) = zqsat |