| [5] | 1 | ## $Header$ |
|---|
| [6] | 2 | # |
|---|
| 3 | ## Planete: |
|---|
| 4 | planet_type=titan |
|---|
| 5 | # |
|---|
| [5] | 6 | ## nombre de pas par jour (multiple de iperiod) ( ici pour dt ~ 2 min ) |
|---|
| [122] | 7 | day_step=8000 |
|---|
| [5] | 8 | ## periode pour le pas Matsuno (en pas) |
|---|
| 9 | iperiod=5 |
|---|
| [1056] | 10 | ## periode de la dissipation (en pas) PAR DEFAUT: 0 ie calcul autom. |
|---|
| 11 | ## dissip_period=5 |
|---|
| [5] | 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 |
|---|
| [6] | 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 |
|---|
| [5] | 34 | iflag_phys=1 |
|---|
| [6] | 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 |
|---|
| [5] | 39 | ## periode de la physique (en pas) |
|---|
| [122] | 40 | iphysiq=40 |
|---|
| [5] | 41 | ## avec ou sans traceurs |
|---|
| 42 | iflag_trac=1 |
|---|
| [108] | 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 |
|---|
| [1056] | 52 | dissip_pupstart=1.e2 |
|---|
| [108] | 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 |
|---|
| [495] | 65 | ## Maree gravitationnelle ou non |
|---|
| 66 | tidal=y |
|---|
| 67 | |
|---|
| [5] | 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 |
|---|
