--- trunk/libf/phylmd/YOMCST.f90	2011/01/06 17:52:19	38
+++ trunk/phylmd/YOMCST.f	2015/02/06 15:00:28	125
@@ -4,11 +4,27 @@
 
   ! From phylmd/YOMCST.h, version 1.2 2005/06/06 13:16:33
 
-  ! A1.1.bis Constantes concernant l'orbite de la Terre:
-  REAL R_ecc ! excentricité
-  real R_peri ! équinoxe
-  real R_incl ! inclinaison
+  ! 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
+
+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