## $Header$
#
## Planete:
planet_type=titan 
# 
## nombre de pas par jour (multiple de iperiod) ( ici pour  dt ~ 2 min )      
day_step=12000
## periode pour le pas Matsuno (en pas)
iperiod=5
## periode de la dissipation (en pas)
idissip=5
## choix de l'operateur de dissipation (star ou  non star )
lstardis=y
## 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=2.e5
## temps de dissipation des plus petites long.d ondes pour u,v(nxgradrot)
tetagrot=2.e5
## temps de dissipation des plus petites long.d ondes pour  h ( divgrad) 
tetatemp=2.e5
## coefficient pour gamdissip                                            
coefdis=0.
## choix du shema d'integration temporelle (Matsuno ou Matsuno-leapfrog) 
purmats=n
## avec ou sans physique
## 0: pas de physique (e.g. en mode Shallow Water)
## 1: avec physique (e.g. physique phylmd)
## 2: avec rappel newtonien dans la dynamique                                         
iflag_phys=1
## avec ou sans fichiers de demarrage (start.nc, startphy.nc) ?
## (sans fichiers de demarrage, initialisation des champs par iniacademic
##  dans la dynamique) PAS AU POINT POUR TITAN
read_start=y
## periode de la physique (en pas)                                       
iphysiq=20
## avec ou sans traceurs                                                 
iflag_trac=1
##  Avec ou sans strato // i.e. Couche eponge et second palier pour dissip horiz.
ok_strato=y
## Dissipation horizontale
dissip_fac_mid=2.
dissip_fac_up=10.
# deltaz et hdelta en km
dissip_deltaz=10.
dissip_hdelta=5.
# pupstart en Pa
dissip_pupstart=1.e3
## Couche eponge 
#   1: dans les 4 derniers niveaux
#   2: dans les couches de pression plus faible que 100 fois la pression de la derniere couche
iflag_top_bound=1
## Mode Couche eponge 
#   mode = 0 : pas de sponge
#   mode = 1 : u et v -> 0
#   mode = 2 : u et v -> moyenne zonale
#   mode = 3 : u, v et h -> moyenne zonale
mode_top_bound=1
#  Coefficient pour la couche eponge (valeur derniere couche)
tau_top_bound=4.e-5
#
## 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.0
## facteur de grossissement du zoom ,selon latitude                      
grossismy=1.0
##  Fonction  f(y)  hyperbolique  si = .true.  , sinon  sinusoidale         
fxyhypb=y
## extension en longitude  de la zone du zoom  ( fraction de la zone totale)
dzoomx=0.0
## extension en latitude de la zone  du zoom  ( fraction de la zone totale)
dzoomy=0.0
##raideur du zoom en  X
taux=3.
##raideur du zoom en  Y
tauy=3.
##  Fonction  f(y) avec y = Sin(latit.) si = .true. , sinon y = latit.         
ysinus=y