--- trunk/libf/dyn3d/conf_gcm.f90	2011/02/22 13:49:36	40
+++ trunk/dyn3d/conf_gcm.f	2015/02/27 16:44:07	131
@@ -20,21 +20,29 @@
   integer:: iphysiq = 5
   ! number of time steps of dynamics between calls to physics
 
-  integer:: dayref = 1 ! jour de l'année de l'état initial
-  ! (= 350 si 20 décembre par exemple)
-
-  integer:: anneeref = 1998 ! Annee de l'etat initial (avec 4 chiffres)
-
   logical:: raz_date = .false.
-  ! (remise a zero de la date initiale, prise en compte de la date de
-  ! gcm.def avec remise a zero des compteurs de pas de temps)
-  ! (pas de remise a zero: on garde la date du fichier restart)
+  ! prise en compte de la date initiale de la namelist et remise à
+  ! zéro des compteurs de pas de temps (sinon on garde la date du
+  ! fichier restart)
 
-  REAL:: periodav = 1. ! time interval between outputs to "histmoy" (in days) 
+  integer:: periodav = 1 
+  ! time interval between outputs in the dynamical part, in days
 
   logical:: offline = .FALSE.
   ! permet de mettre en route la nouvelle parametrisation de l'eau liquide
 
+  integer:: prt_level = 0 ! niveau d'impression souhaité (0 = minimum)
+
+  LOGICAL:: purmats= .FALSE.
+  ! Help = Choix du schema d'integration temporel.
+  ! y = pure Matsuno sinon c'est du Matsuno-leapfrog
+
+  INTEGER:: iflag_phys = 1
+  ! contrôle l'appel à la physique :
+  ! 0 : pas de physique
+  ! 1 : physique normale (appel à phylmd, phymars...) (default)
+  ! 2 : rappel Newtonien pour la température + friction au sol
+
 contains
 
   SUBROUTINE conf_gcm
@@ -42,61 +50,45 @@
     ! Auteurs : L. Fairhead, P. Le Van
     ! Version du 29/04/97
 
-    ! On ne compare pas les valeurs des paramètres du zoom, grossismx, 
-    ! grossismy, clon, clat, fxyhypb lues sur le fichier start avec
-    ! celles passées par run.def, au début du gcm. 
-    ! Ces paramètres définissent entre autres la grille et doivent être
-    ! cohérents, sinon il y aura divergence du gcm.
-
     use abort_gcm_m, only: abort_gcm
     use comdissnew, only: read_comdissnew
-    use logic, only: read_logic
-    use serre, only: clon, clat, grossismx, grossismy, alphax, alphay, &
-         dzoomx, dzoomy, taux, tauy
-    use iniprint, only: read_iniprint
-
-    namelist /conf_gcm_nml/dayref, anneeref, raz_date, nday, day_step, &
-         iperiod, iapp_tracvl, iconser, iecri, periodav, &
-         iphysiq, clon, clat, grossismx, grossismy, dzoomx, dzoomy, taux, &
-         tauy, offline
+    use unit_nml_m, only: unit_nml
+
+    namelist /conf_gcm_nml/ raz_date, nday, day_step, iperiod, iapp_tracvl, &
+         iconser, iecri, periodav, iphysiq, offline
+
+    namelist /iniprint_nml/ prt_level
+
+    namelist /logic_nml/ purmats, iflag_phys
 
     !------------------------------------
 
     print *, "Call sequence information: conf_gcm"
 
-    call read_iniprint
-    call read_logic
+    print *, "Enter namelist 'iniprint_nml'."
+    read(unit=*, nml=iniprint_nml)
+    write(unit_nml, nml=iniprint_nml)
+
+    print *, "Enter namelist 'logic_nml'."
+    read(unit=*, nml=logic_nml)
+    write(unit_nml, nml=logic_nml)
+
     call read_comdissnew
 
     print *, "Enter namelist 'conf_gcm_nml'."
     read(unit=*, nml=conf_gcm_nml)
-    write(unit=*, nml=conf_gcm_nml)
+    write(unit_nml, nml=conf_gcm_nml)
 
     IF (MOD(day_step, iperiod) /= 0) call abort_gcm(modname = "conf_gcm", &
          message = &
          'Il faut choisir un nombre de pas par jour multiple de "iperiod".', &
          ierr = 1)
 
-    IF (MOD(day_step,iphysiq)/= 0) call abort_gcm(modname = "conf_gcm", &
+    IF (MOD(day_step, iphysiq)/= 0) call abort_gcm(modname = "conf_gcm", &
          message = &
          'Il faut choisir un nombre de pas par jour multiple de "iphysiq".', &
          ierr = 1)
 
-    IF (grossismx < 1.) THEN
-       PRINT *, 'Error: grossismx < 1'
-       STOP 1
-    ELSE
-       alphax = 1. - 1. / grossismx
-    ENDIF
-    IF (grossismy < 1.) THEN
-       PRINT *, 'Error: grossismy < 1'
-       STOP 1
-    ELSE
-       alphay = 1. - 1. / grossismy
-    ENDIF
-    PRINT *, 'alphax = ', alphax
-    PRINT *, 'alphay = ', alphay
-
   END SUBROUTINE conf_gcm
 
 end module conf_gcm_m