--- trunk/phylmd/Conflx/flxflux.f 2014/07/15 13:43:24 102 +++ trunk/phylmd/Conflx/flxflux.f 2014/08/29 13:00:05 103 @@ -49,7 +49,7 @@ ! Local: REAL cevapcu(klev) - REAL ztmsmlt, zdelta, zqsat + REAL ztmsmlt, zqsat !jq The variable maxpdmfdp(klon) has been introduced by Olivier Boucher !jq 14/11/00 to fix the problem with the negative precipitation. @@ -160,8 +160,7 @@ zsnmlt=MIN(pmflxs(i, k), zfac*(pten(i, k)-ztmelp2)) pdpmel(i, k)=zsnmlt ztmsmlt=pten(i, k)-zsnmlt/zfac - zdelta=MAX(0., SIGN(1., RTT-ztmsmlt)) - zqsat = R2ES * FOEEW(ztmsmlt, zdelta) / pap(i, k) + zqsat = R2ES * FOEEW(ztmsmlt, RTT >= ztmsmlt) / pap(i, k) zqsat = MIN(0.5, zqsat) zqsat = zqsat / (1. - RETV * zqsat) pqsen(i, k) = zqsat