/[lmdze]/trunk/phylmd/Interface_surf/fonte_neige.f
ViewVC logotype

Diff of /trunk/phylmd/Interface_surf/fonte_neige.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 102 by guez, Mon Jul 7 17:45:21 2014 UTC revision 103 by guez, Fri Aug 29 13:00:05 2014 UTC
# Line 78  contains Line 78  contains
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
# Line 96  contains Line 97  contains
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)
Legend:
Removed from v.102  
changed lines
  Added in v.103
  ViewVC Help
Powered by ViewVC 1.1.21