--- trunk/Sources/phylmd/suphec.f	2016/05/31 16:17:35	198
+++ trunk/phylmd/suphec.f	2018/02/05 10:39:38	254
@@ -42,12 +42,14 @@
   ! specific ideal gas constant for dry air, in J K-1 kg-1
   ! (factor 1e3: conversion from g to kg)
 
-  real, parameter:: RCPV = 4. * RV ! Cpv, gas phase
+  real, parameter:: RCPV = 4. * RV 
+  ! specific heat capacity at constant pressure of water vapor, in J K-1 kg-1
 
-  real, save:: RCVD, RCVV
+  real, save:: RCVD
+  real, save:: RCVV
 
   real, parameter:: RCPD = 7. / 2 * RD 
-  ! specific heat capacity for dry air, in J K-1 kg-1
+  ! specific heat capacity at constant pressure of dry air, in J K-1 kg-1
 
   real, parameter:: RMO3 = 47.9942
   REAL, parameter:: RKAPPA = RD/RCPD
@@ -61,7 +63,10 @@
   ! A1.7 Thermodynamic transition of phase
   REAL, save:: RLMLT
   real, parameter:: RTT = 273.16
-  real, parameter:: RLVTT = 2.5008E+6
+
+  real, parameter:: RLVTT = 2.5008E+6 
+  ! specific latent heat of vaporization of water at triple point, in J kg-1
+
   real, parameter:: RLSTT = 2.8345E+6
   real, parameter:: RATM = 1e5
 
@@ -107,7 +112,6 @@
     RCVV = RCPV-RV
     RETV = RV / RD - 1.
     print *, 'Thermodynamics, gas'
-    print '('' Ozone mass = '', e13.7)', RMO3
     print *, "rd = ", RD, "J K-1 kg-1"
     print *, "rv = ", RV, "J K-1 kg-1"
     print '('' Cpd = '', e13.7)', RCPD
@@ -127,7 +131,6 @@
     RLMLT = RLSTT-RLVTT
     print *, 'Thermodynamic, transition of phase:'
     print '('' Fusion point = '', E13.7)', RTT
-    print '('' RLvTt = '', E13.7)', RLVTT
     print '('' RLsTt = '', E13.7)', RLSTT
     print '('' RLMlt = '', E13.7)', RLMLT
     print '('' Normal pressure = '', E13.7)', RATM