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

Diff of /trunk/phylmd/hbtm.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 198  contains Line 198  contains
198      REAL zm(klon)      REAL zm(klon)
199      ! current level height + one level up      ! current level height + one level up
200      REAL zp(klon)      REAL zp(klon)
201      REAL zcor, zdelta, zcvm5      REAL zcor, zcvm5
202    
203      REAL fac, pblmin, zmzp, term      REAL fac, pblmin, zmzp, term
204    
# Line 520  contains Line 520  contains
520               if (.not. Zsat(i)) then               if (.not. Zsat(i)) then
521                  T2 = T2m(i) * s(i, k)                  T2 = T2m(i) * s(i, k)
522                  ! thermodyn functions                  ! thermodyn functions
523                  zdelta=MAX(0., SIGN(1., RTT - T2))                  qqsat= r2es * FOEEW(T2, RTT >= T2) / pplay(i, k)
                 qqsat= r2es * FOEEW(T2, zdelta) / pplay(i, k)  
524                  qqsat=MIN(0.5, qqsat)                  qqsat=MIN(0.5, qqsat)
525                  zcor=1./(1.-retv*qqsat)                  zcor=1./(1.-retv*qqsat)
526                  qqsat=qqsat*zcor                  qqsat=qqsat*zcor
Legend:
Removed from v.102  
changed lines
  Added in v.103
  ViewVC Help
Powered by ViewVC 1.1.21