--- trunk/phylmd/hbtm.f	2014/07/15 13:43:24	102
+++ trunk/phylmd/hbtm.f	2014/08/29 13:00:05	103
@@ -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