! $Id: conf_gcm.F90 4519 2023-04-24 16:11:11Z evignon $

  USE control_mod

#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 infotrac, ONLY : type_trac

  use assert_m, only: assert

  USE comconst_mod, ONLY: dissip_deltaz, dissip_factz, dissip_zref, &

                          iflag_top_bound, mode_top_bound, tau_top_bound, &

                          ngroup, maxlatfilter

  USE logic_mod, ONLY: fxyhypb, iflag_phys, ok_etat0, ok_gradsfile, &

                       ok_guide, ok_limit, ok_strato, purmats, read_start, &

                       ysinus, read_orop

  USE serre_mod, ONLY: clon,clat,grossismx,grossismy,dzoomx,dzoomy, &

                       alphax,alphay,taux,tauy

  USE temps_mod, ONLY: calend, year_len

  IMPLICIT NONE

  !-----------------------------------------------------------------------

  !     Auteurs :   L. Fairhead , P. Le Van  .

  !     Arguments :

  !     tapedef   :

  !     etatinit  :     = TRUE   , on ne  compare pas les valeurs des para-

  !     -metres  du zoom  avec  celles lues sur le fichier start .

  LOGICAL,INTENT(IN) :: etatinit

  INTEGER,INTENT(IN) :: tapedef

  !   Declarations :

  !   --------------

  include "dimensions.h"

  include "paramet.h"

  include "comdissnew.h"

  include "iniprint.h"

  !   local:

  !   ------

  REAL clonn,clatt,grossismxx,grossismyy

  REAL dzoomxx,dzoomyy, tauxx,tauyy

  LOGICAL  fxyhypbb, ysinuss

  !  -------------------------------------------------------------------

  !       .........     Version  du 29/04/97       ..........

  !   Nouveaux parametres nitergdiv,nitergrot,niterh,tetagdiv,tetagrot,

  !      tetatemp   ajoutes  pour la dissipation   .

  !   Autre parametre ajoute en fin de liste de tapedef : ** fxyhypb **

  !  Si fxyhypb = .TRUE. , choix de la fonction a derivee tangente hyperb.

  !    Sinon , choix de fxynew  , a derivee sinusoidale  ..

  !   ......  etatinit = . TRUE. si defrun  est appele dans ETAT0_LMD  ou

  !         LIMIT_LMD  pour l'initialisation de start.dat (dic) et

  !                de limit.dat ( dic)                        ...........

  !           Sinon  etatinit = . FALSE .

  !   Donc etatinit = .F.  si on veut comparer les valeurs de  grossismx ,

  !    grossismy,clon,clat, fxyhypb  lues sur  le fichier  start  avec

  !   celles passees  par run.def ,  au debut du gcm, apres l'appel a

  !    lectba .

  !   Ces parmetres definissant entre autres la grille et doivent etre

  !   pareils et coherents , sinon il y aura  divergence du gcm .

  !-----------------------------------------------------------------------

  !   initialisations:

  !   ----------------

  !Config  Key  = lunout

  !Config  Desc = unite de fichier pour les impressions

  !Config  Def  = 6

  !Config  Help = unite de fichier pour les impressions

  !Config         (defaut sortie standard = 6)

     OPEN(UNIT=lunout,FILE='lmdz.out',ACTION='write',                      &

  ENDIF

  !Config  Key  = prt_level

  !Config  Desc = niveau d'impressions de d\'ebogage

  !Config  Def  = 0

  !Config  Help = Niveau d'impression pour le d\'ebogage

  !Config         (0 = minimum d'impression)

  !-----------------------------------------------------------------------

  !  Parametres de controle du run:

  !-----------------------------------------------------------------------

  !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

  !Config  Key  = calend

  !Config  Desc = type de calendrier utilise

  !Config  Def  = earth_360d

  !Config  Help = valeur possible: earth_360d, earth_365d, earth_366d

  !Config

! initialize year_len for aquaplanets and 1D

  ELSE

  ENDIF

  !Config  Key  = dayref

  !Config  Desc = Jour de l'etat initial

  !Config  Def  = 1

  !Config  Help = Jour de l'etat initial ( = 350  si 20 Decembre ,

  !Config         par expl. ,comme ici ) ... A completer

  !Config  Key  = anneeref

  !Config  Desc = Annee de l'etat initial

  !Config  Def  = 1998

  !Config  Help = Annee de l'etat  initial

  !Config         (   avec  4  chiffres   ) ... A completer

  !Config  Key  = raz_date

  !Config  Desc = Remise a zero de la date initiale

  !Config  Def  = 0 (pas de remise a zero)

  !Config  Help = Remise a zero de la date initiale

  !Config         0 pas de remise a zero, on garde la date du fichier restart

  !Config         1 prise en compte de la date de gcm.def avec remise a zero

  !Config         des compteurs de pas de temps

  !Config  Key  = resetvarc

  !Config  Desc = Reinit des variables de controle

  !Config  Def  = n

  !Config  Help = Reinit des variables de controle

  !Config  Key  = nday

  !Config  Desc = Nombre de jours d'integration

  !Config  Def  = 10

  !Config  Help = Nombre de jours d'integration

  !Config         ... On pourait aussi permettre des mois ou des annees !

  !Config  Key  = starttime

  !Config  Desc = Heure de depart de la simulation

  !Config  Def  = 0

  !Config  Help = Heure de depart de la simulation

  !Config         en jour

  !Config  Key  = day_step

  !Config  Desc = nombre de pas par jour

  !Config  Def  = 240

  !Config  Help = nombre de pas par jour (multiple de iperiod) (

  !Config          ici pour  dt = 1 min )

  !Config  Key  = nsplit_phys

  !Config  Key  = iperiod

  !Config  Desc = periode pour le pas Matsuno

  !Config  Def  = 5

  !Config  Help = periode pour le pas Matsuno (en pas de temps)

  !Config  Key  = iapp_tracvl

  !Config  Desc = frequence du groupement des flux

  !Config  Def  = iperiod

  !Config  Help = frequence du groupement des flux (en pas de temps)

  !Config  Key  = iconser

  !Config  Desc = periode de sortie des variables de controle

  !Config  Def  = 240

  !Config  Help = periode de sortie des variables de controle

  !Config         (En pas de temps)

  !Config  Key  = iecri

  !Config  Desc = periode d'ecriture du fichier histoire

  !Config  Def  = 1

  !Config  Help = periode d'ecriture du fichier histoire (en jour)

  !Config  Key  = periodav

  !Config  Desc = periode de stockage fichier histmoy

  !Config  Def  = 1

  !Config  Help = periode de stockage fichier histmoy (en jour)

  !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

  !Config  Key  = dissip_period

  !Config  Desc = periode de la dissipation

  !Config  Def  = 0

  !Config  Help = periode de la dissipation

  !Config  dissip_period=0 => la valeur sera calcule dans inidissip

  !Config  dissip_period>0 => on prend cette valeur

  !cc  ....   P. Le Van , modif le 29/04/97 .pour la dissipation  ...

  !cc

  !Config  Key  = lstardis

  !Config  Desc = choix de l'operateur de dissipation

  !Config  Def  = y

  !Config  Help = choix de l'operateur de dissipation

  !Config         'y' si on veut star et 'n' si on veut non-start !

  !Config         Moi y en a pas comprendre !

  !Config  Key  = nitergdiv

  !Config  Desc = Nombre d'iteration de gradiv

  !Config  Def  = 1

  !Config  Help = nombre d'iterations de l'operateur de dissipation

  !Config         gradiv

  !Config  Key  = nitergrot

  !Config  Desc = nombre d'iterations de nxgradrot

  !Config  Def  = 2

  !Config  Help = nombre d'iterations de l'operateur de dissipation

  !Config         nxgradrot

  !Config  Key  = niterh

  !Config  Desc = nombre d'iterations de divgrad

  !Config  Def  = 2

  !Config  Help = nombre d'iterations de l'operateur de dissipation

  !Config         divgrad

  !Config  Key  = tetagdiv

  !Config  Desc = temps de dissipation pour div

  !Config  Def  = 7200

  !Config  Help = temps de dissipation des plus petites longeur

  !Config         d'ondes pour u,v (gradiv)

  !Config  Key  = tetagrot

  !Config  Desc = temps de dissipation pour grad

  !Config  Def  = 7200

  !Config  Help = temps de dissipation des plus petites longeur

  !Config         d'ondes pour u,v (nxgradrot)

  !Config  Key  = tetatemp

  !Config  Desc = temps de dissipation pour h

  !Config  Def  = 7200

  !Config  Help =  temps de dissipation des plus petites longeur

  !Config         d'ondes pour h (divgrad)

  ! Parametres controlant la variation sur la verticale des constantes de

  ! dissipation.

  ! Pour le moment actifs uniquement dans la version a 39 niveaux

  ! avec ok_strato=y

  ! maxlatfilter

  ! ngroup

  ! top_bound sponge: only active if ok_strato=.true. and iflag_top_bound!=0

  !                   iflag_top_bound=0 for no sponge

  !                   iflag_top_bound=1 for sponge over 4 topmost layers

  !                   iflag_top_bound=2 for sponge from top to ~1% of top layer pressure

  ! mode_top_bound : fields towards which sponge relaxation will be done:

  !                  mode_top_bound=0: no relaxation

  !                  mode_top_bound=1: u and v relax towards 0

  !                  mode_top_bound=2: u and v relax towards their zonal mean

  !                  mode_top_bound=3: u,v and pot. temp. relax towards their zonal mean

  ! top_bound sponge : inverse of charactericstic relaxation time scale for sponge

  !Config  Key  = coefdis

  !Config  Desc = coefficient pour gamdissip

  !Config  Def  = 0

  !Config  Help = coefficient pour gamdissip

  !Config  Key  = purmats

  !Config  Desc = Schema d'integration

  !Config  Def  = n

  !Config  Help = Choix du schema d'integration temporel.

  !Config         y = pure Matsuno sinon c'est du Matsuno-leapfrog

  !Config  Key  = ok_guide

  !Config  Desc = Guidage

  !Config  Def  = n

  !Config  Help = Guidage

  !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

  !Config  Key  = iflag_phys

  !Config  Desc = Avec ls physique

  !Config  Def  = 1

  !Config  Help = Permet de faire tourner le modele sans

  !Config         physique.

  !Config  Key  =  iphysiq

  !Config  Desc = Periode de la physique

  !Config  Def  = 5

  !Config  Help = Periode de la physique en pas de temps de la dynamique.

  !Config  Key  = ip_ebil_dyn

  !Config  Desc = PRINT level for energy conserv. diag.

  !Config  Def  = 0

  !Config  Help = PRINT level for energy conservation diag. ;

  !               les options suivantes existent :

  !Config         0 pas de print

  !Config         1 pas de print

  !Config         2 print,

  !cc  ....   P. Le Van , ajout  le 7/03/95 .pour le zoom ...

  !     .........   (  modif  le 17/04/96 )   .........

     !Config  Key  = clon

     !Config  Desc = centre du zoom, longitude

     !Config  Def  = 0

     !Config  Help = longitude en degres du centre

     !Config         du zoom

     !Config  Key  = clat

     !Config  Desc = centre du zoom, latitude

     !Config  Def  = 0

     !Config  Help = latitude en degres du centre du zoom

     !Config

     ENDIF

     !Config  Key  = grossismx

     !Config  Desc = zoom en longitude

     !Config  Def  = 1.0

     !Config  Help = facteur de grossissement du zoom,

     !Config         selon la longitude

     ENDIF

     !Config  Key  = grossismy

     !Config  Desc = zoom en latitude

     !Config  Def  = 1.0

     !Config  Help = facteur de grossissement du zoom,

     !Config         selon la latitude

     ENDIF

        write(lunout,*) &

     ELSE

     ENDIF

        write(lunout,*) &

     ELSE

     ENDIF

     !    alphax et alphay sont les anciennes formulat. des grossissements

     !Config  Key  = fxyhypb

     !Config  Desc = Fonction  hyperbolique

     !Config  Def  = y

     !Config  Help = Fonction  f(y)  hyperbolique  si = .true.

     !Config         sinon  sinusoidale

        ENDIF

     ELSE

        ENDIF

     ENDIF

     !Config  Key  = dzoomx

     !Config  Desc = extension en longitude

     !Config  Def  = 0

     !Config  Help = extension en longitude  de la zone du zoom

     !Config         ( fraction de la zone totale)

        ENDIF

     ENDIF

     !Config  Key  = dzoomy

     !Config  Desc = extension en latitude

     !Config  Def  = 0

     !Config  Help = extension en latitude de la zone  du zoom

     !Config         ( fraction de la zone totale)

        ENDIF

     ENDIF

     !Config  Key  = taux

     !Config  Desc = raideur du zoom en  X

     !Config  Def  = 3

     !Config  Help = raideur du zoom en  X

        ENDIF

     ENDIF

     !Config  Key  = tauyy

     !Config  Desc = raideur du zoom en  Y

     !Config  Def  = 3

     !Config  Help = raideur du zoom en  Y

        ENDIF

     ENDIF

     !c

        !Config  Key  = ysinus

        !Config  IF   = !fxyhypb

        !Config  Desc = Fonction en Sinus

        !Config  Def  = y

        !Config  Help = Fonction  f(y) avec y = Sin(latit.) si = .true.

        !Config         sinon y = latit.

           ENDIF

        ELSE

           ENDIF

        ENDIF

     ENDIF ! of IF( .NOT.fxyhypb  )

     !Config  Key  = offline

     !Config  Desc = Nouvelle eau liquide

     !Config  Def  = n

     !Config  Help = Permet de mettre en route la

     !Config         nouvelle parametrisation de l'eau liquide !

     !Config  Key  = type_trac

     !Config  Desc = Choix de couplage avec model de chimie INCA ou REPROBUS

     !Config  Def  = lmdz

     !Config  Help =

     !Config         'lmdz' = pas de couplage, pur LMDZ

     !Config         'inca' = model de chime INCA

     !Config         'repr' = model de chime REPROBUS

     !Config         'inco' = INCA + CO2i (temporaire)

     !Config  Key  = ok_dynzon

     !Config  Desc = calcul et sortie des transports

     !Config  Def  = n

     !Config  Help = Permet de mettre en route le calcul des transports

     !Config

     !Config  Key  = ok_dyn_ins

     !Config  Desc = sorties instantanees dans la dynamique

     !Config  Def  = n

     !Config  Help =

     !Config

     !Config  Key  = ok_dyn_ave

     !Config  Desc = sorties moyennes dans la dynamique

     !Config  Def  = n

     !Config  Help =

     !Config

  ELSE

     !Config  Key  = clon

     !Config  Desc = centre du zoom, longitude

     !Config  Def  = 0

     !Config  Help = longitude en degres du centre

     !Config         du zoom

     !Config  Key  = clat

     !Config  Desc = centre du zoom, latitude

     !Config  Def  = 0

     !Config  Help = latitude en degres du centre du zoom

     !Config

     !Config  Key  = grossismx

     !Config  Desc = zoom en longitude

     !Config  Def  = 1.0

     !Config  Help = facteur de grossissement du zoom,

     !Config         selon la longitude

     !Config  Key  = grossismy

     !Config  Desc = zoom en latitude

     !Config  Def  = 1.0

     !Config  Help = facteur de grossissement du zoom,

     !Config         selon la latitude

        write(lunout,*) &

     ELSE

     ENDIF

     ELSE

     ENDIF

     !    alphax et alphay sont les anciennes formulat. des grossissements

     !Config  Key  = fxyhypb

     !Config  Desc = Fonction  hyperbolique

     !Config  Def  = y

     !Config  Help = Fonction  f(y)  hyperbolique  si = .true.

     !Config         sinon  sinusoidale

     !Config  Key  = dzoomx

     !Config  Desc = extension en longitude

     !Config  Def  = 0

     !Config  Help = extension en longitude  de la zone du zoom

     !Config         ( fraction de la zone totale)

     !Config  Key  = dzoomy

     !Config  Desc = extension en latitude

     !Config  Def  = 0

     !Config  Help = extension en latitude de la zone  du zoom

     !Config         ( fraction de la zone totale)

     !Config  Key  = taux

     !Config  Desc = raideur du zoom en  X

     !Config  Def  = 3

     !Config  Help = raideur du zoom en  X

     !Config  Key  = tauy

     !Config  Desc = raideur du zoom en  Y

     !Config  Def  = 3

     !Config  Help = raideur du zoom en  Y

     !Config  Key  = ysinus

     !Config  IF   = !fxyhypb

     !Config  Desc = Fonction en Sinus

     !Config  Def  = y

     !Config  Help = Fonction  f(y) avec y = Sin(latit.) si = .true.

     !Config         sinon y = latit.

     !Config  Key  = offline

     !Config  Desc = Nouvelle eau liquide

     !Config  Def  = n

     !Config  Help = Permet de mettre en route la

     !Config         nouvelle parametrisation de l'eau liquide !

     !Config  Key  = type_trac

     !Config  Desc = Choix de couplage avec model de chimie INCA ou REPROBUS

     !Config  Def  = lmdz

     !Config  Help =

     !Config         'lmdz' = pas de couplage, pur LMDZ

     !Config         'inca' = model de chime INCA

     !Config         'repr' = model de chime REPROBUS

     !Config         'inco' = INCA + CO2i (temporaire)

     !Config  Key  = ok_dynzon

     !Config  Desc = sortie des transports zonaux dans la dynamique

     !Config  Def  = n

     !Config  Help = Permet de mettre en route le calcul des transports

     !Config

     !Config  Key  = ok_dyn_ins

     !Config  Desc = sorties instantanees dans la dynamique

     !Config  Def  = n

     !Config  Help =

     !Config

     !Config  Key  = ok_dyn_ave

     !Config  Desc = sorties moyennes dans la dynamique

     !Config  Def  = n

     !Config  Help =

     !Config

     !Config key = ok_strato

     !Config  Desc = activation de la version strato

     !Config  Def  = .FALSE.

     !Config  Help = active la version stratosph\'erique de LMDZ de F. Lott

     call assert(vert_prof_dissip == 0 .or. vert_prof_dissip ==  1, &

     !Config  Key  = ok_gradsfile

     !Config  Desc = activation des sorties grads du guidage

     !Config  Def  = n

     !Config  Help = active les sorties grads du guidage

     !Config  Key  = ok_limit

     !Config  Desc = creation des fichiers limit dans create_etat0_limit

     !Config  Def  = y

     !Config  Help = production du fichier limit.nc requise

     !Config  Key  = ok_etat0

     !Config  Desc = creation des fichiers etat0 dans create_etat0_limit

     !Config  Def  = y

     !Config  Help = production des fichiers start.nc, startphy.nc requise

     !Config  Key  = read_orop

     !Config  Desc = lecture du fichier de params orographiques sous maille

     !Config  Def  = f

     !Config  Help = lecture fichier plutot que grid_noro

  ENDIF test_etatinit