/[lmdze]/trunk/phylmd/Interface_surf/hbtm.f90
ViewVC logotype

Diff of /trunk/phylmd/Interface_surf/hbtm.f90

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

trunk/libf/phylmd/hbtm.f90 revision 49 by guez, Wed Aug 24 11:43:14 2011 UTC trunk/phylmd/hbtm.f revision 103 by guez, Fri Aug 29 13:00:05 2014 UTC
# Line 29  contains Line 29  contains
29      ! mais on garde la possibilité de changer si besoin est (jusqu'à présent      ! mais on garde la possibilité de changer si besoin est (jusqu'à présent
30      ! la forme de HB avec le 1er niveau modele etait conservee)      ! la forme de HB avec le 1er niveau modele etait conservee)
31    
32      USE dimphy, ONLY: klev, klon, max      USE dimphy, ONLY: klev, klon
33      USE suphec_m, ONLY: rcpd, rd, retv, rg, rkappa, rlvtt, rtt, rv      USE suphec_m, ONLY: rcpd, rd, retv, rg, rkappa, rlvtt, rtt, rv
34      USE yoethf_m, ONLY: r2es, rvtmp2      USE yoethf_m, ONLY: r2es, rvtmp2
35      USE fcttre, ONLY: foeew      USE fcttre, ONLY: foeew
# 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.49  
changed lines
  Added in v.103
  ViewVC Help
Powered by ViewVC 1.1.21