--- trunk/libf/phylmd/YOMCST.f90	2008/02/27 13:16:39	3
+++ trunk/phylmd/YOMCST.f	2018/12/11 14:12:45	317
@@ -4,32 +4,27 @@
 
   ! From phylmd/YOMCST.h, version 1.2 2005/06/06 13:16:33
 
-  ! A1.0 Fundamental constants
-  REAL RPI,RCLUM,RHPLA,RKBOL,RNAVO
-  ! A1.1 Astronomical constants
-  REAL RDAY,REA,REPSM,RSIYEA,RSIDAY,ROMEGA
-
-  ! A1.1.bis Constantes concernant l'orbite de la Terre:
-  REAL R_ecc ! excentricite
-  real R_peri ! equinoxe
-  real R_incl ! inclinaison
-
-  ! A1.2 Geoide
-  REAL RA,RG,R1SA
-  ! A1.3 Radiation
-  !     REAL RSIGMA,RI0
-  REAL RSIGMA
-  ! A1.4 Thermodynamic gas phase
-  REAL R,RMD,RMO3,RMV,RD,RV,RCPD,RCPV,RCVD,RCVV
-  REAL RKAPPA,RETV
-  ! A1.5,6 Thermodynamic liquid,solid phases
-  REAL RCW,RCS
-  ! A1.7 Thermodynamic transition of phase
-  REAL RLVTT,RLSTT,RLMLT,RTT,RATM
-  ! A1.8 Curve of saturation
-  REAL RESTT,RALPW,RBETW,RGAMW,RALPS,RBETS,RGAMS
-  REAL RALPD,RBETD,RGAMD
+  ! A1.1.bis Constantes concernant l'orbite de la Terre
+  ! Default values from AMIP II
 
-  save
+  REAL:: R_ecc = 0.016715 ! excentricité
+  real:: R_peri = 102.7 ! équinoxe
+  real:: R_incl = 23.441 ! inclinaison, en degr\'es
+
+contains
+
+  subroutine read_YOMCST
+
+    use unit_nml_m, only: unit_nml
+
+    namelist /YOMCST_nml/ R_ecc, R_peri, R_incl
+
+    !------------------------------------
+
+    print *, "Enter namelist 'YOMCST_nml'."
+    read(unit=*, nml=YOMCST_nml)
+    write(unit_nml, nml=YOMCST_nml)
+
+  end subroutine read_YOMCST
 
 end module YOMCST