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

# Planete
planet_type = titan

# Nombre de jours d'integration                                         
nday = 11

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

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

# Periode de sortie des variables de controle (en pas)
iconser = 1000

# Periode d'ecriture du fichier histoire (en jour)
iecri = 100

# Periode de stockage fichier histmoy (en jour)
periodav = 60.

# Periode de la dissipation (en pas)
dissip_period = 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 = 25000.

# Temps de dissipation des plus petites long. d'ondes pour u,v (nxgradrot)
tetagrot = 25000.

# Temps de dissipation des plus petites long. d'ondes pour h (divgrad)
tetatemp = 25000.

# Facteur d'amplification de la dissip en haut du modele
dissip_fac_up = 30.

# Intervalle (km) pour le changement mid / up de la dissip
# dissip_deltaz = 30.

# Pression (Pa) au bas la transition mid / up
dissip_pupstart = 1.e1

# Coefficient pour gamdissip
coefdis = 0.

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

# Avec ou sans physique
physic = .true.

# Periode de la physique (en pas)
iphysiq = 40

# Titan : marees gravitationnelles
tidal = .true.

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

# Frequence (en pas) de l'ecriture du fichier diagfi
ecritphy = 8000

# Longitude en degres du centre du zoom
clon = 0.

# Latitude en degres du centre du zoom
clat = 0.

# Facteur de grossissement du zoom, selon longitude
grossismx = 1.

# Facteur de grossissement du zoom, selon latitude
grossismy = 1.

# Fonction f(y) hyperbolique si = .true., sinon sinusoidale
fxyhypb = .false.

# Extension en longitude de la zone du zoom (fraction de la zone totale)
dzoomx = 0.

# Extension en latitude de la zone du zoom (fraction de la zone totale)
dzoomy = 0.

# Raideur du zoom en X
taux = 2.

# Raideur du zoom en Y
tauy = 2.

# Fonction f(y) avec y = Sin(latit.), si = .true., sinon y = latit.
ysinus = .false.

# Avec sponge layer
callsponge = .true.

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

# Sponge : nsponge
nsponge = 4

# Sponge : tetasponge (secondes)
tetasponge = 10000

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