Changeset 1140

Timestamp:

Mar 30, 2009, 4:46:54 PM (15 years ago)

Author:

Ehouarn Millour

Message:

Premiere vaque de modifications pour l'unification des dynamiques (planetes-Terre) et un peu de netoyage ...

• modified 'makegcm' and 'makegcm_fcm' to remove 'CPP_PHYS' key and add 'CPP_EARTH' preprocessing key instead
• updated 'diagedyn.F' (in dyn3d and dyn3dpar) to use 'CPP_EARTH' key
• added 'ioipsl_getincom.F90' and 'ioipsl_stringop.F90' to 'dyn3d' and 'dyn3dpar' for future possibility of running without IOIPSL library
• modified conf_gcm.F ( in d'yn3d' and 'dyn3dpar') to read in flag 'planet_type' (default=='earth') (flag added in 'control.h')
• modified 'gcm.F' (in 'dyn3d' and 'dyn3dpar') so that flags so that 'read_start' and 'iflag_phys' (known from conf_gcm.F) are used
• added flag 'output_grads_dyn' (read by conf_gcm.F, stored in 'control.h') to write grads outputs from 'leapfrog.F' and 'leapfrog_p.F'
• removed 'comdiss.h' from 'dyn3d' and 'dyn3dpar' (it is not used)
• removed variable 'lstardis' from 'comdissip.h' (it is also in

'comdissnew.h'), in dyn3d as well as in dyn3dpar

• adapted 'dyn3d/iniacademic.F' to not use 'inicons0.F' but 'iniconst.F'
• updated 'dyn3d/etat0_netcdf.F' to not use 'inicons0' but 'iniconst' (added prerequisite pa=50000 instruction) and added #ifdef CPP_EARTH keys
• removed 'inicons0.F' and 'disvert0.F' (not used any more)

Location:

LMDZ4/branches/LMDZ4-dev

Files:

• libf/dyn3d/comdiss.h (deleted)
• libf/dyn3d/comdissip.h (modified) (1 diff)
• libf/dyn3d/conf_gcm.F (modified) (24 diffs)
• libf/dyn3d/control.h (modified) (2 diffs)
• libf/dyn3d/diagedyn.F (modified) (3 diffs)
• libf/dyn3d/disvert0.F (deleted)
• libf/dyn3d/etat0_netcdf.F (modified) (7 diffs)
• libf/dyn3d/gcm.F (modified) (16 diffs)
• libf/dyn3d/iniacademic.F (modified) (5 diffs)
• libf/dyn3d/inicons0.F (deleted)
• libf/dyn3d/ioipsl_getincom.F90 (added)
• libf/dyn3d/ioipsl_stringop.F90 (added)
• libf/dyn3d/leapfrog.F (modified) (16 diffs)
• libf/dyn3dpar/comdiss.h (deleted)
• libf/dyn3dpar/comdissip.h (modified) (1 diff)
• libf/dyn3dpar/conf_gcm.F (modified) (30 diffs)
• libf/dyn3dpar/control.h (modified) (2 diffs)
• libf/dyn3dpar/diagedyn.F (modified) (3 diffs)
• libf/dyn3dpar/disvert0.F (deleted)
• libf/dyn3dpar/etat0_netcdf.F (modified) (7 diffs)
• libf/dyn3dpar/gcm.F (modified) (12 diffs)
• libf/dyn3dpar/iniacademic.F (modified) (5 diffs)
• libf/dyn3dpar/inicons0.F (deleted)
• libf/dyn3dpar/ioipsl_getincom.F90 (added)
• libf/dyn3dpar/ioipsl_stringop.F90 (added)
• libf/dyn3dpar/leapfrog_p.F (modified) (20 diffs)
• makegcm (modified) (7 diffs)
• makelmdz_fcm (modified) (5 diffs)

LMDZ4/branches/LMDZ4-dev/libf/dyn3d/comdissip.h

@@ -2 +2 @@
 ! $Header$
 !
-c-----------------------------------------------------------------------
-c INCLUDE dissip.h
+!-----------------------------------------------------------------------
+! INCLUDE comdissip.h
 
-      COMMON/comdissip/
-     $    lstardis,niterdis,coefdis,tetavel,tetatemp,gamdissip
+      COMMON/comdissip/                                                 &
+     &    niterdis,coefdis,tetavel,tetatemp,gamdissip
 
 
-      LOGICAL lstardis
       INTEGER niterdis
 
       REAL tetavel,tetatemp,coefdis,gamdissip
 
-c-----------------------------------------------------------------------
+!-----------------------------------------------------------------------

LMDZ4/branches/LMDZ4-dev/libf/dyn3d/conf_gcm.F

@@ -6 +6 @@
       SUBROUTINE conf_gcm( tapedef, etatinit, clesphy0 )
 c
+#ifdef CPP_IOIPSL
       use IOIPSL
+#else
+! if not using IOIPSL, we still need to use (a local version of) getin
+      use ioipsl_getincom
+#endif
       IMPLICIT NONE
 c-----------------------------------------------------------------------
@@ -99 +104 @@
 c  Parametres de controle du run:
 c-----------------------------------------------------------------------
+!Config  Key  = planet_type
+!Config  Desc = planet type ("earth", "mars", "venus", ...)
+!Config  Def  = earth
+!Config  Help = this flag sets the type of atymosphere that is considered
+      planet_type="earth"
+      CALL getin('planet_type',planet_type)
 
 !Config  Key  = dayref
@@ -179 +190 @@
        CALL getin('periodav',periodav)
 
+!Config  Key  = output_grads_dyn
+!Config  Desc = output dynamics diagnostics in 'dyn.dat' file
+!Config  Def  = n
+!Config  Help = output dynamics diagnostics in Grads-readable 'dyn.dat' file
+       output_grads_dyn=.false.
+       CALL getin('output_grads_dyn',output_grads_dyn)
+
 !Config  Key  = idissip
 !Config  Desc = periode de la dissipation 
@@ -274 +292 @@
 c    ...............................................................
 
+!Config  Key  =  read_start
+!Config  Desc = Initialize model using a 'start.nc' file
+!Config  Def  = y
+!Config  Help = y: intialize dynamical fields using a 'start.nc' file
+!               n: fields are initialized by 'iniacademic' routine
+       read_start= .true.
+       CALL getin('read_start',read_start)
+
 !Config  Key  = iflag_phys
 !Config  Desc = Avec ls physique 
@@ -330 +356 @@
 c
       IF( ABS(clat - clatt).GE. 0.001 )  THEN
-        PRINT *,' La valeur de clat passee par run.def est differente de 
-     *  celle lue sur le fichier  start '
+        write(lunout,*)'conf_gcm: La valeur de clat passee par run.def',
+     &    ' est differente de celle lue sur le fichier  start '
         STOP
       ENDIF
@@ -345 +371 @@
 
       IF( ABS(grossismx - grossismxx).GE. 0.001 )  THEN
-        PRINT *,' La valeur de grossismx passee par run.def est differente  
-     *  de celle lue sur le fichier  start '
+        write(lunout,*)'conf_gcm: La valeur de grossismx passee par ',
+     &  'run.def est differente de celle lue sur le fichier  start '
         STOP
       ENDIF
@@ -359 +385 @@
 
       IF( ABS(grossismy - grossismyy).GE. 0.001 )  THEN
-        PRINT *,' La valeur de grossismy passee par run.def est differen
-     * te de celle lue sur le fichier  start '
+        write(lunout,*)'conf_gcm: La valeur de grossismy passee par ',
+     & 'run.def est differente de celle lue sur le fichier  start '
         STOP
       ENDIF
       
       IF( grossismx.LT.1. )  THEN
-        PRINT *,' ***  ATTENTION !! grossismx < 1 .   *** '
+        write(lunout,*)
+     &       'conf_gcm: ***  ATTENTION !! grossismx < 1 .   *** '
          STOP
       ELSE
@@ -373 +400 @@
 
       IF( grossismy.LT.1. )  THEN
-        PRINT *,' ***  ATTENTION !! grossismy < 1 .   *** '
+        write(lunout,*)
+     &       'conf_gcm: ***  ATTENTION !! grossismy < 1 .   *** '
          STOP
       ELSE
@@ -379 +407 @@
       ENDIF
 
-      PRINT *,' alphax alphay defrun ',alphax,alphay
+      write(lunout,*)'conf_gcm: alphax alphay',alphax,alphay
 c
 c    alphax et alphay sont les anciennes formulat. des grossissements
@@ -394 +422 @@
 
       IF( .NOT.fxyhypb )  THEN
-           IF( fxyhypbb )     THEN
-              PRINT *,' ********  PBS DANS  DEFRUN  ******** '
-              PRINT *,' *** fxyhypb lu sur le fichier start est F ',
-     *       'alors  qu il est  T  sur  run.def  ***'
+         IF( fxyhypbb )     THEN
+            write(lunout,*)' ********  PBS DANS  CONF_GCM  ******** '
+            write(lunout,*)' *** fxyhypb lu sur le fichier start est ',
+     *       'F alors  qu il est  T  sur  run.def  ***'
               STOP
-           ENDIF
+         ENDIF
       ELSE
-           IF( .NOT.fxyhypbb )   THEN
-              PRINT *,' ********  PBS DANS  DEFRUN  ******** '
-              PRINT *,' ***  fxyhypb lu sur le fichier start est T ',
-     *        'alors  qu il est  F  sur  run.def  ****  '
+         IF( .NOT.fxyhypbb )   THEN
+            write(lunout,*)' ********  PBS DANS  CONF_GCM  ******** '
+            write(lunout,*)' ***  fxyhypb lu sur le fichier start est ',
+     *        'T alors  qu il est  F  sur  run.def  ****  '
               STOP
-           ENDIF
+         ENDIF
       ENDIF
 c
@@ -419 +447 @@
       IF( fxyhypb )  THEN
        IF( ABS(dzoomx - dzoomxx).GE. 0.001 )  THEN
-        PRINT *,' La valeur de dzoomx passee par run.def est differente 
-     *  de celle lue sur le fichier  start '
+        write(lunout,*)'conf_gcm: La valeur de dzoomx passee par ',
+     *  'run.def est differente de celle lue sur le fichier  start '
         STOP
        ENDIF
@@ -435 +463 @@
       IF( fxyhypb )  THEN
        IF( ABS(dzoomy - dzoomyy).GE. 0.001 )  THEN
-        PRINT *,' La valeur de dzoomy passee par run.def est differente 
-     * de celle lue sur le fichier  start '
+        write(lunout,*)'conf_gcm: La valeur de dzoomy passee par ',
+     * 'run.def est differente de celle lue sur le fichier  start '
         STOP
        ENDIF
@@ -450 +478 @@
       IF( fxyhypb )  THEN
        IF( ABS(taux - tauxx).GE. 0.001 )  THEN
-        PRINT *,' La valeur de taux passee par run.def est differente
-     * de celle lue sur le fichier  start '
+        write(lunout,*)'conf_gcm: La valeur de taux passee par ',
+     * 'run.def est differente de celle lue sur le fichier  start '
         STOP
        ENDIF
@@ -465 +493 @@
       IF( fxyhypb )  THEN
        IF( ABS(tauy - tauyy).GE. 0.001 )  THEN
-        PRINT *,' La valeur de tauy passee par run.def est differente
-     * de celle lue sur le fichier  start '
+        write(lunout,*)'conf_gcm: La valeur de tauy passee par ',
+     * 'run.def est differente de celle lue sur le fichier  start '
         STOP
        ENDIF
@@ -484 +512 @@
 
         IF( .NOT.ysinus )  THEN
-           IF( ysinuss )     THEN
-              PRINT *,' ********  PBS DANS  DEFRUN  ******** '
-              PRINT *,' *** ysinus lu sur le fichier start est F ',
-     *       'alors  qu il est  T  sur  run.def  ***'
+          IF( ysinuss )     THEN
+            write(lunout,*)' ********  PBS DANS  CONF_GCM  ******** '
+            write(lunout,*)' *** ysinus lu sur le fichier start est F',
+     *       ' alors  qu il est  T  sur  run.def  ***'
+            STOP
+          ENDIF
+        ELSE
+          IF( .NOT.ysinuss )   THEN
+            write(lunout,*)' ********  PBS DANS  CONF_GCM  ******** '
+            write(lunout,*)' *** ysinus lu sur le fichier start est T',
+     *        ' alors  qu il est  F  sur  run.def  ****  '
               STOP
-           ENDIF
-        ELSE
-           IF( .NOT.ysinuss )   THEN
-              PRINT *,' ********  PBS DANS  DEFRUN  ******** '
-              PRINT *,' ***  ysinus lu sur le fichier start est T ',
-     *        'alors  qu il est  F  sur  run.def  ****  '
-              STOP
-           ENDIF
+          ENDIF
         ENDIF
-      ENDIF
+      ENDIF ! of IF( .NOT.fxyhypb  )
 c
 !Config  Key  = offline
@@ -521 +549 @@
       write(lunout,*)' #########################################'
       write(lunout,*)' Configuration des parametres du gcm: '
+      write(lunout,*)' planet_type = ', planet_type
       write(lunout,*)' dayref = ', dayref
       write(lunout,*)' anneeref = ', anneeref
@@ -529 +558 @@
       write(lunout,*)' iecri = ', iecri
       write(lunout,*)' periodav = ', periodav 
+      write(lunout,*)' output_grads_dyn = ', output_grads_dyn
       write(lunout,*)' idissip = ', idissip
       write(lunout,*)' lstardis = ', lstardis
@@ -539 +569 @@
       write(lunout,*)' coefdis = ', coefdis
       write(lunout,*)' purmats = ', purmats
+      write(lunout,*)' read_start = ', read_start
       write(lunout,*)' iflag_phys = ', iflag_phys
       write(lunout,*)' iphysiq = ', iphysiq
@@ -590 +621 @@
 
       IF( grossismx.LT.1. )  THEN
-        PRINT *,' ***  ATTENTION !! grossismx < 1 .   *** '
+        write(lunout,*)
+     &   'conf_gcm: ***  ATTENTION !! grossismx < 1 .   *** '
          STOP
       ELSE
@@ -598 +630 @@
 
       IF( grossismy.LT.1. )  THEN
-        PRINT *,' ***  ATTENTION !! grossismy < 1 .   *** '
+        write(lunout,*)
+     &  'conf_gcm: ***  ATTENTION !! grossismy < 1 .   *** '
          STOP
       ELSE
@@ -604 +637 @@
       ENDIF
 
-      PRINT *,' alphax alphay defrun ',alphax,alphay
+      write(lunout,*)'conf_gcm: alphax alphay ',alphax,alphay
 c
 c    alphax et alphay sont les anciennes formulat. des grossissements
@@ -693 +726 @@
       write(lunout,*)' #########################################'
       write(lunout,*)' Configuration des parametres du gcm: '
+      write(lunout,*)' planet_type = ', planet_type
       write(lunout,*)' dayref = ', dayref
       write(lunout,*)' anneeref = ', anneeref
@@ -701 +735 @@
       write(lunout,*)' iecri = ', iecri
       write(lunout,*)' periodav = ', periodav 
+      write(lunout,*)' output_grads_dyn = ', output_grads_dyn
       write(lunout,*)' idissip = ', idissip
       write(lunout,*)' lstardis = ', lstardis
@@ -711 +746 @@
       write(lunout,*)' coefdis = ', coefdis
       write(lunout,*)' purmats = ', purmats
+      write(lunout,*)' read_start = ', read_start
       write(lunout,*)' iflag_phys = ', iflag_phys
       write(lunout,*)' iphysiq = ', iphysiq
-      write(lunout,*)' clonn = ', clonn 
-      write(lunout,*)' clatt = ', clatt
+      write(lunout,*)' clon = ', clon
+      write(lunout,*)' clat = ', clat
       write(lunout,*)' grossismx = ', grossismx
       write(lunout,*)' grossismy = ', grossismy
-      write(lunout,*)' fxyhypbb = ', fxyhypbb
+      write(lunout,*)' fxyhypb = ', fxyhypb
       write(lunout,*)' dzoomx = ', dzoomx
       write(lunout,*)' dzoomy = ', dzoomy

LMDZ4/branches/LMDZ4-dev/libf/dyn3d/control.h

@@ -14 +14 @@
      &              iperiod,iapp_tracvl,iconser,iecri,idissip,iphysiq , &
      &              periodav,iecrimoy,dayref,anneeref,                  &
-     &              raz_date,offline,ip_ebil_dyn,config_inca
+     &              raz_date,offline,ip_ebil_dyn,config_inca,           &
+     &              planet_type,output_grads_dyn
 
       INTEGER   nday,day_step,iperiod,iapp_tracvl,iconser,iecri,        &
@@ -22 +23 @@
       logical offline
       CHARACTER (len=4) :: config_inca
+      CHARACTER(len=10) :: planet_type ! planet type ('earth','mars',...)
+      LOGICAL :: output_grads_dyn ! output dynamics diagnostics in
+                                  ! binary grads file 'dyn.dat' (y/n)
 !-----------------------------------------------------------------------

LMDZ4/branches/LMDZ4-dev/libf/dyn3d/diagedyn.F

@@ -58 +58 @@
 #include "paramet.h"
 #include "comgeom.h"
-
-#ifdef CPP_PHYS
+#include "iniprint.h"
+
+#ifdef CPP_EARTH
 #include "../phylmd/YOMCST.h"
 #include "../phylmd/YOETHF.h"
@@ -139 +140 @@
 
 
-#ifdef CPP_PHYS
+#ifdef CPP_EARTH
 c======================================================================
 C     Compute Kinetic enrgy
@@ -314 +315 @@
 C
 #else
-      print*,'Pour l instant diagedyn a besoin de la physique'
+      write(lunout,*),'diagedyn: Needs Earth physics to function'
 #endif
+! #endif of #ifdef CPP_EARTH 
       RETURN 
       END 

LMDZ4/branches/LMDZ4-dev/libf/dyn3d/etat0_netcdf.F

@@ -5 +5 @@
 c
       SUBROUTINE etat0_netcdf (interbar, masque)
-    
+#ifdef CPP_EARTH        
       USE startvar
       USE ioipsl
@@ -14 +14 @@
       USE phys_state_var_mod
       USE filtreg_mod
+#endif
+!#endif of #ifdef CPP_EARTH
       !
       IMPLICIT NONE
@@ -25 +27 @@
 !     .KLON=KFDIA-KIDIA+1,KLEV=llm
       !
+#ifdef CPP_EARTH    
 #include "comgeom2.h"
 #include "comvert.h"
@@ -31 +34 @@
 #include "dimsoil.h"
 #include "temps.h"
-      !
+#endif
+!#endif of #ifdef CPP_EARTH
+      ! arguments:
       LOGICAL interbar
+      REAL :: masque(iip1,jjp1)
+
+#ifdef CPP_EARTH
+      ! local variables:
       REAL :: latfi(klon), lonfi(klon)
-      REAL :: orog(iip1,jjp1), rugo(iip1,jjp1), masque(iip1,jjp1)
+      REAL :: orog(iip1,jjp1), rugo(iip1,jjp1)
       REAL :: psol(iip1, jjp1), phis(iip1, jjp1)
       REAL :: p3d(iip1, jjp1, llm+1)
@@ -143 +152 @@
       !
       preff     = 101325.
+      pa        =  50000.
       unskap = 1./kappa
       !
@@ -168 +178 @@
 
 
-      CALL inicons0()
+      CALL iniconst()
       CALL inigeom()
 
@@ -759 +769 @@
       DEALLOCATE(q3d)
 
+#endif 
+!#endif of #ifdef CPP_EARTH
       RETURN
       !

LMDZ4/branches/LMDZ4-dev/libf/dyn3d/gcm.F

@@ -8 +8 @@
 #ifdef CPP_IOIPSL
       USE IOIPSL
+#else
+! if not using IOIPSL, we still need to use (a local version of) getin
+      USE ioipsl_getincom
 #endif
 
@@ -17 +20 @@
 ! A nettoyer. On ne veut qu'une ou deux routines d'interface 
 ! dynamique -> physique pour l'initialisation
-#ifdef CPP_PHYS
+! Ehouarn: for now these only apply to Earth:
+#ifdef CPP_EARTH
       USE dimphy
       USE comgeomphy
@@ -158 +162 @@
 
 
-c--------------------------------------------------------------------------
-c   Iflag_phys controle l'appel a la physique :
-c   -------------------------------------------
-c      0 : pas de physique
-c      1 : Normale (appel a phylmd, phymars ...)
-c      2 : rappel Newtonien pour la temperature + friction au sol
-      iflag_phys=1
-
-c--------------------------------------------------------------------------
-c   Lecture de l'etat initial :
-c   ---------------------------
-c     T : on lit start.nc
-c     F : le modele s'autoinitialise avec un cas academique (iniacademic)
-      read_start=.true.
-#ifdef CPP_IOIPSL
-#else
-      read_start=.false.
-#endif
-#ifdef CPP_PHYS
-#else
-      read_start=.false.
-#endif
 
 c-----------------------------------------------------------------------
@@ -196 +178 @@
 c  ---------------------------------------
 c
-#ifdef CPP_IOIPSL
+! Ehouarn: dump possibility of using defrun
+!#ifdef CPP_IOIPSL
       CALL conf_gcm( 99, .TRUE. , clesphy0 )
-#else
-      CALL defrun( 99, .TRUE. , clesphy0 )
-#endif
+!#else
+!      CALL defrun( 99, .TRUE. , clesphy0 )
+!#endif
 
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
@@ -206 +189 @@
 ! A nettoyer. On ne veut qu'une ou deux routines d'interface 
 ! dynamique -> physique pour l'initialisation
-#ifdef CPP_PHYS
+! Ehouarn : temporarily (?) keep this only for Earth
+      if (planet_type.eq."earth") then
+#ifdef CPP_EARTH
       CALL Init_Phys_lmdz(iim,jjp1,llm,1,(jjm-1)*iim+2)
       call InitComgeomphy
 #endif
+      endif
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 
@@ -242 +228 @@
 c  lecture du fichier start.nc
       if (read_start) then
-#ifdef CPP_IOIPSL
+      ! we still need to run iniacademic to initialize some
+      ! constants & fields, if we run the 'newtonian' case:
+        if (iflag_phys.eq.2) then
+          CALL iniacademic(vcov,ucov,teta,q,masse,ps,phis,time_0)
+        endif
+!#ifdef CPP_IOIPSL
+        if (planet_type.eq."earth") then
+#ifdef CPP_EARTH
+! Load an Earth-format start file
          CALL dynetat0("start.nc",vcov,ucov,
      .              teta,q,masse,ps,phis, time_0)
+#endif
+        endif ! of if (planet_type.eq."earth")
 c       write(73,*) 'ucov',ucov
 c       write(74,*) 'vcov',vcov
@@ -251 +247 @@
 c       write(77,*) 'q',q
 
-#endif
-      endif
+      endif ! of if (read_start)
 
       IF (config_inca /= 'none') THEN
@@ -264 +259 @@
 c le cas echeant, creation d un etat initial
       IF (prt_level > 9) WRITE(lunout,*)
-     .                 'AVANT iniacademic AVANT AVANT AVANT AVANT'
+     .              'GCM: AVANT iniacademic AVANT AVANT AVANT AVANT'
       if (.not.read_start) then
          CALL iniacademic(vcov,ucov,teta,q,masse,ps,phis,time_0)
@@ -298 +293 @@
       if (annee_ref .ne. anneeref .or. day_ref .ne. dayref) then
         write(lunout,*)
-     .  ' Attention les dates initiales lues dans le fichier'
+     .  'GCM: Attention les dates initiales lues dans le fichier'
         write(lunout,*)
      .  ' restart ne correspondent pas a celles lues dans '
@@ -304 +299 @@
         if (raz_date .ne. 1) then
           write(lunout,*)
-     .    ' On garde les dates du fichier restart'
+     .    'GCM: On garde les dates du fichier restart'
         else
           annee_ref = anneeref
@@ -313 +308 @@
           time_0 = 0.
           write(lunout,*)
-     .   ' On reinitialise a la date lue dans gcm.def'
+     .   'GCM: On reinitialise a la date lue dans gcm.def'
         endif
       ELSE
@@ -350 +345 @@
 c   Initialisation de la physique :
 c   -------------------------------
-#ifdef CPP_PHYS
-      IF (call_iniphys.and.iflag_phys.eq.1) THEN
+
+      IF (call_iniphys.and.(iflag_phys.eq.1)) THEN
          latfi(1)=rlatu(1)
          lonfi(1)=0.
@@ -370 +365 @@
          CALL gr_dyn_fi(1,iip1,jjp1,ngridmx,aire,airefi)
          WRITE(lunout,*)
-     .           'WARNING!!! vitesse verticale nulle dans la physique'
+     .       'GCM: WARNING!!! vitesse verticale nulle dans la physique'
+! Earth:
+         if (planet_type.eq."earth") then
+#ifdef CPP_EARTH
          CALL iniphysiq(ngridmx,llm,daysec,day_ini,dtphys ,
      ,                latfi,lonfi,airefi,zcufi,zcvfi,rad,g,r,cpp     )
+#endif
+         endif ! of if (planet_type.eq."earth")
          call_iniphys=.false.
-      ENDIF
-#endif
+      ENDIF ! of IF (call_iniphys.and.(iflag_phys.eq.1))
+!#endif
 
 c  numero de stockage pour les fichiers de redemarrage:
@@ -386 +386 @@
       day_end = day_ini + nday
       WRITE(lunout,300)day_ini,day_end
+ 300  FORMAT('1'/,15x,'run du jour',i7,2x,'au jour',i7//)
+
+      if (planet_type.eq."earth") then
+#ifdef CPP_EARTH
+      CALL dynredem0("restart.nc", day_end, phis)
+#endif
+      endif
+
+      ecripar = .TRUE.
 
 #ifdef CPP_IOIPSL
-      CALL dynredem0("restart.nc", day_end, phis)
-
-      ecripar = .TRUE.
-
       if ( 1.eq.1) then
       time_step = zdtvr
@@ -409 +414 @@
 
 #endif
+! #endif of #ifdef CPP_IOIPSL
 
 c  Choix des frequences de stokage pour le offline
@@ -432 +438 @@
      .              time_0)
 
-
-
- 300  FORMAT('1'/,15x,'run du pas',i7,2x,'au pas',i7,2x,
-     . 'c''est a dire du jour',i7,3x,'au jour',i7//)
       END
 

LMDZ4/branches/LMDZ4-dev/libf/dyn3d/iniacademic.F

@@ -45 +45 @@
 #include "temps.h"
 #include "control.h"
+#include "iniprint.h"
 
 c   Arguments:
@@ -57 +58 @@
       REAL ps(ip1jmp1)                       ! pression  au sol
       REAL masse(ip1jmp1,llm)                ! masse d'air
+      REAL phis(ip1jmp1)                     ! geopotentiel au sol
+
+c   Local:
+c   ------
+
       REAL p (ip1jmp1,llmp1  )               ! pression aux interfac.des couches
       REAL pks(ip1jmp1)                      ! exner au  sol
       REAL pk(ip1jmp1,llm)                   ! exner au milieu des couches
       REAL pkf(ip1jmp1,llm)                  ! exner filt.au milieu des couches
-      REAL phis(ip1jmp1)                     ! geopotentiel au sol
       REAL phi(ip1jmp1,llm)                  ! geopotentiel
-
-
-
-
-
-c   Local:
-c   ------
-
       REAL ddsin,tetarappelj,tetarappell,zsig
       real tetajl(jjp1,llm)
@@ -81 +78 @@
 
 c-----------------------------------------------------------------------
+! 1. Initializations for Earth-like case
+! --------------------------------------
+      if (planet_type=="earth") then
+c
+        time_0=0.
 
-c
-      time_0=0.
+        im         = iim
+        jm         = jjm
+        day_ini    = 0
+        omeg       = 4.*asin(1.)/86400.
+        rad    = 6371229.
+        g      = 9.8
+        daysec = 86400.
+        dtvr    = daysec/FLOAT(day_step)
+        zdtvr=dtvr
+        kappa  = 0.2857143
+        cpp    = 1004.70885
+        preff     = 101325.
+        pa        =  50000.
+        etot0      = 0.
+        ptot0      = 0.
+        ztot0      = 0.
+        stot0      = 0.
+        ang0       = 0.
 
-      im         = iim
-      jm         = jjm
-      day_ini    = 0
-      omeg       = 4.*asin(1.)/86400.
-      rad    = 6371229.
-      g      = 9.8
-      daysec = 86400.
-      dtvr    = daysec/FLOAT(day_step)
-      zdtvr=dtvr
-      kappa  = 0.2857143
-      cpp    = 1004.70885
-      preff     = 101325.
-      pa        =  50 000.
-      etot0      = 0.
-      ptot0      = 0.
-      ztot0      = 0.
-      stot0      = 0.
-      ang0       = 0.
-      pa         = 0.
+        CALL iniconst
+        CALL inigeom
+        CALL inifilr
 
-      CALL inicons0
-      CALL inigeom
-      CALL inifilr
-
-      ps=0.
-      phis=0.
+        ps=0.
+        phis=0.
 c---------------------------------------------------------------------
 
-      taurappel=10.*daysec
+        taurappel=10.*daysec
 
 c---------------------------------------------------------------------
@@ -119 +117 @@
 c   --------------------------------------
 
-      DO l=1,llm
-       zsig=ap(l)/preff+bp(l)
-       if (zsig.gt.0.3) then
-         lsup=l
-         tetarappell=1./8.*(-log(zsig)-.5)
-         DO j=1,jjp1
-            ddsin=sin(rlatu(j))-sin(pi/20.)
-            tetajl(j,l)=300.*(1+1./18.*(1.-3.*ddsin*ddsin)+tetarappell)
-         ENDDO
-        else
+        DO l=1,llm
+         zsig=ap(l)/preff+bp(l)
+         if (zsig.gt.0.3) then
+           lsup=l
+           tetarappell=1./8.*(-log(zsig)-.5)
+           DO j=1,jjp1
+             ddsin=sin(rlatu(j))-sin(pi/20.)
+             tetajl(j,l)=300.*(1+1./18.*(1.-3.*ddsin*ddsin)+tetarappell)
+           ENDDO
+          else
 c   Choix isotherme au-dessus de 300 mbar
-         do j=1,jjp1
-            tetajl(j,l)=tetajl(j,lsup)*(0.3/zsig)**kappa
-         enddo
-        endif
-      ENDDO
+           do j=1,jjp1
+             tetajl(j,l)=tetajl(j,lsup)*(0.3/zsig)**kappa
+           enddo
+          endif ! of if (zsig.gt.0.3)
+        ENDDO ! of DO l=1,llm
 
-      do l=1,llm
-         do j=1,jjp1
-            do i=1,iip1
-               ij=(j-1)*iip1+i
-               tetarappel(ij,l)=tetajl(j,l)
-            enddo
-         enddo
-      enddo
+        do l=1,llm
+           do j=1,jjp1
+              do i=1,iip1
+                 ij=(j-1)*iip1+i
+                 tetarappel(ij,l)=tetajl(j,l)
+              enddo
+           enddo
+        enddo
 
-c     call dump2d(jjp1,llm,tetajl,'TEQ   ')
+c       call dump2d(jjp1,llm,tetajl,'TEQ   ')
 
-      ps=1.e5
-      phis=0.
-      CALL pression ( ip1jmp1, ap, bp, ps, p       )
-      CALL exner_hyb( ip1jmp1, ps, p,alpha,beta, pks, pk, pkf )
-      CALL massdair(p,masse)
+        ps=1.e5
+        phis=0.
+        CALL pression ( ip1jmp1, ap, bp, ps, p       )
+        CALL exner_hyb( ip1jmp1, ps, p,alpha,beta, pks, pk, pkf )
+        CALL massdair(p,masse)
 
 c  intialisation du vent et de la temperature
-      teta(:,:)=tetarappel(:,:)
-      CALL geopot(ip1jmp1,teta,pk,pks,phis,phi)
-      call ugeostr(phi,ucov)
-      vcov=0.
-      q(:,:,1   )=1.e-10
-      q(:,:,2   )=1.e-15
-      q(:,:,3:nqtot)=0.
+        teta(:,:)=tetarappel(:,:)
+        CALL geopot(ip1jmp1,teta,pk,pks,phis,phi)
+        call ugeostr(phi,ucov)
+        vcov=0.
+        q(:,:,1   )=1.e-10
+        q(:,:,2   )=1.e-15
+        q(:,:,3:nqtot)=0.
 
 
-c   perturbation al\351atoire sur la temp\351rature
-      idum  = -1
-      zz = ran1(idum)
-      idum  = 0
-      do l=1,llm
-         do ij=iip2,ip1jm
-            teta(ij,l)=teta(ij,l)*(1.+0.005*ran1(idum))
-         enddo
-      enddo
+c   perturbation aleatoire sur la temperature
+        idum  = -1
+        zz = ran1(idum)
+        idum  = 0
+        do l=1,llm
+           do ij=iip2,ip1jm
+              teta(ij,l)=teta(ij,l)*(1.+0.005*ran1(idum))
+           enddo
+        enddo
 
-      do l=1,llm
-         do ij=1,ip1jmp1,iip1
-            teta(ij+iim,l)=teta(ij,l)
-         enddo
-      enddo
+        do l=1,llm
+           do ij=1,ip1jmp1,iip1
+              teta(ij+iim,l)=teta(ij,l)
+           enddo
+        enddo
 
 
@@ -187 +185 @@
 
 c   initialisation d'un traceur sur une colonne
-      j=jjp1*3/4
-      i=iip1/2
-      ij=(j-1)*iip1+i
-      q(ij,:,3)=1.
-
+        j=jjp1*3/4
+        i=iip1/2
+        ij=(j-1)*iip1+i
+        q(ij,:,3)=1.
+      
+      else
+        write(lunout,*)"iniacademic: planet types other than earth",
+     &                 " not implemented (yet)."
+        stop
+      endif ! of if (planet_type=="earth")
       return
       END

LMDZ4/branches/LMDZ4-dev/libf/dyn3d/leapfrog.F

@@ -7 +7 @@
 
 cIM : pour sortir les param. du modele dans un fis. netcdf 110106
-      USE IOIPSL 
+#ifdef CPP_IOIPSL
+      use IOIPSL
+#endif
       USE infotrac
 
@@ -161 +163 @@
 
       character*80 dynhist_file, dynhistave_file
-      character*20 modname
+      character(len=20) :: modname
       character*80 abort_message
 
@@ -217 +219 @@
         call guide(itau,ucov,vcov,teta,q,masse,ps)
       else
-        IF(prt_level>9)WRITE(*,*)'attention on ne guide pas les ',
-     .    '6 dernieres heures'
+        IF(prt_level>9)WRITE(lunout,*)'leapfrog: attention on ne ',
+     .    'guide pas les 6 dernieres heures'
       endif
 #endif
@@ -227 +229 @@
 c     ENDIF 
 c
+
+! Save fields obtained at previous time step as '...m1'
       CALL SCOPY( ijmllm ,vcov , 1, vcovm1 , 1 )
       CALL SCOPY( ijp1llm,ucov , 1, ucovm1 , 1 )
@@ -242 +246 @@
       CALL filtreg ( finvmaold ,jjp1, llm, -2,2, .TRUE., 1 )
 
-      call minmax(ijp1llm,q(:,:,3),zqmin,zqmax)
+! Ehouarn: what is this for? zqmin & zqmax are not used anyway ...
+!      call minmax(ijp1llm,q(:,:,3),zqmin,zqmax)
 
    2  CONTINUE
@@ -302 +307 @@
 
 
-         ENDIF
-c
-      ENDIF
+         ENDIF ! of IF (offline)
+c
+      ENDIF ! of IF( forward. OR . leapf )
 
 
@@ -350 +355 @@
 c   -----------------------------------------------------
 
-#ifdef CPP_PHYS
 c+jld
 
 c  Diagnostique de conservation de l'énergie : initialisation
-      IF (ip_ebil_dyn.ge.1 ) THEN 
+         IF (ip_ebil_dyn.ge.1 ) THEN 
           ztit='bil dyn'
-          CALL diagedyn(ztit,2,1,1,dtphys
-     e  , ucov    , vcov , ps, p ,pk , teta , q(:,:,1), q(:,:,2))
-      ENDIF 
+! Ehouarn: be careful, diagedyn is Earth-specific (includes ../phylmd/..)!
+           IF (planet_type.eq."earth") THEN
+            CALL diagedyn(ztit,2,1,1,dtphys
+     &    , ucov    , vcov , ps, p ,pk , teta , q(:,:,1), q(:,:,2))
+           ENDIF
+         ENDIF ! of IF (ip_ebil_dyn.ge.1 )
 c-jld
+#ifdef CPP_IOIPSL
 cIM : pour sortir les param. du modele dans un fis. netcdf 110106
-      IF (first) THEN
-       first=.false.
+         IF (first) THEN
+          first=.false.
 #include "ini_paramLMDZ_dyn.h"
-      ENDIF
+         ENDIF
 c
 #include "write_paramLMDZ_dyn.h"
 c
-
-        CALL calfis( lafin ,rdayvrai,time  ,
+#endif
+! #endif of #ifdef CPP_IOIPSL
+         CALL calfis( lafin ,rdayvrai,time  ,
      $               ucov,vcov,teta,q,masse,ps,p,pk,phis,phi ,
      $               du,dv,dteta,dq,
@@ -375 +384 @@
      $               clesphy0, dufi,dvfi,dtetafi,dqfi,dpfi  )
 
-       IF (ok_strato) THEN
-         CALL top_bound( vcov,ucov,teta, dufi,dvfi,dtetafi)
-       ENDIF
+         IF (ok_strato) THEN
+           CALL top_bound( vcov,ucov,teta, dufi,dvfi,dtetafi)
+         ENDIF
        
 c      ajout des tendances physiques:
@@ -386 +395 @@
 c
 c  Diagnostique de conservation de l'énergie : difference
-      IF (ip_ebil_dyn.ge.1 ) THEN 
+         IF (ip_ebil_dyn.ge.1 ) THEN 
           ztit='bil phys'
-          CALL diagedyn(ztit,2,1,1,dtphys
-     e  , ucov    , vcov , ps, p ,pk , teta , q(:,:,1), q(:,:,2))
-      ENDIF 
-#endif
+          IF (planet_type.eq."earth") THEN
+           CALL diagedyn(ztit,2,1,1,dtphys
+     &     , ucov    , vcov , ps, p ,pk , teta , q(:,:,1), q(:,:,2))
+          ENDIF
+         ENDIF ! of IF (ip_ebil_dyn.ge.1 )
+
        ENDIF ! of IF( apphys )
 
-      IF(iflag_phys.EQ.2) THEN ! "Newtonian physics" case
+      IF(iflag_phys.EQ.2) THEN ! "Newtonian" case
 c   Calcul academique de la physique = Rappel Newtonien + friction 
 c   --------------------------------------------------------------
@@ -472 +483 @@
 
 
-      END IF
+      END IF ! of IF(apdiss)
 
 c ajout debug
@@ -545 +556 @@
 c           IF( MOD(itau,iecri*day_step).EQ.0) THEN
 
-               nbetat = nbetatdem
-       CALL geopot  ( ip1jmp1, teta  , pk , pks,  phis  , phi        )
-        unat=0.
-        do l=1,llm
-           unat(iip2:ip1jm,l)=ucov(iip2:ip1jm,l)/cu(iip2:ip1jm)
-           vnat(:,l)=vcov(:,l)/cv(:)
-        enddo
-#ifdef CPP_IOIPSL
-c        CALL writehist(histid,histvid,itau,vcov, 
-c     s                       ucov,teta,phi,q,masse,ps,phis)
-#else
+              nbetat = nbetatdem
+              CALL geopot(ip1jmp1,teta,pk,pks,phis,phi)
+              unat=0.
+              do l=1,llm
+                unat(iip2:ip1jm,l)=ucov(iip2:ip1jm,l)/cu(iip2:ip1jm)
+                vnat(:,l)=vcov(:,l)/cv(:)
+              enddo
+#ifdef CPP_IOIPSL
+c             CALL writehist(histid,histvid,itau,vcov, 
+c     &                      ucov,teta,phi,q,masse,ps,phis)
+#endif
+! For some Grads outputs of fields
+             if (output_grads_dyn) then
 #include "write_grads_dyn.h"
-#endif
-
-
-            ENDIF
+             endif
+
+            ENDIF ! of IF(MOD(itau,iecri).EQ.0)
 
             IF(itau.EQ.itaufin) THEN
 
 
-#ifdef CPP_IOIPSL
-       CALL dynredem1("restart.nc",0.0,
-     ,                     vcov,ucov,teta,q,masse,ps)
-#endif
+              if (planet_type.eq."earth") then
+#ifdef CPP_EARTH
+! Write an Earth-format restart file
+                CALL dynredem1("restart.nc",0.0,
+     &                         vcov,ucov,teta,q,masse,ps)
+#endif
+              endif ! of if (planet_type.eq."earth")
 
               CLOSE(99)
-            ENDIF
+            ENDIF ! of IF (itau.EQ.itaufin)
 
 c-----------------------------------------------------------------------
@@ -593 +608 @@
                         leapf =  .TRUE.
                         dt  =  2.*dtvr
-                        GO TO 2
-                   END IF
+                        GO TO 2 
+                   END IF ! of IF (forward)
             ELSE
 
@@ -602 +617 @@
                  dt  = 2.*dtvr
                  GO TO 2
-            END IF
-
-      ELSE
+            END IF ! of IF (MOD(itau,iperiod).EQ.0)
+                   !    ELSEIF (MOD(itau-1,iperiod).EQ.0)
+
+      ELSE ! of IF (.not.purmats)
 
 c       ........................................................
@@ -627 +643 @@
                GO TO 2
 
-            ELSE
-
-            IF(MOD(itau,iperiod).EQ.0 .OR. itau.EQ.itaufin) THEN
+            ELSE ! of IF(forward)
+
+              IF(MOD(itau,iperiod).EQ.0 .OR. itau.EQ.itaufin) THEN
                IF(itau.EQ.itaufin) THEN
                   iav=1
@@ -636 +652 @@
                ENDIF
 #ifdef CPP_IOIPSL
-              CALL writedynav(histaveid, itau,vcov ,
+               CALL writedynav(histaveid, itau,vcov ,
      ,                          ucov,teta,pk,phi,q,masse,ps,phis)
                call bilan_dyn (2,dtvr*iperiod,dtvr*day_step*periodav,
@@ -642 +658 @@
 #endif
 
-            ENDIF
-
-               IF(MOD(itau,iecri         ).EQ.0) THEN
+              ENDIF ! of IF(MOD(itau,iperiod).EQ.0 .OR. itau.EQ.itaufin)
+
+              IF(MOD(itau,iecri         ).EQ.0) THEN
 c              IF(MOD(itau,iecri*day_step).EQ.0) THEN
-                  nbetat = nbetatdem
-       CALL geopot  ( ip1jmp1, teta  , pk , pks,  phis  , phi       )
-        unat=0.
-        do l=1,llm
-           unat(iip2:ip1jm,l)=ucov(iip2:ip1jm,l)/cu(iip2:ip1jm)
-           vnat(:,l)=vcov(:,l)/cv(:)
-        enddo
-#ifdef CPP_IOIPSL
-c       CALL writehist( histid, histvid, itau,vcov , 
-c    ,                           ucov,teta,phi,q,masse,ps,phis)
-#else
+                nbetat = nbetatdem
+                CALL geopot(ip1jmp1,teta,pk,pks,phis,phi)
+                unat=0.
+                do l=1,llm
+                  unat(iip2:ip1jm,l)=ucov(iip2:ip1jm,l)/cu(iip2:ip1jm)
+                  vnat(:,l)=vcov(:,l)/cv(:)
+                enddo
+#ifdef CPP_IOIPSL
+c               CALL writehist( histid, histvid, itau,vcov , 
+c    &                           ucov,teta,phi,q,masse,ps,phis)
+#endif
+! For some Grads outputs
+                if (output_grads_dyn) then
 #include "write_grads_dyn.h"
-#endif
-
-
-               ENDIF
-
-#ifdef CPP_IOIPSL
-                 IF(itau.EQ.itaufin)
-     . CALL dynredem1("restart.nc",0.0,
-     .                     vcov,ucov,teta,q,masse,ps)
-#endif
-
-                 forward = .TRUE.
-                 GO TO  1
-
-            ENDIF
-
-      END IF
+                endif
+
+              ENDIF ! of IF(MOD(itau,iecri         ).EQ.0) 
+
+              IF(itau.EQ.itaufin) THEN
+                if (planet_type.eq."earth") then
+#ifdef CPP_EARTH
+                  CALL dynredem1("restart.nc",0.0,
+     &                           vcov,ucov,teta,q,masse,ps)
+#endif
+                endif ! of if (planet_type.eq."earth")
+              ENDIF ! of IF(itau.EQ.itaufin)
+
+              forward = .TRUE.
+              GO TO  1
+
+            ENDIF ! of IF (forward)
+
+      END IF ! of IF(.not.purmats)
 
       STOP

LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/comdissip.h

@@ -2 +2 @@
 ! $Header$
 !
-c-----------------------------------------------------------------------
-c INCLUDE dissip.h
+!-----------------------------------------------------------------------
+! INCLUDE comdissip.h
 
-      COMMON/comdissip/
-     $    lstardis,niterdis,coefdis,tetavel,tetatemp,gamdissip
+      COMMON/comdissip/                                                 &
+     &    niterdis,coefdis,tetavel,tetatemp,gamdissip
 
 
-      LOGICAL lstardis
       INTEGER niterdis
 
       REAL tetavel,tetatemp,coefdis,gamdissip
 
-c-----------------------------------------------------------------------
+!-----------------------------------------------------------------------

LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/conf_gcm.F

@@ -6 +6 @@
       SUBROUTINE conf_gcm( tapedef, etatinit, clesphy0 )
 c
+#ifdef CPP_IOIPSL
       use IOIPSL
+#else
+! if not using IOIPSL, we still need to use (a local version of) getin
+      use ioipsl_getincom
+#endif
       use misc_mod
       use mod_filtre_fft, ONLY : use_filtre_fft
@@ -99 +104 @@
 
 !Config  Key  = prt_level
-!Config  Desc = niveau d'impressions de d�bogage
-!Config  Def  = 0
-!Config  Help = Niveau d'impression pour le d�bogage
+!Config  Desc = niveau d'impressions de débogage
+!Config  Def  = 0
+!Config  Help = Niveau d'impression pour le débogage
 !Config         (0 = minimum d'impression)
       prt_level = 0
@@ -109 +114 @@
 c  Parametres de controle du run:
 c-----------------------------------------------------------------------
+!Config  Key  = planet_type
+!Config  Desc = planet type ("earth", "mars", "venus", ...)
+!Config  Def  = earth
+!Config  Help = this flag sets the type of atymosphere that is considered
+      planet_type="earth"
+      CALL getin('planet_type',planet_type)
 
 !Config  Key  = dayref
@@ -189 +200 @@
        CALL getin('periodav',periodav)
 
+!Config  Key  = output_grads_dyn
+!Config  Desc = output dynamics diagnostics in 'dyn.dat' file
+!Config  Def  = n
+!Config  Help = output dynamics diagnostics in Grads-readable 'dyn.dat' file
+       output_grads_dyn=.false.
+       CALL getin('output_grads_dyn',output_grads_dyn)
+
 !Config  Key  = idissip
 !Config  Desc = periode de la dissipation 
@@ -284 +302 @@
 c    ...............................................................
 
+!Config  Key  =  read_start
+!Config  Desc = Initialize model using a 'start.nc' file
+!Config  Def  = y
+!Config  Help = y: intialize dynamical fields using a 'start.nc' file
+!               n: fields are initialized by 'iniacademic' routine
+       read_start= .true.
+       CALL getin('read_start',read_start)
+
 !Config  Key  = iflag_phys
 !Config  Desc = Avec ls physique 
@@ -341 +367 @@
 c
       IF( ABS(clat - clatt).GE. 0.001 )  THEN
-        PRINT *,' La valeur de clat passee par run.def est differente de 
-     *  celle lue sur le fichier  start '
+        write(lunout,*)'conf_gcm: La valeur de clat passee par run.def',
+     &    ' est differente de celle lue sur le fichier  start '
         STOP
       ENDIF
@@ -356 +382 @@
 
       IF( ABS(grossismx - grossismxx).GE. 0.001 )  THEN
-        PRINT *,' La valeur de grossismx passee par run.def est differente  
-     *  de celle lue sur le fichier  start '
+        write(lunout,*)'conf_gcm: La valeur de grossismx passee par ',
+     &  'run.def est differente de celle lue sur le fichier  start '
         STOP
       ENDIF
@@ -370 +396 @@
 
       IF( ABS(grossismy - grossismyy).GE. 0.001 )  THEN
-        PRINT *,' La valeur de grossismy passee par run.def est differen
-     * te de celle lue sur le fichier  start '
+        write(lunout,*)'conf_gcm: La valeur de grossismy passee par ',
+     & 'run.def est differente de celle lue sur le fichier  start '
         STOP
       ENDIF
       
       IF( grossismx.LT.1. )  THEN
-        PRINT *,' ***  ATTENTION !! grossismx < 1 .   *** '
+        write(lunout,*)
+     &       'conf_gcm: ***  ATTENTION !! grossismx < 1 .   *** '
          STOP
       ELSE
@@ -384 +411 @@
 
       IF( grossismy.LT.1. )  THEN
-        PRINT *,' ***  ATTENTION !! grossismy < 1 .   *** '
+        write(lunout,*)
+     &       'conf_gcm: ***  ATTENTION !! grossismy < 1 .   *** '
          STOP
       ELSE
@@ -390 +418 @@
       ENDIF
 
-      PRINT *,' alphax alphay defrun ',alphax,alphay
+      write(lunout,*)'conf_gcm: alphax alphay',alphax,alphay
 c
 c    alphax et alphay sont les anciennes formulat. des grossissements
@@ -405 +433 @@
 
       IF( .NOT.fxyhypb )  THEN
-           IF( fxyhypbb )     THEN
-              PRINT *,' ********  PBS DANS  DEFRUN  ******** '
-              PRINT *,' *** fxyhypb lu sur le fichier start est F ',
-     *       'alors  qu il est  T  sur  run.def  ***'
+         IF( fxyhypbb )     THEN
+            write(lunout,*)' ********  PBS DANS  CONF_GCM  ******** '
+            write(lunout,*)' *** fxyhypb lu sur le fichier start est ',
+     *       'F alors  qu il est  T  sur  run.def  ***'
               STOP
-           ENDIF
+         ENDIF
       ELSE
-           IF( .NOT.fxyhypbb )   THEN
-              PRINT *,' ********  PBS DANS  DEFRUN  ******** '
-              PRINT *,' ***  fxyhypb lu sur le fichier start est T ',
-     *        'alors  qu il est  F  sur  run.def  ****  '
+         IF( .NOT.fxyhypbb )   THEN
+            write(lunout,*)' ********  PBS DANS  CONF_GCM  ******** '
+            write(lunout,*)' ***  fxyhypb lu sur le fichier start est ',
+     *        'T alors  qu il est  F  sur  run.def  ****  '
               STOP
-           ENDIF
+         ENDIF
       ENDIF
 c
@@ -430 +458 @@
       IF( fxyhypb )  THEN
        IF( ABS(dzoomx - dzoomxx).GE. 0.001 )  THEN
-        PRINT *,' La valeur de dzoomx passee par run.def est differente 
-     *  de celle lue sur le fichier  start '
+        write(lunout,*)'conf_gcm: La valeur de dzoomx passee par ',
+     *  'run.def est differente de celle lue sur le fichier  start '
         STOP
        ENDIF
@@ -446 +474 @@
       IF( fxyhypb )  THEN
        IF( ABS(dzoomy - dzoomyy).GE. 0.001 )  THEN
-        PRINT *,' La valeur de dzoomy passee par run.def est differente 
-     * de celle lue sur le fichier  start '
+        write(lunout,*)'conf_gcm: La valeur de dzoomy passee par ',
+     * 'run.def est differente de celle lue sur le fichier  start '
         STOP
        ENDIF
@@ -461 +489 @@
       IF( fxyhypb )  THEN
        IF( ABS(taux - tauxx).GE. 0.001 )  THEN
-        PRINT *,' La valeur de taux passee par run.def est differente
-     * de celle lue sur le fichier  start '
+        write(lunout,*)'conf_gcm: La valeur de taux passee par ',
+     * 'run.def est differente de celle lue sur le fichier  start '
         STOP
        ENDIF
@@ -476 +504 @@
       IF( fxyhypb )  THEN
        IF( ABS(tauy - tauyy).GE. 0.001 )  THEN
-        PRINT *,' La valeur de tauy passee par run.def est differente
-     * de celle lue sur le fichier  start '
+        write(lunout,*)'conf_gcm: La valeur de tauy passee par ',
+     * 'run.def est differente de celle lue sur le fichier  start '
         STOP
        ENDIF
@@ -495 +523 @@
 
         IF( .NOT.ysinus )  THEN
-           IF( ysinuss )     THEN
-              PRINT *,' ********  PBS DANS  DEFRUN  ******** '
-              PRINT *,' *** ysinus lu sur le fichier start est F ',
-     *       'alors  qu il est  T  sur  run.def  ***'
+          IF( ysinuss )     THEN
+            write(lunout,*)' ********  PBS DANS  CONF_GCM  ******** '
+            write(lunout,*)' *** ysinus lu sur le fichier start est F',
+     *       ' alors  qu il est  T  sur  run.def  ***'
+            STOP
+          ENDIF
+        ELSE
+          IF( .NOT.ysinuss )   THEN
+            write(lunout,*)' ********  PBS DANS  CONF_GCM  ******** '
+            write(lunout,*)' *** ysinus lu sur le fichier start est T',
+     *        ' alors  qu il est  F  sur  run.def  ****  '
               STOP
-           ENDIF
-        ELSE
-           IF( .NOT.ysinuss )   THEN
-              PRINT *,' ********  PBS DANS  DEFRUN  ******** '
-              PRINT *,' ***  ysinus lu sur le fichier start est T ',
-     *        'alors  qu il est  F  sur  run.def  ****  '
-              STOP
-           ENDIF
+          ENDIF
         ENDIF
-      ENDIF
+      ENDIF ! of IF( .NOT.fxyhypb  )
 c
 !Config  Key  = offline
@@ -532 +560 @@
       write(lunout,*)' #########################################'
       write(lunout,*)' Configuration des parametres du gcm: '
+      write(lunout,*)' planet_type = ', planet_type
       write(lunout,*)' dayref = ', dayref
       write(lunout,*)' anneeref = ', anneeref
@@ -540 +569 @@
       write(lunout,*)' iecri = ', iecri
       write(lunout,*)' periodav = ', periodav 
+      write(lunout,*)' output_grads_dyn = ', output_grads_dyn
       write(lunout,*)' idissip = ', idissip
       write(lunout,*)' lstardis = ', lstardis
@@ -550 +580 @@
       write(lunout,*)' coefdis = ', coefdis
       write(lunout,*)' purmats = ', purmats
+      write(lunout,*)' read_start = ', read_start
       write(lunout,*)' iflag_phys = ', iflag_phys
       write(lunout,*)' clonn = ', clonn 
@@ -600 +631 @@
 
       IF( grossismx.LT.1. )  THEN
-        PRINT *,' ***  ATTENTION !! grossismx < 1 .   *** '
+        write(lunout,*)
+     &   'conf_gcm: ***  ATTENTION !! grossismx < 1 .   *** '
          STOP
       ELSE
@@ -608 +640 @@
 
       IF( grossismy.LT.1. )  THEN
-        PRINT *,' ***  ATTENTION !! grossismy < 1 .   *** '
+        write(lunout,*)
+     &  'conf_gcm: ***  ATTENTION !! grossismy < 1 .   *** '
          STOP
       ELSE
@@ -614 +647 @@
       ENDIF
 
-      PRINT *,' alphax alphay defrun ',alphax,alphay
+      write(lunout,*)'conf_gcm: alphax alphay ',alphax,alphay
 c
 c    alphax et alphay sont les anciennes formulat. des grossissements
@@ -697 +730 @@
         write(lunout,*)"Le zoom en longitude est incompatible",
      &                 " avec l'utilisation du filtre FFT ",
-     &                 "---> filtre FFT désactivé "
+     &                 "---> filtre FFT désactivé "
        use_filtre_fft=.FALSE.
       ENDIF
@@ -704 +737 @@
       
 !Config  Key  = use_mpi_alloc
-!Config  Desc = Utilise un buffer MPI en m�moire globale
+!Config  Desc = Utilise un buffer MPI en m�moire globale
 !Config  Def  = false
 !Config  Help = permet d'activer l'utilisation d'un buffer MPI
-!Config         en m�moire globale a l'aide de la fonction MPI_ALLOC.
-!Config         Cela peut am�liorer la bande passante des transferts MPI
+!Config         en m�moire globale a l'aide de la fonction MPI_ALLOC.
+!Config         Cela peut am�liorer la bande passante des transferts MPI
 !Config         d'un facteur 2  
       use_mpi_alloc=.FALSE.
@@ -716 +749 @@
 !Config  Desc = taille des blocs openmp
 !Config  Def  = 1
-!Config  Help = defini la taille des packets d'it�ration openmp
-!Config         distribu�e � chaque t�che lors de l'entr�e dans une
-!Config         boucle parall�lis�e
+!Config  Help = defini la taille des packets d'it�ration openmp
+!Config         distribu�e � chaque t�che lors de l'entr�e dans une
+!Config         boucle parall�lis�e
   
       omp_chunk=1
@@ -726 +759 @@
 !Config  Desc = activation de la version strato
 !Config  Def  = .FALSE.
-!Config  Help = active la version stratosph�rique de LMDZ de F. Lott
+!Config  Help = active la version stratosph�rique de LMDZ de F. Lott
 
       ok_strato=.FALSE.
@@ -741 +774 @@
       write(lunout,*)' #########################################'
       write(lunout,*)' Configuration des parametres du gcm: '
+      write(lunout,*)' planet_type = ', planet_type
       write(lunout,*)' dayref = ', dayref
       write(lunout,*)' anneeref = ', anneeref
@@ -749 +783 @@
       write(lunout,*)' iecri = ', iecri
       write(lunout,*)' periodav = ', periodav 
+      write(lunout,*)' output_grads_dyn = ', output_grads_dyn
       write(lunout,*)' idissip = ', idissip
       write(lunout,*)' lstardis = ', lstardis
@@ -759 +794 @@
       write(lunout,*)' coefdis = ', coefdis
       write(lunout,*)' purmats = ', purmats
+      write(lunout,*)' read_start = ', read_start
       write(lunout,*)' iflag_phys = ', iflag_phys
       write(lunout,*)' clon = ', clon 
@@ -764 +800 @@
       write(lunout,*)' grossismx = ', grossismx
       write(lunout,*)' grossismy = ', grossismy
-      write(lunout,*)' fxyhypbb = ', fxyhypbb
+      write(lunout,*)' fxyhypb = ', fxyhypb
       write(lunout,*)' dzoomx = ', dzoomx
       write(lunout,*)' dzoomy = ', dzoomy

LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/control.h

@@ -14 +14 @@
      &              iperiod,iapp_tracvl,iconser,iecri,idissip,iphysiq , &
      &              periodav,iecrimoy,dayref,anneeref,                  &
-     &              raz_date,offline,ip_ebil_dyn,config_inca
+     &              raz_date,offline,ip_ebil_dyn,config_inca,           &
+     &              planet_type,output_grads_dyn
 
       INTEGER   nday,day_step,iperiod,iapp_tracvl,iconser,iecri,        &
@@ -21 +22 @@
       REAL periodav
       logical offline
-      CHARACTER*4 config_inca
+      CHARACTER (len=4) :: config_inca
+      CHARACTER(len=10) :: planet_type ! planet type ('earth','mars',...)
+      LOGICAL :: output_grads_dyn ! output dynamics diagnostics in
+                                  ! binary grads file 'dyn.dat' (y/n)
 !-----------------------------------------------------------------------

LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/diagedyn.F

@@ -58 +58 @@
 #include "paramet.h"
 #include "comgeom.h"
-
-#ifdef CPP_PHYS
+#include "iniprint.h"
+
+#ifdef CPP_EARTH
 #include "../phylmd/YOMCST.h"
 #include "../phylmd/YOETHF.h"
@@ -139 +140 @@
 
 
-#ifdef CPP_PHYS
+#ifdef CPP_EARTH
 c======================================================================
 C     Compute Kinetic enrgy
@@ -314 +315 @@
 C
 #else
-      print*,'Pour l instant diagedyn a besoin de la physique'
+      write(lunout,*),'diagedyn: Needs Earth physics to function'
 #endif
+! #endif of #ifdef CPP_EARTH 
       RETURN 
       END 

LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/etat0_netcdf.F

@@ -5 +5 @@
 c
       SUBROUTINE etat0_netcdf (interbar, masque)
-    
+#ifdef CPP_EARTH        
       USE startvar
       USE ioipsl
@@ -14 +14 @@
       USE filtreg_mod
       USE infotrac
+#endif
+!#endif of #ifdef CPP_EARTH
       !
       IMPLICIT NONE
@@ -25 +27 @@
 !     .KLON=KFDIA-KIDIA+1,KLEV=llm
       !
+#ifdef CPP_EARTH    
 #include "comgeom2.h"
 #include "comvert.h"
@@ -31 +34 @@
 #include "dimsoil.h"
 #include "temps.h"
-      !
+#endif
+!#endif of #ifdef CPP_EARTH
+      ! arguments:
       LOGICAL interbar
+      REAL :: masque(iip1,jjp1)
+
+#ifdef CPP_EARTH
+      ! local variables:
       REAL :: latfi(klon), lonfi(klon)
-      REAL :: orog(iip1,jjp1), rugo(iip1,jjp1), masque(iip1,jjp1),
+      REAL :: orog(iip1,jjp1), rugo(iip1,jjp1),
      . psol(iip1, jjp1), phis(iip1, jjp1)
       REAL :: p3d(iip1, jjp1, llm+1)
@@ -144 +153 @@
       !
       preff     = 101325.
+      pa        =  50000.
       unskap = 1./kappa
       !
@@ -167 +177 @@
       print*,'dtvr',dtvr
 
-      CALL inicons0()
+      CALL iniconst()
       CALL inigeom()
       !
@@ -754 +764 @@
       print*,'entree histclo'
       CALL histclo
+
+#endif 
+!#endif of #ifdef CPP_EARTH
       RETURN
       !

LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/gcm.F

@@ -9 +9 @@
       USE IOIPSL
 #endif
+
       USE mod_const_mpi, ONLY: init_const_mpi
       USE parallel
       USE mod_phys_lmdz_para, ONLY : klon_mpi_para_nb
+      USE infotrac
+      USE mod_interface_dyn_phys
+      USE mod_hallo
+      USE Bands
+
+      USE filtreg_mod
+
+! Ehouarn: for now these only apply to Earth:
+#ifdef CPP_EARTH
       USE mod_grid_phy_lmdz
       USE mod_phys_lmdz_omp_data, ONLY: klon_omp 
       USE dimphy
-      USE infotrac
-      USE mod_interface_dyn_phys
       USE comgeomphy
-      USE mod_hallo
-      USE Bands
-
-      USE filtreg_mod
-
+#endif
       IMPLICIT NONE
 
@@ -160 +164 @@
       dynhistave_file = 'dyn_hist_ave'
 
-c--------------------------------------------------------------------------
-c   Iflag_phys controle l'appel a la physique :
-c   -------------------------------------------
-c      0 : pas de physique
-c      1 : Normale (appel a phylmd, phymars ...)
-c      2 : rappel Newtonien pour la temperature + friction au sol
-      iflag_phys=1
-
-c--------------------------------------------------------------------------
-c   Lecture de l'etat initial :
-c   ---------------------------
-c     T : on lit start.nc
-c     F : le modele s'autoinitialise avec un cas academique (iniacademic)
-#ifdef CPP_IOIPSL
-      read_start=.true.
-#else
-      read_start=.false.
-#endif
 
 c-----------------------------------------------------------------------
@@ -194 +180 @@
 c  ---------------------------------------
 c
-#ifdef CPP_IOIPSL
+! Ehouarn: dump possibility of using defrun
+!#ifdef CPP_IOIPSL
       CALL conf_gcm( 99, .TRUE. , clesphy0 )
-#else
-      CALL defrun( 99, .TRUE. , clesphy0 )
-#endif
+!#else
+!      CALL defrun( 99, .TRUE. , clesphy0 )
+!#endif
 c
 c
@@ -208 +195 @@
       call init_parallel
       call Read_Distrib
-      CALL Init_Phys_lmdz(iim,jjp1,llm,mpi_size,distrib_phys)
+! Ehouarn : temporarily (?) keep this only for Earth
+      if (planet_type.eq."earth") then
+#ifdef CPP_EARTH
+        CALL Init_Phys_lmdz(iim,jjp1,llm,mpi_size,distrib_phys)
+#endif
+      endif ! of if (planet_type.eq."earth")
       CALL set_bands
       CALL Init_interface_dyn_phys
@@ -220 +212 @@
 c$OMP END PARALLEL
 
+! Ehouarn : temporarily (?) keep this only for Earth
+      if (planet_type.eq."earth") then
+#ifdef CPP_EARTH
 c$OMP PARALLEL
       call InitComgeomphy
 c$OMP END PARALLEL 
+#endif
+      endif ! of if (planet_type.eq."earth")
 
       IF (config_inca /= 'none') THEN
@@ -252 +249 @@
 c  lecture du fichier start.nc
       if (read_start) then
-#ifdef CPP_IOIPSL
+      ! we still need to run iniacademic to initialize some
+      ! constants & fields, if we run the 'newtonian' case:
+        if (iflag_phys.eq.2) then
+          CALL iniacademic(vcov,ucov,teta,q,masse,ps,phis,time_0)
+        endif
+!#ifdef CPP_IOIPSL
+        if (planet_type.eq."earth") then
+#ifdef CPP_EARTH
+! Load an Earth-format start file
          CALL dynetat0("start.nc",vcov,ucov,
      .              teta,q,masse,ps,phis, time_0)
+#endif
+        endif ! of if (planet_type.eq."earth")
 c       write(73,*) 'ucov',ucov
 c       write(74,*) 'vcov',vcov
@@ -261 +268 @@
 c       write(77,*) 'q',q
 
-#endif
-      endif
+      endif ! of if (read_start)
 
 c le cas echeant, creation d un etat initial
       IF (prt_level > 9) WRITE(lunout,*)
-     .                 'AVANT iniacademic AVANT AVANT AVANT AVANT'
+     .              'GCM: AVANT iniacademic AVANT AVANT AVANT AVANT'
       if (.not.read_start) then
          CALL iniacademic(vcov,ucov,teta,q,masse,ps,phis,time_0)
@@ -351 +357 @@
 c   Initialisation de la physique :
 c   -------------------------------
-#ifdef CPP_PHYS
       IF (call_iniphys.and.iflag_phys.eq.1) THEN
          latfi(1)=rlatu(1)
@@ -372 +377 @@
 
          WRITE(lunout,*)
-     .           'WARNING!!! vitesse verticale nulle dans la physique'
-
+     .       'GCM: WARNING!!! vitesse verticale nulle dans la physique'
+! Earth:
+         if (planet_type.eq."earth") then
+#ifdef CPP_EARTH
          CALL iniphysiq(ngridmx,llm,daysec,day_ini,dtphys ,
      ,                latfi,lonfi,airefi,zcufi,zcvfi,rad,g,r,cpp     )
-
+#endif
+         endif ! of if (planet_type.eq."earth")
          call_iniphys=.false.
-
-      ENDIF
-#endif
+      ENDIF ! of IF (call_iniphys.and.(iflag_phys.eq.1))
 
 
@@ -402 +408 @@
       day_end = day_ini + nday
       WRITE(lunout,300)day_ini,day_end
+ 300  FORMAT('1'/,15x,'run du jour',i7,2x,'au jour',i7//)
+
+!#ifdef CPP_IOIPSL
+      if (planet_type.eq."earth") then
+#ifdef CPP_EARTH
+      CALL dynredem0_p("restart.nc", day_end, phis)
+#endif
+      endif
+
+      ecripar = .TRUE.
 
 #ifdef CPP_IOIPSL
-      CALL dynredem0_p("restart.nc", day_end, phis)
-
-      ecripar = .TRUE.
-
       if ( 1.eq.1) then
       time_step = zdtvr
@@ -425 +437 @@
 
 #endif
+! #endif of #ifdef CPP_IOIPSL
 
 c  Choix des frequences de stokage pour le offline
@@ -450 +463 @@
 
 
- 300  FORMAT('1'/,15x,'run du pas',i7,2x,'au pas',i7,2x,
-     . 'c''est a dire du jour',i7,3x,'au jour',i7//)
       END
 

LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/iniacademic.F

@@ -45 +45 @@
 #include "temps.h"
 #include "control.h"
+#include "iniprint.h"
 
 c   Arguments:
@@ -57 +58 @@
       REAL ps(ip1jmp1)                       ! pression  au sol
       REAL masse(ip1jmp1,llm)                ! masse d'air
+      REAL phis(ip1jmp1)                     ! geopotentiel au sol
+
+c   Local:
+c   ------
+
       REAL p (ip1jmp1,llmp1  )               ! pression aux interfac.des couches
       REAL pks(ip1jmp1)                      ! exner au  sol
       REAL pk(ip1jmp1,llm)                   ! exner au milieu des couches
       REAL pkf(ip1jmp1,llm)                  ! exner filt.au milieu des couches
-      REAL phis(ip1jmp1)                     ! geopotentiel au sol
       REAL phi(ip1jmp1,llm)                  ! geopotentiel
-
-
-
-
-
-c   Local:
-c   ------
-
       REAL ddsin,tetarappelj,tetarappell,zsig
       real tetajl(jjp1,llm)
@@ -81 +78 @@
 
 c-----------------------------------------------------------------------
+! 1. Initializations for Earth-like case
+! --------------------------------------
+      if (planet_type=="earth") then
+c
+        time_0=0.
 
-c
-      time_0=0.
+        im         = iim
+        jm         = jjm
+        day_ini    = 0
+        omeg       = 4.*asin(1.)/86400.
+        rad    = 6371229.
+        g      = 9.8
+        daysec = 86400.
+        dtvr    = daysec/FLOAT(day_step)
+        zdtvr=dtvr
+        kappa  = 0.2857143
+        cpp    = 1004.70885
+        preff     = 101325.
+        pa        =  50000.
+        etot0      = 0.
+        ptot0      = 0.
+        ztot0      = 0.
+        stot0      = 0.
+        ang0       = 0.
 
-      im         = iim
-      jm         = jjm
-      day_ini    = 0
-      omeg       = 4.*asin(1.)/86400.
-      rad    = 6371229.
-      g      = 9.8
-      daysec = 86400.
-      dtvr    = daysec/FLOAT(day_step)
-      zdtvr=dtvr
-      kappa  = 0.2857143
-      cpp    = 1004.70885
-      preff     = 101325.
-      pa        =  50 000.
-      etot0      = 0.
-      ptot0      = 0.
-      ztot0      = 0.
-      stot0      = 0.
-      ang0       = 0.
-      pa         = 0.
+        CALL iniconst
+        CALL inigeom
+        CALL inifilr
 
-      CALL inicons0
-      CALL inigeom
-      CALL inifilr
-
-      ps=0.
-      phis=0.
+        ps=0.
+        phis=0.
 c---------------------------------------------------------------------
 
-      taurappel=10.*daysec
+        taurappel=10.*daysec
 
 c---------------------------------------------------------------------
@@ -119 +117 @@
 c   --------------------------------------
 
-      DO l=1,llm
-       zsig=ap(l)/preff+bp(l)
-       if (zsig.gt.0.3) then
-         lsup=l
-         tetarappell=1./8.*(-log(zsig)-.5)
-         DO j=1,jjp1
-            ddsin=sin(rlatu(j))-sin(pi/20.)
-            tetajl(j,l)=300.*(1+1./18.*(1.-3.*ddsin*ddsin)+tetarappell)
-         ENDDO
-        else
+        DO l=1,llm
+         zsig=ap(l)/preff+bp(l)
+         if (zsig.gt.0.3) then
+           lsup=l
+           tetarappell=1./8.*(-log(zsig)-.5)
+           DO j=1,jjp1
+             ddsin=sin(rlatu(j))-sin(pi/20.)
+             tetajl(j,l)=300.*(1+1./18.*(1.-3.*ddsin*ddsin)+tetarappell)
+           ENDDO
+          else
 c   Choix isotherme au-dessus de 300 mbar
-         do j=1,jjp1
-            tetajl(j,l)=tetajl(j,lsup)*(0.3/zsig)**kappa
-         enddo
-        endif
-      ENDDO
+           do j=1,jjp1
+             tetajl(j,l)=tetajl(j,lsup)*(0.3/zsig)**kappa
+           enddo
+          endif ! of if (zsig.gt.0.3)
+        ENDDO ! of DO l=1,llm
 
-      do l=1,llm
-         do j=1,jjp1
-            do i=1,iip1
-               ij=(j-1)*iip1+i
-               tetarappel(ij,l)=tetajl(j,l)
-            enddo
-         enddo
-      enddo
+        do l=1,llm
+           do j=1,jjp1
+              do i=1,iip1
+                 ij=(j-1)*iip1+i
+                 tetarappel(ij,l)=tetajl(j,l)
+              enddo
+           enddo
+        enddo
 
-c     call dump2d(jjp1,llm,tetajl,'TEQ   ')
+c       call dump2d(jjp1,llm,tetajl,'TEQ   ')
 
-      ps=1.e5
-      phis=0.
-      CALL pression ( ip1jmp1, ap, bp, ps, p       )
-      CALL exner_hyb( ip1jmp1, ps, p,alpha,beta, pks, pk, pkf )
-      CALL massdair(p,masse)
+        ps=1.e5
+        phis=0.
+        CALL pression ( ip1jmp1, ap, bp, ps, p       )
+        CALL exner_hyb( ip1jmp1, ps, p,alpha,beta, pks, pk, pkf )
+        CALL massdair(p,masse)
 
 c  intialisation du vent et de la temperature
-      teta(:,:)=tetarappel(:,:)
-      CALL geopot(ip1jmp1,teta,pk,pks,phis,phi)
-      call ugeostr(phi,ucov)
-      vcov=0.
-      q(:,:,1   )=1.e-10
-      q(:,:,2   )=1.e-15
-      q(:,:,3:nqtot)=0.
+        teta(:,:)=tetarappel(:,:)
+        CALL geopot(ip1jmp1,teta,pk,pks,phis,phi)
+        call ugeostr(phi,ucov)
+        vcov=0.
+        q(:,:,1   )=1.e-10
+        q(:,:,2   )=1.e-15
+        q(:,:,3:nqtot)=0.
 
 
-c   perturbation al\351atoire sur la temp\351rature
-      idum  = -1
-      zz = ran1(idum)
-      idum  = 0
-      do l=1,llm
-         do ij=iip2,ip1jm
-            teta(ij,l)=teta(ij,l)*(1.+0.005*ran1(idum))
-         enddo
-      enddo
+c   perturbation aleatoire sur la temperature
+        idum  = -1
+        zz = ran1(idum)
+        idum  = 0
+        do l=1,llm
+           do ij=iip2,ip1jm
+              teta(ij,l)=teta(ij,l)*(1.+0.005*ran1(idum))
+           enddo
+        enddo
 
-      do l=1,llm
-         do ij=1,ip1jmp1,iip1
-            teta(ij+iim,l)=teta(ij,l)
-         enddo
-      enddo
+        do l=1,llm
+           do ij=1,ip1jmp1,iip1
+              teta(ij+iim,l)=teta(ij,l)
+           enddo
+        enddo
 
 
@@ -187 +185 @@
 
 c   initialisation d'un traceur sur une colonne
-      j=jjp1*3/4
-      i=iip1/2
-      ij=(j-1)*iip1+i
-      q(ij,:,3)=1.
-
+        j=jjp1*3/4
+        i=iip1/2
+        ij=(j-1)*iip1+i
+        q(ij,:,3)=1.
+      
+      else
+        write(lunout,*)"iniacademic: planet types other than earth",
+     &                 " not implemented (yet)."
+        stop
+      endif ! of if (planet_type=="earth")
       return
       END

LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/leapfrog_p.F

@@ -4 +4 @@
 c
 c
-#define IO_DEBUG
-
-#undef CPP_IOIPSL
-#define CPP_IOIPSL
 
       SUBROUTINE leapfrog_p(ucov,vcov,teta,ps,masse,phis,q,clesphy0,
@@ -285 +281 @@
 c$OMP BARRIER
        else
-          
+! Save fields obtained at previous time step as '...m1'
          ijb=ij_begin
          ije=ij_end
@@ -312 +308 @@
      .                    llm, -2,2, .TRUE., 1 )
 
-       endif
+       endif ! of if (FirstCaldyn)
        
       forward = .TRUE.
@@ -356 +352 @@
          IF( MOD(itau,idissip ).EQ.0.AND..NOT.forward ) apdiss = .TRUE.
          IF( MOD(itau,iphysiq ).EQ.0.AND..NOT.forward 
-     s          .and. iflag_phys.NE.0                 ) apphys = .TRUE.
+     s          .and. iflag_phys.EQ.1                 ) apphys = .TRUE.
       ELSE
          IF( MOD(itau   ,iconser) .EQ. 0              ) conser = .TRUE.
          IF( MOD(itau+1,idissip)  .EQ. 0              ) apdiss = .TRUE.
-         IF( MOD(itau+1,iphysiq).EQ.0.AND.iflag_phys.NE.0) apphys=.TRUE.
+         IF( MOD(itau+1,iphysiq).EQ.0.AND.iflag_phys.EQ.1) apphys=.TRUE.
       END IF
 
@@ -537 +533 @@
 
 c$OMP MASTER
-      IF (prt_level>9)WRITE(lunout,*)"Iteration No",True_itau
+      IF (prt_level>9) THEN
+        WRITE(lunout,*)"leapfrog_p: Iteration No",True_itau
+      ENDIF
 
 
@@ -594 +592 @@
 
 
-         ENDIF
-c
-      ENDIF
+         ENDIF ! of IF (offline)
+c
+      ENDIF ! of IF( forward. OR . leapf )
 cc$OMP END PARALLEL
 
@@ -673 +671 @@
 c$OMP MASTER
          call suspend_timer(timer_caldyn)
-         WRITE(lunout,*)'Entree dans la physique : Iteration No ',      &
-     &                   true_itau
+
+         write(lunout,*)
+     &   'leapfrog_p: Entree dans la physique : Iteration No ',true_itau
 c$OMP END MASTER
 
@@ -693 +692 @@
 c   -----------------------------------------------------
 
-#ifdef CPP_PHYS
 c+jld
 
@@ -699 +697 @@
       IF (ip_ebil_dyn.ge.1 ) THEN 
           ztit='bil dyn'
-          CALL diagedyn(ztit,2,1,1,dtphys
-     e  , ucov    , vcov , ps, p ,pk , teta , q(:,:,1), q(:,:,2))
+! Ehouarn: be careful, diagedyn is Earth-specific (includes ../phylmd/..)!
+           IF (planet_type.eq."earth") THEN
+            CALL diagedyn(ztit,2,1,1,dtphys
+     &    , ucov    , vcov , ps, p ,pk , teta , q(:,:,1), q(:,:,2))
+           ENDIF
       ENDIF 
 c-jld
@@ -787 +788 @@
           dpfi_tmp(1:iip1)     = dpfi(ijb:ijb+iim)  
 c$OMP END MASTER
-        endif
+        endif ! of if ( .not. pole_nord)
 
 c$OMP BARRIER
@@ -843 +844 @@
 c$OMP END MASTER
           
-        endif
+        endif ! of if (.not. pole_nord)
 c$OMP BARRIER
 cc$OMP MASTER        
@@ -932 +933 @@
 cc$OMP END MASTER
 
-#else
-
+
+c-jld
+c$OMP MASTER
+         call resume_timer(timer_caldyn)
+         if (FirstPhysic) then
+           ok_start_timer=.TRUE.
+           FirstPhysic=.false.
+         endif
+c$OMP END MASTER
+       ENDIF ! of IF( apphys )
+
+      IF(iflag_phys.EQ.2) THEN ! "Newtonian" case
 c   Calcul academique de la physique = Rappel Newtonien + fritcion 
 c   --------------------------------------------------------------
@@ -949 +960 @@
 
        call friction_p(ucov,vcov,iphysiq*dtvr)
-
-#endif
-
-c-jld
-c$OMP MASTER
-         call resume_timer(timer_caldyn)
-         if (FirstPhysic) then
-           ok_start_timer=.TRUE.
-           FirstPhysic=.false.
-         endif
-c$OMP END MASTER
-       ENDIF
+      ENDIF ! of IF(iflag_phys.EQ.2)
+
 
         CALL pression_p ( ip1jmp1, ap, bp, ps, p                  )
@@ -1312 +1313 @@
 c      IF( MOD(itau,iecri         ).EQ.0) THEN
 
-           IF( MOD(itau,iecri*day_step).EQ.0) THEN
-c$OMP BARRIER
-c$OMP MASTER
-               nbetat = nbetatdem
-        CALL geopot_p  ( ip1jmp1, teta  , pk , pks,  phis  , phi )
+            IF( MOD(itau,iecri*day_step).EQ.0) THEN
+c$OMP BARRIER
+c$OMP MASTER
+              nbetat = nbetatdem
+              CALL geopot_p(ip1jmp1,teta,pk,pks,phis,phi)
        
 cym        unat=0.
         
-        ijb=ij_begin
-        ije=ij_end
-        
-        if (pole_nord) then
-          ijb=ij_begin+iip1
-          unat(1:iip1,:)=0.
-        endif
-        
-        if (pole_sud) then 
-          ije=ij_end-iip1
-          unat(ij_end-iip1+1:ij_end,:)=0.
-        endif
+              ijb=ij_begin
+              ije=ij_end
+        
+              if (pole_nord) then
+                ijb=ij_begin+iip1
+                unat(1:iip1,:)=0.
+              endif
+        
+              if (pole_sud) then 
+                ije=ij_end-iip1
+                unat(ij_end-iip1+1:ij_end,:)=0.
+              endif
             
-        do l=1,llm
-           unat(ijb:ije,l)=ucov(ijb:ije,l)/cu(ijb:ije)
-        enddo
-
-        ijb=ij_begin
-        ije=ij_end
-        if (pole_sud) ije=ij_end-iip1
-        
-        do l=1,llm
-           vnat(ijb:ije,l)=vcov(ijb:ije,l)/cv(ijb:ije)
-        enddo
+              do l=1,llm
+                unat(ijb:ije,l)=ucov(ijb:ije,l)/cu(ijb:ije)
+              enddo
+
+              ijb=ij_begin
+              ije=ij_end
+              if (pole_sud) ije=ij_end-iip1
+        
+              do l=1,llm
+                vnat(ijb:ije,l)=vcov(ijb:ije,l)/cv(ijb:ije)
+              enddo
         
 #ifdef CPP_IOIPSL
  
-        CALL writehist_p(histid,histvid, itau,vcov, 
-     s                       ucov,teta,phi,q,masse,ps,phis)
+              CALL writehist_p(histid,histvid, itau,vcov, 
+     &                         ucov,teta,phi,q,masse,ps,phis)
 
 #endif
-c$OMP END MASTER
-           ENDIF
+! For some Grads outputs of fields
+              if (output_grads_dyn) then
+! Ehouarn: hope this works the way I think it does:
+                  call Gather_Field(unat,ip1jmp1,llm,0)
+                  call Gather_Field(vnat,ip1jm,llm,0)
+                  call Gather_Field(teta,ip1jmp1,llm,0)
+                  call Gather_Field(ps,ip1jmp1,1,0)
+                  do iq=1,nqtot
+                    call Gather_Field(q(1,1,iq),ip1jmp1,llm,0)
+                  enddo
+                  if (mpi_rank==0) then
+#include "write_grads_dyn.h"
+                  endif
+              endif ! of if (output_grads_dyn)
+c$OMP END MASTER
+            ENDIF ! of IF(MOD(itau,iecri).EQ.0)
 
             IF(itau.EQ.itaufin) THEN
@@ -1359 +1374 @@
 c$OMP MASTER
 
-c#ifdef CPP_IOIPSL
-
-       CALL dynredem1_p("restart.nc",0.0,
-     ,                     vcov,ucov,teta,q,masse,ps)
-c#endif
+              if (planet_type.eq."earth") then
+#ifdef CPP_EARTH
+! Write an Earth-format restart file
+                CALL dynredem1_p("restart.nc",0.0,
+     &                           vcov,ucov,teta,q,masse,ps)
+
+#endif
+              endif ! of if (planet_type.eq."earth")
 
               CLOSE(99)
 c$OMP END MASTER
-            ENDIF
+            ENDIF ! of IF (itau.EQ.itaufin)
 
 c-----------------------------------------------------------------------
@@ -1398 +1416 @@
                  dt  = 2.*dtvr
                  GO TO 2
-            END IF
-
-      ELSE
+            END IF ! of IF (MOD(itau,iperiod).EQ.0)
+                   !    ELSEIF (MOD(itau-1,iperiod).EQ.0)
+
+
+      ELSE ! of IF (.not.purmats)
 
 c       ........................................................
@@ -1427 +1447 @@
                GO TO 2
 
-            ELSE
-
-            IF(MOD(itau,iperiod).EQ.0 .OR. itau.EQ.itaufin) THEN
+            ELSE ! of IF(forward)
+
+              IF(MOD(itau,iperiod).EQ.0 .OR. itau.EQ.itaufin) THEN
                IF(itau.EQ.itaufin) THEN
                   iav=1
@@ -1438 +1458 @@
 c$OMP BARRIER
 
-              call Register_Hallo(vcov,ip1jm,llm,1,0,0,1,TestRequest)
-              call SendRequest(TestRequest)
-c$OMP BARRIER
-              call WaitRequest(TestRequest)
-
-c$OMP BARRIER
-c$OMP MASTER
-              CALL writedynav_p(histaveid, itau,vcov ,
+               call Register_Hallo(vcov,ip1jm,llm,1,0,0,1,TestRequest)
+               call SendRequest(TestRequest)
+c$OMP BARRIER
+               call WaitRequest(TestRequest)
+
+c$OMP BARRIER
+c$OMP MASTER
+               CALL writedynav_p(histaveid, itau,vcov ,
      ,                          ucov,teta,pk,phi,q,masse,ps,phis)
                call bilan_dyn_p (2,dtvr*iperiod,dtvr*day_step*periodav,
@@ -1451 +1471 @@
 c$OMP END MASTER
 #endif
-            ENDIF
+              ENDIF ! of IF(MOD(itau,iperiod).EQ.0 .OR. itau.EQ.itaufin)
+
 
 c               IF(MOD(itau,iecri         ).EQ.0) THEN
@@ -1457 +1478 @@
 c$OMP BARRIER
 c$OMP MASTER
-                  nbetat = nbetatdem
-       CALL geopot_p( ip1jmp1, teta  , pk , pks,  phis  , phi   )
+                nbetat = nbetatdem
+                CALL geopot_p(ip1jmp1,teta,pk,pks,phis,phi)
 
 cym        unat=0.
-        ijb=ij_begin
-        ije=ij_end
-        
-        if (pole_nord) then
-          ijb=ij_begin+iip1
-          unat(1:iip1,:)=0.
-        endif
-        
-        if (pole_sud) then 
-          ije=ij_end-iip1
-          unat(ij_end-iip1+1:ij_end,:)=0.
-        endif
+                ijb=ij_begin
+                ije=ij_end
+        
+                if (pole_nord) then
+                  ijb=ij_begin+iip1
+                  unat(1:iip1,:)=0.
+                endif
+        
+                if (pole_sud) then 
+                  ije=ij_end-iip1
+                  unat(ij_end-iip1+1:ij_end,:)=0.
+                endif
             
-        do l=1,llm
-           unat(ijb:ije,l)=ucov(ijb:ije,l)/cu(ijb:ije)
-        enddo
-
-        ijb=ij_begin
-        ije=ij_end
-        if (pole_sud) ije=ij_end-iip1
-        
-        do l=1,llm
-           vnat(ijb:ije,l)=vcov(ijb:ije,l)/cv(ijb:ije)
-        enddo
+                do l=1,llm
+                  unat(ijb:ije,l)=ucov(ijb:ije,l)/cu(ijb:ije)
+                enddo
+
+                ijb=ij_begin
+                ije=ij_end
+                if (pole_sud) ije=ij_end-iip1
+        
+                do l=1,llm
+                  vnat(ijb:ije,l)=vcov(ijb:ije,l)/cv(ijb:ije)
+                enddo
 
 #ifdef CPP_IOIPSL
 
-       CALL writehist_p( histid, histvid, itau,vcov , 
-     ,                           ucov,teta,phi,q,masse,ps,phis)
-c#else
-c      call Gather_Field(unat,ip1jmp1,llm,0)
-c      call Gather_Field(vnat,ip1jm,llm,0)
-c      call Gather_Field(teta,ip1jmp1,llm,0)
-c      call Gather_Field(ps,ip1jmp1,1,0)
-c      do iq=1,nqtot
-c        call Gather_Field(q(1,1,iq),ip1jmp1,llm,0)
-c      enddo
+                CALL writehist_p(histid, histvid, itau,vcov , 
+     &                           ucov,teta,phi,q,masse,ps,phis)
+#endif
+! For some Grads output (but does it work?)
+                if (output_grads_dyn) then
+                  call Gather_Field(unat,ip1jmp1,llm,0)
+                  call Gather_Field(vnat,ip1jm,llm,0)
+                  call Gather_Field(teta,ip1jmp1,llm,0)
+                  call Gather_Field(ps,ip1jmp1,1,0)
+                  do iq=1,nqtot
+                    call Gather_Field(q(1,1,iq),ip1jmp1,llm,0)
+                  enddo
 c      
-c      if (mpi_rank==0) then
-c#include "write_grads_dyn.h"
-c      endif
-#endif
-
-c$OMP END MASTER
-               ENDIF
-
-                 IF(itau.EQ.itaufin) THEN
+                  if (mpi_rank==0) then
+#include "write_grads_dyn.h"
+                  endif
+                endif ! of if (output_grads_dyn)
+
+c$OMP END MASTER
+              ENDIF ! of IF(MOD(itau,iecri*day_step).EQ.0)
+
+              IF(itau.EQ.itaufin) THEN
+                if (planet_type.eq."earth") then
+#ifdef CPP_EARTH
 c$OMP MASTER
                    CALL dynredem1_p("restart.nc",0.0,
      .                               vcov,ucov,teta,q,masse,ps)
 c$OMP END MASTER
-                 ENDIF
-                 forward = .TRUE.
-                 GO TO  1
-
-            ENDIF
-
-      END IF
-c$OMP MASTER
-        call finalize_parallel
-c$OMP END MASTER
-        RETURN
+#endif
+                endif ! of if (planet_type.eq."earth")
+              ENDIF ! of IF(itau.EQ.itaufin)
+
+              forward = .TRUE.
+              GO TO  1
+
+            ENDIF ! of IF (forward)
+
+      END IF ! of IF(.not.purmats)
+c$OMP MASTER
+      call finalize_parallel
+c$OMP END MASTER
+      RETURN
       END

LMDZ4/branches/LMDZ4-dev/makegcm

@@ -206 +206 @@
    set optim90="-Wp'-P' -DCRAY -p$IOIPSLDIR "'-p$(LIBO) -eiv '
    set optim="$optim90"
-   set oplink="-Wl'-DSTACK=128 -f indef' -L$IOIPSLDIR -lioipsl  -L$NCDFLIB -lnetcdf "
+   if ( $io == "ioipsl" ) then
+     set oplink="-Wl'-DSTACK=128 -f indef' -L$IOIPSLDIR -lioipsl  -L$NCDFLIB -lnetcdf "
+   else
+     set oplink="-Wl'-DSTACK=128 -f indef' -L$IOIPSLDIR -L$NCDFLIB -lnetcdf "
+   endif
    set mod_loc_dir=" "
    set mod_suffix=" "
@@ -216 +220 @@
    set optim90=" -fast -fixed "
    set optimtru90=" -fast -free "
-   set opt_link="-lf77compat -L$MODIPSLDIR -lsechiba -lparameters -lstomate -lioipsl -L$NCDFLIB -lnetcdf "
+   if ( $io == "ioipsl" ) then
+     set opt_link="-lf77compat -L$MODIPSLDIR -lsechiba -lparameters -lstomate -lioipsl -L$NCDFLIB -lnetcdf "
+   else
+     set opt_link="-lf77compat -L$MODIPSLDIR -lsechiba -lparameters -lstomate -L$NCDFLIB -lnetcdf "
+   endif
    set mod_loc_dir=$localdir
    set mod_suffix=mod
@@ -270 +278 @@
    set optim="$optim90"
    set optimbis=" "
-   set opt_link=" -C hopt -float0 -ew -P static -L$MODIPSLDIR -lioipsl  $NCDFLIB -lnetcdf_i8r8_v "
+   if ( $io == "ioipsl" ) then
+     set opt_link=" -C hopt -float0 -ew -P static -L$MODIPSLDIR -lioipsl  $NCDFLIB -lnetcdf_i8r8_v "
+   else
+     set opt_link=" -C hopt -float0 -ew -P static -L$MODIPSLDIR $NCDFLIB -lnetcdf_i8r8_v "
+   endif
    set mod_loc_dir="."
    set mod_suffix="mod"
@@ -533 +545 @@
    set phys="L_PHY= LIBPHY="
 else
-   set cppflags="$cppflags -DCPP_PHYS"
+   #Default planet type is Earth
+   set cppflags="$cppflags -DCPP_EARTH"
 endif
 
@@ -685 +698 @@
 ######################################################################
 
-if ( -f $libf/phy$physique/raddim.h ) then
- if ( -f $libf/phy$physique/raddim.$dimh.h ) then
-  \rm -f $libf/phy$physique/raddim.h
-  cp -p $libf/phy$physique/raddim.$dimh.h $libf/phy$physique/raddim.h
-  echo $libf/phy$physique/raddim.$dimh.h 
-  cat $libf/phy$physique/raddim.$dimh.h 
-  cat $libf/phy$physique/raddim.h
- else
-  echo On peut diminuer la taille de l executable en creant
-  echo le fichier $libf/phy$physique/raddim.$dimh.h
-  \cp -p $libf/phy$physique/raddim.defaut.h $libf/phy$physique/raddim.h
- endif
-endif
+#if ( -f $libf/phy$physique/raddim.h ) then
+# if ( -f $libf/phy$physique/raddim.$dimh.h ) then
+#  \rm -f $libf/phy$physique/raddim.h
+#  cp -p $libf/phy$physique/raddim.$dimh.h $libf/phy$physique/raddim.h
+#  echo $libf/phy$physique/raddim.$dimh.h 
+#  cat $libf/phy$physique/raddim.$dimh.h 
+#  cat $libf/phy$physique/raddim.h
+# else
+#  echo On peut diminuer la taille de l executable en creant
+#  echo le fichier $libf/phy$physique/raddim.$dimh.h
+#  \cp -p $libf/phy$physique/raddim.defaut.h $libf/phy$physique/raddim.h
+# endif
+#endif
 
 ######################################################################
@@ -793 +806 @@
    set f90=$FC_LINUX
    if ( $FC_LINUX == 'pgf90' ) then
-   set opt_link=" -L$MODIPSLDIR $link_veget -L$NCDFLIB -lioipsl -lnetcdf "
+     if ( $io == "ioipsl" ) then
+       set opt_link=" -L$MODIPSLDIR $link_veget -L$NCDFLIB -lioipsl -lnetcdf "
+     else
+       set opt_link=" -L$MODIPSLDIR $link_veget -L$NCDFLIB -lnetcdf "
+     endif
    else if ($FC_LINUX == 'g95') then
-     set opt_link="-L$MODIPSLDIR $link_veget -lioipsl -L$NCDFLIB -lnetcdf -lioipsl -lnetcdf "
+     if ( $io == "ioipsl" ) then
+       set opt_link="-L$MODIPSLDIR $link_veget -lioipsl -L$NCDFLIB -lnetcdf -lioipsl -lnetcdf "
+     else
+       set opt_link="-L$MODIPSLDIR $link_veget -lioipsl -L$NCDFLIB -lnetcdf -lnetcdf "
+     endif
    else
      set opt_link=" "
@@ -804 +825 @@
    set f77=f90
    set f90=f90
-   set opt_link="-lf77compat -L$MODIPSLDIR $link_veget -lioipsl -L$NCDFLIB -lnetcdf "
+   if ( $io == "ioipsl" ) then
+     set opt_link="-lf77compat -L$MODIPSLDIR $link_veget -lioipsl -L$NCDFLIB -lnetcdf "
+   else
+     set opt_link="-lf77compat -L$MODIPSLDIR $link_veget -L$NCDFLIB -lnetcdf "
+   endif
 #################
 else if $NEC then

LMDZ4/branches/LMDZ4-dev/makelmdz_fcm

@@ -3 +3 @@
 # $Header$
 #
-# FH : on ne cr�e plus le fichier arch.mk qui est suppos� exister par
+# FH : on ne crée plus le fichier arch.mk qui est supposé exister par
 # FH : ailleurs.
-# FH : ult�rieurement, ce fichier sera pr�-existant pour une s�rie
-# FH : de configurations en versions optimis�es et debug qui seront
-# FH : li�s (ln -s) avec arch.mk en fonction de l'architecture.
-# FH : Pour le moment, cette version est en test et on peut cr�er les
-# FH : arch.mk en lan�ant une premi�re fois makegcm7.
+# FH : ulterieurement, ce fichier sera pré-existant pour une série
+# FH : de configurations en versions optimisées et debug qui seront
+# FH : liés (ln -s) avec arch.mk en fonction de l'architecture.
+# FH : Pour le moment, cette version est en test et on peut créer les
+# FH : arch.mk en lançant une première fois makegcm.
 #
 set verbose echo
@@ -58 +58 @@
 Usage :
 makegcm [options] -m arch exec
-[-h]                       : manuel abrégé
+[-h]                       : manuel abrégé
 [-d [[IMx]JMx]LM]          : IM, JM, LM sont les dims en x, y, z (def: $dim)
 [-p PHYS]                  : compilation avec la physique libf/phyPHYS, (def: lmd)
 [-prod / -dev / -debug]    : compilation en mode production (default) / developpement / debug .
-[-c false/MPI1/MPI2]       : couplé océan : MPI1/MPI2/false (def: false)
-[-v false/true]            : avec ou sans végétation (def: false)
+[-c false/MPI1/MPI2]       : couplé océan : MPI1/MPI2/false (def: false)
+[-v false/true]            : avec ou sans végétation (def: false)
 [-chimie INCA/false]       : avec ou sans model de chimie INCA (def: false)
 [-parallel none/mpi/omp/mpi_omp] : parallelisation (default: none) : mpi, openmp ou mixte mpi_openmp
@@ -75 +75 @@
 [-fcm_path path]           : chemin pour fcm (def: le chemin est suppose deja exister dans le PATH)
  -arch nom_arch            : nom de l'architecture cible
- exec                      : exécutable généré
+ exec                      : exécutable généré
 fin
         exit
@@ -157 +157 @@
 
 ###############################################################
-# lecture des chemins propres a l'architecture de la machine #
+# lecture des chemins propres à l'architecture de la machine #
 ###############################################################
 
@@ -180 +180 @@
    
 else
-   set CPP_KEY="$CPP_KEY CPP_PHYS"
+   #Default planet type is Earth
+   set CPP_KEY="$CPP_KEY CPP_EARTH"
 endif