--- trunk/libf/phylmd/hbtm.f90	2011/08/24 11:43:14	49
+++ trunk/phylmd/hbtm.f	2014/08/29 13:00:05	103
@@ -29,7 +29,7 @@
     ! mais on garde la possibilité de changer si besoin est (jusqu'à présent
     ! la forme de HB avec le 1er niveau modele etait conservee)
 
-    USE dimphy, ONLY: klev, klon, max
+    USE dimphy, ONLY: klev, klon
     USE suphec_m, ONLY: rcpd, rd, retv, rg, rkappa, rlvtt, rtt, rv
     USE yoethf_m, ONLY: r2es, rvtmp2
     USE fcttre, ONLY: foeew
@@ -198,7 +198,7 @@
     REAL zm(klon)
     ! current level height + one level up
     REAL zp(klon)
-    REAL zcor, zdelta, zcvm5
+    REAL zcor, zcvm5
 
     REAL fac, pblmin, zmzp, term
 
@@ -520,8 +520,7 @@
              if (.not. Zsat(i)) then
                 T2 = T2m(i) * s(i, k)
                 ! thermodyn functions
-                zdelta=MAX(0., SIGN(1., RTT - T2))
-                qqsat= r2es * FOEEW(T2, zdelta) / pplay(i, k)
+                qqsat= r2es * FOEEW(T2, RTT >= T2) / pplay(i, k)
                 qqsat=MIN(0.5, qqsat)
                 zcor=1./(1.-retv*qqsat)
                 qqsat=qqsat*zcor