#
#-----------------------------------------------------------------------
#Parametres de controle du run:
#------------------------------

 planet_type = generic

# Nombre de jours d'integration
 nday = 57

# nombre de pas par jour (multiple de iperiod) ( ici pour dt = 1 min )
 day_step = 4800

 day_ini = 0

# periode de la physique (en pas)
  iphysiq=10

# periode pour le pas Matsuno (en pas)
 iperiod=5

# periode de stockage fichier histmoy (en jour)
 periodav=70.

# periode de la dissipation (en pas)
  idissip=5

# choix de l'operateur de dissipation (star ou  non star )
 lstardis=.true.

# avec ou sans coordonnee hybrides
 hybrid=.true.

# nombre d'iterations de l'operateur de dissipation   gradiv
nitergdiv=1

# nombre d'iterations de l'operateur de dissipation  nxgradrot
nitergrot=2

# nombre d'iterations de l'operateur de dissipation  divgrad
niterh=2

# temps de dissipation des plus petites long.d ondes pour u,v (gradiv)
tetagdiv=80000. # weak 50000 strong 5000

# temps de dissipation des plus petites long.d ondes pour u,v(nxgradrot)
tetagrot=200000. # weak 20000 strong 2000

# temps de dissipation des plus petites long.d ondes pour  h ( divgrad)
tetatemp=200000. # weak 20000 strong 2000

# Values of dissip_factz are all defaults. Used ?
dissip_factz =    2.00000000000000

# Values of dissip_deltaz are all defaults. Used ?
dissip_deltaz =    10.

# Values of dissip_zref are all defaults.
dissip_zref =    100.

# Values of dissip_fac_mid are all defaults.
dissip_fac_mid =    1

# Values of dissip_fac_up are all defaults.
dissip_fac_up =    8.

# Values of dissip_hdelta are all defaults. Higher values = faster transition between mid and up
dissip_hdelta =    10.

# Values of dissip_pupstart are all defaults. Where is the transition (in Pa)
dissip_pupstart =   1000.

# coefficient pour gamdissip
  coefdis=0.

# choix du shema d'integration temporelle (Matsuno ou Matsuno-leapfrog)
  purmats=.false.

# avec ou sans physique
   physic=.true.

# choix d'une grille reguliere
  grireg=.true.

# Avec sponge layer
  callsponge  = .true.

# Sponge:  mode0(u=v=0), mode1(u=umoy,v=0), mode2(u=umoy,v=vmoy)
  mode_sponge= 2

# Sponge:  number of sponge layers
  nsponge= 3

# Sponge:  tetasponge (secondes)
  tetasponge = 50000

# some definitions for the physics, in file 'callphys.def'
INCLUDEDEF=callphys.def

