--- trunk/libf/phylmd/FCTTRE.f90	2010/12/21 15:45:48	37
+++ trunk/phylmd/FCTTRE.f90	2013/11/15 18:45:49	76
@@ -2,11 +2,14 @@
 
   ! From phylmd/FCTTRE.inc, version 1.2 2004/06/22 11:45:32
 
-  ! This COMDECK includes the thermodynamical functions for the cycle
-  ! 39 ECMWF Physics package. Consistent with YOMCST basic physics
-  ! constants, assuming the partial pressure of water vapour is given
-  ! by a first order Taylor expansion of Qs(T) with respect to
-  ! temperature, using constants in yoethf.
+  ! This module includes the thermodynamical functions for the cycle
+  ! 39 ECMWF Physics package. Consistent with "SUPHEC_M" basic
+  ! physical constants, assuming the partial pressure of water vapour
+  ! is given by a first order Taylor expansion of "Qs(T)" with respect
+  ! to temperature, using constants in "yoethf_m".
+
+  ! Probably from Buck, 1981, Journal of Applied Meteorology, volume
+  ! 20, number 12, page 1527.
 
   implicit none
 
@@ -16,10 +19,11 @@
 
   REAL function FOEEW(T, DEL)
 
-    use yoethf, only: R3LES, R3IES, R4LES, R4IES
-    use YOMCST, only: rtt
+    use yoethf_m, only: R3LES, R3IES, R4LES, R4IES
+    use SUPHEC_M, only: rtt
 
-    REAL, intent(in):: T, DEL
+    REAL, intent(in):: T
+    REAL, intent(in):: DEL ! 1 for ice, 0 for liquid
 
     !-----------------------
 
@@ -32,7 +36,7 @@
 
   REAL function FOEDE(T, DEL, P5ARG, QS, PCOARG)
 
-    use yoethf, only: R4LES, R4IES
+    use yoethf_m, only: R4LES, R4IES
 
     REAL, intent(in):: T, DEL
     real, intent(in):: P5ARG, QS, PCOARG
@@ -74,7 +78,7 @@
 
   REAL function dqsats(t, qs)
 
-    use YOMCST, only: RLVTT, rcpd
+    use SUPHEC_M, only: RLVTT, rcpd
 
     REAL, intent(in):: T, qs
 
@@ -89,7 +93,7 @@
 
   REAL function dqsatl(t, qs)
 
-    use YOMCST, only: RLVTT, rcpd
+    use SUPHEC_M, only: RLVTT, rcpd
 
     REAL, intent(in):: T, qs