34 |
! Variables locales: |
! Variables locales: |
35 |
INTEGER i, k, kb, invb(klon) |
INTEGER i, k, kb, invb(klon) |
36 |
REAL zqs, zrhb, zcll, zdthmin(klon), zdthdp |
REAL zqs, zrhb, zcll, zdthmin(klon), zdthdp |
37 |
REAL zdelta, zcor |
REAL zcor |
38 |
|
|
39 |
!----------------------------------------------------------- |
!----------------------------------------------------------- |
40 |
|
|
67 |
DO i = 1, klon |
DO i = 1, klon |
68 |
kb=invb(i) |
kb=invb(i) |
69 |
IF (thermcep) THEN |
IF (thermcep) THEN |
70 |
zdelta=MAX(0., SIGN(1., RTT-t(i, kb))) |
zqs= R2ES*FOEEW(t(i, kb), RTT >= t(i, kb))/pplay(i, kb) |
|
zqs= R2ES*FOEEW(t(i, kb), zdelta)/pplay(i, kb) |
|
71 |
zqs=MIN(0.5, zqs) |
zqs=MIN(0.5, zqs) |
72 |
zcor=1./(1.-RETV*zqs) |
zcor=1./(1.-RETV*zqs) |
73 |
zqs=zqs*zcor |
zqs=zqs*zcor |