--- trunk/Sources/phylmd/suphec.f	2015/09/14 17:13:16	169
+++ trunk/Sources/phylmd/suphec.f	2016/06/06 17:42:15	201
@@ -11,10 +11,14 @@
   real, parameter:: NAVO = 6.0221367E23 ! Avogadro number, in mol-1
 
   ! A1.1 Astronomical constants
-  REAL RSIYEA, RSIDAY, ROMEGA
+
+  REAL ROMEGA
   real, parameter:: RDAY = 86400.
-  real, parameter:: REA = 149597870000.
-  real, parameter:: REPSM = 0.409093
+
+  REAL, parameter:: RSIYEA = 365.25 * RDAY * 2. * PI / 6.283076
+  ! sideral year, in s
+
+  REAL, parameter:: RSIDAY = RDAY / (1. + RDAY / RSIYEA) ! sideral day, in s
 
   ! A1.2 Geoide
   real, parameter:: RG = 9.80665 ! acceleration of gravity, in m s-2
@@ -38,22 +42,31 @@
   ! specific ideal gas constant for dry air, in J K-1 kg-1
   ! (factor 1e3: conversion from g to kg)
 
-  real, save:: RCPV, RCVD, RCVV
+  real, parameter:: RCPV = 4. * RV 
+  ! specific heat capacity at constant pressure of water vapor, in J K-1 kg-1
+
+  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
   real, save:: RETV
 
   ! A1.5, 6 Thermodynamic liquid, solid phases
-  REAL, save:: RCW, RCS
+
+  REAL, parameter:: RCW = RCPV ! LIQUID PHASE Cw
+  real, save:: RCS
 
   ! 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
 
@@ -77,16 +90,9 @@
 
     ! 2. DEFINE ASTRONOMICAL CONSTANTS
 
-    RSIYEA = 365.25*RDAY*2.*PI/6.283076
-    RSIDAY = RDAY/(1.+RDAY/RSIYEA)
     ROMEGA = 2.*PI/RSIDAY
 
     print *, 'Astronomical constants '
-    print '('' day = '', E13.7, '' s'')', RDAY
-    print '('' half g. axis = '', E13.7, '' m'')', REA
-    print '('' mean anomaly = '', E13.7, '' -'')', REPSM
-    print '('' sideral year = '', E13.7, '' s'')', RSIYEA
-    print '('' sideral day = '', E13.7, '' s'')', RSIDAY
     print '('' omega = '', E13.7, '' s-1'')', ROMEGA
 
     ! 3. DEFINE GEOIDE.
@@ -103,7 +109,6 @@
     ! 5. DEFINE THERMODYNAMIC CONSTANTS, GAS PHASE.
 
     RCVD = RCPD-RD
-    RCPV = 4. * RV
     RCVV = RCPV-RV
     RETV = RV / RD - 1.
     print *, 'Thermodynamics, gas'
@@ -112,17 +117,10 @@
     print *, "rv = ", RV, "J K-1 kg-1"
     print '('' Cpd = '', e13.7)', RCPD
     print '('' Cvd = '', e13.7)', RCVD
-    print '('' Cpv = '', e13.7)', RCPV
     print '('' Cvv = '', e13.7)', RCVV
     print '('' Rd/Cpd = '', e13.7)', RKAPPA
     print '('' Rv / Rd - 1 = '', e13.7)', RETV
 
-    ! 6. DEFINE THERMODYNAMIC CONSTANTS, LIQUID PHASE.
-
-    RCW = RCPV
-    print *, 'Thermodynamic, liquid '
-    print '('' Cw = '', E13.7)', RCW
-
     ! 7. DEFINE THERMODYNAMIC CONSTANTS, SOLID PHASE.
 
     RCS = RCPV
@@ -134,7 +132,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