/[lmdze]/trunk/phylmd/fisrtilp.f
ViewVC logotype

Diff of /trunk/phylmd/fisrtilp.f

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

revision 102 by guez, Wed Mar 5 14:57:53 2014 UTC revision 103 by guez, Fri Aug 29 13:00:05 2014 UTC
# Line 87  contains Line 87  contains
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)
# Line 225  contains Line 226  contains
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
# Line 264  contains Line 264  contains
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))
Legend:
Removed from v.102  
changed lines
  Added in v.103
  ViewVC Help
Powered by ViewVC 1.1.21