78 |
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 |
79 |
real, dimension(klon):: zx_sl, zx_k1 |
real, dimension(klon):: zx_sl, zx_k1 |
80 |
real, dimension(klon):: d_ts |
real, dimension(klon):: d_ts |
81 |
real zdelta, zcvm5, zx_qs, zcor, zx_dq_s_dh |
logical zdelta |
82 |
|
real zcvm5, zx_qs, zcor, zx_dq_s_dh |
83 |
real fq_fonte |
real fq_fonte |
84 |
REAL bil_eau_s(knon) ! in kg m-2 |
REAL bil_eau_s(knon) ! in kg m-2 |
85 |
real snow_evap(klon) ! in kg m-2 s-1 |
real snow_evap(klon) ! in kg m-2 s-1 |
97 |
DO i = 1, knon |
DO i = 1, knon |
98 |
zx_pkh(i) = (ps(i)/ps(i))**RKAPPA |
zx_pkh(i) = (ps(i)/ps(i))**RKAPPA |
99 |
IF (thermcep) THEN |
IF (thermcep) THEN |
100 |
zdelta=MAX(0., SIGN(1., rtt-tsurf(i))) |
zdelta= rtt >= tsurf(i) |
101 |
zcvm5 = R5LES*RLVTT*(1.-zdelta) + R5IES*RLSTT*zdelta |
zcvm5 = merge(R5IES*RLSTT, R5LES*RLVTT, zdelta) |
102 |
zcvm5 = zcvm5 / RCPD / (1. + RVTMP2*q1lay(i)) |
zcvm5 = zcvm5 / RCPD / (1. + RVTMP2*q1lay(i)) |
103 |
zx_qs= r2es * FOEEW(tsurf(i), zdelta)/ps(i) |
zx_qs= r2es * FOEEW(tsurf(i), zdelta)/ps(i) |
104 |
zx_qs=MIN(0.5, zx_qs) |
zx_qs=MIN(0.5, zx_qs) |