87 |
PARAMETER (t_coup=234.0) |
PARAMETER (t_coup=234.0) |
88 |
|
|
89 |
INTEGER i, k, n, kk |
INTEGER i, k, n, kk |
90 |
REAL zqs(klon), zdqs(klon), zdelta, zcor, zcvm5 |
REAL zqs(klon), zdqs(klon), zcor, zcvm5 |
91 |
|
logical zdelta |
92 |
REAL zrfl(klon), zrfln(klon), zqev, zqevt |
REAL zrfl(klon), zrfln(klon), zqev, zqevt |
93 |
REAL zoliq(klon), zcond(klon), zq(klon), zqn(klon), zdelq |
REAL zoliq(klon), zcond(klon), zq(klon), zqn(klon), zdelq |
94 |
REAL ztglace, zt(klon) |
REAL ztglace, zt(klon) |
226 |
DO i = 1, klon |
DO i = 1, klon |
227 |
IF (zrfl(i)>0.) THEN |
IF (zrfl(i)>0.) THEN |
228 |
IF (thermcep) THEN |
IF (thermcep) THEN |
229 |
zdelta = max(0., sign(1., rtt-zt(i))) |
zqs(i) = r2es*foeew(zt(i), rtt >= zt(i))/pplay(i, k) |
|
zqs(i) = r2es*foeew(zt(i), zdelta)/pplay(i, k) |
|
230 |
zqs(i) = min(0.5, zqs(i)) |
zqs(i) = min(0.5, zqs(i)) |
231 |
zcor = 1./(1.-retv*zqs(i)) |
zcor = 1./(1.-retv*zqs(i)) |
232 |
zqs(i) = zqs(i)*zcor |
zqs(i) = zqs(i)*zcor |
264 |
|
|
265 |
IF (thermcep) THEN |
IF (thermcep) THEN |
266 |
DO i = 1, klon |
DO i = 1, klon |
267 |
zdelta = max(0., sign(1., rtt-zt(i))) |
zdelta = rtt >= zt(i) |
268 |
zcvm5 = r5les*rlvtt*(1.-zdelta) + r5ies*rlstt*zdelta |
zcvm5 = merge(r5ies*rlstt, r5les*rlvtt, zdelta) |
269 |
zcvm5 = zcvm5/rcpd/(1.0+rvtmp2*zq(i)) |
zcvm5 = zcvm5/rcpd/(1.0+rvtmp2*zq(i)) |
270 |
zqs(i) = r2es*foeew(zt(i), zdelta)/pplay(i, k) |
zqs(i) = r2es*foeew(zt(i), zdelta)/pplay(i, k) |
271 |
zqs(i) = min(0.5, zqs(i)) |
zqs(i) = min(0.5, zqs(i)) |