1 | ## $Header$ |
---|
2 | # |
---|
3 | ## Planete: |
---|
4 | planet_type=titan |
---|
5 | # |
---|
6 | ## nombre de pas par jour (multiple de iperiod) ( ici pour dt ~ 2 min ) |
---|
7 | day_step=8000 |
---|
8 | ## periode pour le pas Matsuno (en pas) |
---|
9 | iperiod=5 |
---|
10 | ## periode de la dissipation (en pas) |
---|
11 | ## dissip_period=5 PAR DEFAUT: 0 ie calcul autom. (=> 55 pour runs recents) |
---|
12 | ## choix de l'operateur de dissipation (star ou non star ) |
---|
13 | lstardis=y |
---|
14 | ## nombre d'iterations de l'operateur de dissipation gradiv |
---|
15 | nitergdiv=1 |
---|
16 | ## nombre d'iterations de l'operateur de dissipation nxgradrot |
---|
17 | nitergrot=2 |
---|
18 | ## nombre d'iterations de l'operateur de dissipation divgrad |
---|
19 | niterh=2 |
---|
20 | ## temps de dissipation des plus petites long.d ondes pour u,v (gradiv) |
---|
21 | tetagdiv=2.e5 |
---|
22 | ## temps de dissipation des plus petites long.d ondes pour u,v(nxgradrot) |
---|
23 | tetagrot=2.e5 |
---|
24 | ## temps de dissipation des plus petites long.d ondes pour h ( divgrad) |
---|
25 | tetatemp=2.e5 |
---|
26 | ## coefficient pour gamdissip |
---|
27 | coefdis=0. |
---|
28 | ## choix du shema d'integration temporelle (Matsuno ou Matsuno-leapfrog) |
---|
29 | purmats=n |
---|
30 | ## avec ou sans physique |
---|
31 | ## 0: pas de physique (e.g. en mode Shallow Water) |
---|
32 | ## 1: avec physique (e.g. physique phylmd) |
---|
33 | ## 2: avec rappel newtonien dans la dynamique |
---|
34 | iflag_phys=1 |
---|
35 | ## avec ou sans fichiers de demarrage (start.nc, startphy.nc) ? |
---|
36 | ## (sans fichiers de demarrage, initialisation des champs par iniacademic |
---|
37 | ## dans la dynamique) PAS AU POINT POUR TITAN |
---|
38 | read_start=y |
---|
39 | ## periode de la physique (en pas) |
---|
40 | iphysiq=40 |
---|
41 | ## avec ou sans traceurs |
---|
42 | iflag_trac=1 |
---|
43 | ## Avec ou sans strato // i.e. Couche eponge et second palier pour dissip horiz. |
---|
44 | ok_strato=y |
---|
45 | ## Dissipation horizontale |
---|
46 | dissip_fac_mid=2. |
---|
47 | dissip_fac_up=10. |
---|
48 | # deltaz et hdelta en km |
---|
49 | dissip_deltaz=10. |
---|
50 | dissip_hdelta=5. |
---|
51 | # pupstart en Pa |
---|
52 | dissip_pupstart=1.e3 |
---|
53 | ## Couche eponge |
---|
54 | # 1: dans les 4 derniers niveaux |
---|
55 | # 2: dans les couches de pression plus faible que 100 fois la pression de la derniere couche |
---|
56 | iflag_top_bound=1 |
---|
57 | ## Mode Couche eponge |
---|
58 | # mode = 0 : pas de sponge |
---|
59 | # mode = 1 : u et v -> 0 |
---|
60 | # mode = 2 : u et v -> moyenne zonale |
---|
61 | # mode = 3 : u, v et h -> moyenne zonale |
---|
62 | mode_top_bound=1 |
---|
63 | # Coefficient pour la couche eponge (valeur derniere couche) |
---|
64 | tau_top_bound=4.e-5 |
---|
65 | ## Maree gravitationnelle ou non |
---|
66 | tidal=y |
---|
67 | |
---|
68 | ## longitude en degres du centre du zoom |
---|
69 | clon=0. |
---|
70 | ## latitude en degres du centre du zoom |
---|
71 | clat=0. |
---|
72 | ## facteur de grossissement du zoom,selon longitude |
---|
73 | grossismx=1.0 |
---|
74 | ## facteur de grossissement du zoom ,selon latitude |
---|
75 | grossismy=1.0 |
---|
76 | ## Fonction f(y) hyperbolique si = .true. , sinon sinusoidale |
---|
77 | fxyhypb=y |
---|
78 | ## extension en longitude de la zone du zoom ( fraction de la zone totale) |
---|
79 | dzoomx=0.0 |
---|
80 | ## extension en latitude de la zone du zoom ( fraction de la zone totale) |
---|
81 | dzoomy=0.0 |
---|
82 | ##raideur du zoom en X |
---|
83 | taux=3. |
---|
84 | ##raideur du zoom en Y |
---|
85 | tauy=3. |
---|
86 | ## Fonction f(y) avec y = Sin(latit.) si = .true. , sinon y = latit. |
---|
87 | ysinus=y |
---|