75 |
real, dimension(klon) :: zx_pkh, zx_dq_s_dt, zx_qsat, zx_coef |
real, dimension(klon) :: zx_pkh, zx_dq_s_dt, zx_qsat, zx_coef |
76 |
real, dimension(klon) :: zx_sl, zx_k1 |
real, dimension(klon) :: zx_sl, zx_k1 |
77 |
real, dimension(klon) :: zx_q_0 , d_ts |
real, dimension(klon) :: zx_q_0 , d_ts |
78 |
real :: zdelta, zcvm5, zx_qs, zcor, zx_dq_s_dh |
logical zdelta |
79 |
|
real zcvm5, zx_qs, zcor, zx_dq_s_dh |
80 |
real :: bilan_f, fq_fonte |
real :: bilan_f, fq_fonte |
81 |
REAL :: subli, fsno |
REAL :: subli, fsno |
82 |
REAL :: qsat_new, q1_new |
REAL :: qsat_new, q1_new |
121 |
DO i = 1, knon |
DO i = 1, knon |
122 |
zx_pkh(i) = (ps(i)/ps(i))**RKAPPA |
zx_pkh(i) = (ps(i)/ps(i))**RKAPPA |
123 |
IF (thermcep) THEN |
IF (thermcep) THEN |
124 |
zdelta=MAX(0., SIGN(1., rtt-tsurf(i))) |
zdelta= rtt >= tsurf(i) |
125 |
zcvm5 = R5LES*RLVTT*(1.-zdelta) + R5IES*RLSTT*zdelta |
zcvm5 = merge(R5IES*RLSTT, R5LES*RLVTT, zdelta) |
126 |
zcvm5 = zcvm5 / RCPD / (1.0+RVTMP2*q1lay(i)) |
zcvm5 = zcvm5 / RCPD / (1.0+RVTMP2*q1lay(i)) |
127 |
zx_qs= r2es * FOEEW(tsurf(i), zdelta)/ps(i) |
zx_qs= r2es * FOEEW(tsurf(i), zdelta)/ps(i) |
128 |
zx_qs=MIN(0.5, zx_qs) |
zx_qs=MIN(0.5, zx_qs) |