| 1 | ## $Id: gcm.def 1403 2010-07-01 09:02:53Z fairhead $ |
|---|
| 2 | ## nombre de pas par jour (multiple de iperiod) ( ici pour dt = 1 min ) |
|---|
| 3 | day_step=480 |
|---|
| 4 | ## periode pour le pas Matsuno (en pas) |
|---|
| 5 | iperiod=5 |
|---|
| 6 | ## periode de la dissipation (en pas) |
|---|
| 7 | idissip=5 |
|---|
| 8 | ## choix de l'operateur de dissipation (star ou non star ) |
|---|
| 9 | lstardis=y |
|---|
| 10 | ## nombre d'iterations de l'operateur de dissipation gradiv |
|---|
| 11 | nitergdiv=1 |
|---|
| 12 | ## nombre d'iterations de l'operateur de dissipation nxgradrot |
|---|
| 13 | nitergrot=2 |
|---|
| 14 | ## nombre d'iterations de l'operateur de dissipation divgrad |
|---|
| 15 | niterh=2 |
|---|
| 16 | ## temps de dissipation des plus petites long.d ondes pour u,v (gradiv) |
|---|
| 17 | tetagdiv=5400. |
|---|
| 18 | ## temps de dissipation des plus petites long.d ondes pour u,v(nxgradrot) |
|---|
| 19 | #tetagrot=5400. |
|---|
| 20 | tetagrot=5400. |
|---|
| 21 | ## temps de dissipation des plus petites long.d ondes pour h ( divgrad) |
|---|
| 22 | tetatemp=5400. |
|---|
| 23 | ## coefficient pour gamdissip |
|---|
| 24 | coefdis=0. |
|---|
| 25 | ## choix du shema d'integration temporelle (Matsuno:y ou Matsuno-leapfrog:n) |
|---|
| 26 | purmats=n |
|---|
| 27 | ## avec ou sans physique |
|---|
| 28 | ## 0: pas de physique (e.g. en mode Shallow Water) |
|---|
| 29 | ## 1: avec physique (e.g. physique phylmd) |
|---|
| 30 | ## 2: avec rappel newtonien dans la dynamique |
|---|
| 31 | #iflag_phys=1 |
|---|
| 32 | |
|---|
| 33 | ## periode de la physique (en pas dynamiques, n'a de sens que si iflag_phys=1) |
|---|
| 34 | iphysiq=5 |
|---|
| 35 | # |
|---|
| 36 | #### Verticale discretisation |
|---|
| 37 | ## Modified verticale discretisation for stratospheric versions. |
|---|
| 38 | ## The new version can be controled by following parameters |
|---|
| 39 | ## Al the parameters are given in km assuming a given scalehigh |
|---|
| 40 | #vert_sampling=strato_custom |
|---|
| 41 | ## scale hight |
|---|
| 42 | vert_scale=7. |
|---|
| 43 | ## width of first layer |
|---|
| 44 | vert_dzmin=0.017 |
|---|
| 45 | ## dz in the low atmosphere |
|---|
| 46 | vert_dzlow=1. |
|---|
| 47 | ## height at which resolution recches dzlow |
|---|
| 48 | vert_z0low=8.7 |
|---|
| 49 | ## dz in the mid atmsophere |
|---|
| 50 | vert_dzmid=2. |
|---|
| 51 | ## height at which resolution recches dzmid |
|---|
| 52 | vert_z0mid=70. |
|---|
| 53 | ## width of the transition |
|---|
| 54 | vert_h_mid=20. |
|---|
| 55 | ## dz in the high atmsophere |
|---|
| 56 | vert_dzhig=11. |
|---|
| 57 | ## height at which resolution recches dz |
|---|
| 58 | vert_z0hig=75. |
|---|
| 59 | ## width of the transition |
|---|
| 60 | vert_h_hig=20. |
|---|
| 61 | ## vert_prof_dissip used to choose the vertical profile |
|---|
| 62 | ## of horizontal dissipation (values are 0 and 1) |
|---|
| 63 | vert_prof_dissip = 1 |
|---|
| 64 | dissip_zref = 50 |
|---|
| 65 | ## Avec ou sans strato |
|---|
| 66 | ok_strato=y |
|---|
| 67 | ok_hines=y |
|---|
| 68 | # Couche eponge dans les couches de pression plus faible que 100 fois la pression de la derniere couche |
|---|
| 69 | iflag_top_bound=2 |
|---|
| 70 | # Coefficient pour la couche eponge (valeur derniere couche) |
|---|
| 71 | tau_top_bound=5.e-4 |
|---|
| 72 | # |
|---|
| 73 | ### Parametres du zoom |
|---|
| 74 | ## longitude en degres du centre du zoom |
|---|
| 75 | clon=-90. |
|---|
| 76 | ## latitude en degres du centre du zoom |
|---|
| 77 | clat=-15. |
|---|
| 78 | ## facteur de grossissement du zoom,selon longitude |
|---|
| 79 | grossismx=1.0 |
|---|
| 80 | ## facteur de grossissement du zoom ,selon latitude |
|---|
| 81 | grossismy=1.0 |
|---|
| 82 | ## Fonction f(y) hyperbolique si = .true. , sinon sinusoidale |
|---|
| 83 | fxyhypb=y |
|---|
| 84 | ## extension en longitude de la zone du zoom ( fraction de la zone totale) |
|---|
| 85 | dzoomx=0.15 |
|---|
| 86 | ## extension en latitude de la zone du zoom ( fraction de la zone totale) |
|---|
| 87 | dzoomy=0.15 |
|---|
| 88 | ##raideur du zoom en X |
|---|
| 89 | taux=3. |
|---|
| 90 | ##raideur du zoom en Y |
|---|
| 91 | tauy=3. |
|---|
| 92 | ## Fonction f(y) avec y = Sin(latit.) si = .true. , sinon y = latit. |
|---|
| 93 | ysinus=y |
|---|
