| 1 | |
|---|
| 2 | ######### |
|---|
| 3 | ######### THIS IS INTENDED TO WORK BOTH WITH LMDZ.GENERIC + old LMDZ core |
|---|
| 4 | ######### AND LMDZ.GENERIC + new LMDZ core |
|---|
| 5 | ######### |
|---|
| 6 | |
|---|
| 7 | ### A INCLURE ? |
|---|
| 8 | ## calend=venus ## calendrier specifique |
|---|
| 9 | ## resetvarc=n ## reinit des variables de controle |
|---|
| 10 | ## iflag_trac=0 ## avec ou sans traceurs |
|---|
| 11 | |
|---|
| 12 | ################################################################ |
|---|
| 13 | #### LMDZ5 |
|---|
| 14 | #### IN GCM.DEF |
|---|
| 15 | #### -- those parameters are ignored when using old LMDZ core |
|---|
| 16 | ################################################################ |
|---|
| 17 | |
|---|
| 18 | ## period appel a la dissipation (multiple iperiod, 0 si calcul automatique) |
|---|
| 19 | ## (remplace ce qui etait idissip dans l'ancienne dynamique) |
|---|
| 20 | dissip_period = 0 |
|---|
| 21 | dissip_period = 2 |
|---|
| 22 | dissip_period = 5 |
|---|
| 23 | |
|---|
| 24 | ## avec ou sans physique |
|---|
| 25 | ## 0: pas de physique (e.g. en mode Shallow Water) |
|---|
| 26 | ## 1: avec physique (e.g. physique phylmd) |
|---|
| 27 | ## 2: avec rappel newtonien dans la dynamique |
|---|
| 28 | iflag_phys=2 |
|---|
| 29 | iflag_phys=1 |
|---|
| 30 | |
|---|
| 31 | ## avec ou sans fichiers de demarrage (start.nc, startphy.nc) ? |
|---|
| 32 | ## (sans fichiers de demarrage, initialisation des champs par iniacademic |
|---|
| 33 | ## dans la dynamique) |
|---|
| 34 | read_start=n |
|---|
| 35 | read_start=y |
|---|
| 36 | |
|---|
| 37 | ## Avec ou sans strato // i.e. Couche eponge et second palier pour dissip horizontale |
|---|
| 38 | ok_strato=n |
|---|
| 39 | |
|---|
| 40 | ### Dissipation horizontale |
|---|
| 41 | dissip_fac_mid=1. # for no var dissip (see inidissip) |
|---|
| 42 | dissip_fac_up=1. # for no var dissip (see inidissip) |
|---|
| 43 | ## deltaz et hdelta en km |
|---|
| 44 | #dissip_deltaz=10. |
|---|
| 45 | #dissip_hdelta=5. |
|---|
| 46 | ## pupstart en Pa |
|---|
| 47 | #dissip_pupstart=1.e3 |
|---|
| 48 | |
|---|
| 49 | ######################################## |
|---|
| 50 | ## LMDZ.COMMON |
|---|
| 51 | ## iflag_top_bound=0 for no sponge |
|---|
| 52 | ## iflag_top_bound=1 for sponge over 4 topmost layers |
|---|
| 53 | ## iflag_top_bound=2 for sponge from top to ~1% of top layer pressure |
|---|
| 54 | iflag_top_bound=0 |
|---|
| 55 | iflag_top_bound=2 |
|---|
| 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=0 |
|---|
| 63 | mode_top_bound=2 |
|---|
| 64 | mode_top_bound=1 |
|---|
| 65 | ####################################### |
|---|
| 66 | |
|---|
| 67 | # Coefficient pour la couche eponge (valeur derniere couche) |
|---|
| 68 | tau_top_bound=1.e-5 |
|---|
| 69 | #tau_top_bound=5.e-5 |
|---|
| 70 | #tau_top_bound=1.e-4 |
|---|
| 71 | |
|---|
| 72 | ################################################################ |
|---|
| 73 | #### LMDZ5 |
|---|
| 74 | #### IN RUN.DEF |
|---|
| 75 | #### -- those parameters are ignored when using old LMDZ core |
|---|
| 76 | ################################################################ |
|---|
| 77 | |
|---|
| 78 | ## Type de calendrier utilise |
|---|
| 79 | ## valeur possible: earth_360d (defaut), earth_365d, earth_366d |
|---|
| 80 | ip_ebil_phy=1 |
|---|
| 81 | ip_ebil_dyn=1 |
|---|
| 82 | calend=earth_360d |
|---|
| 83 | |
|---|
| 84 | ### Jour de l'etat initial ( = 350 si 20 Decembre ,par expl. ,comme ici ) |
|---|
| 85 | #dayref=10120 |
|---|
| 86 | ### Annee de l'etat initial ( avec 4 chiffres ) |
|---|
| 87 | #anneeref=2000 |
|---|
| 88 | |
|---|
| 89 | ## Remise a zero de la date initiale |
|---|
| 90 | raz_date=0 |
|---|
| 91 | ### for restart runs |
|---|
| 92 | |
|---|
| 93 | ## flag de sortie dynzon |
|---|
| 94 | ok_dynzon=n |
|---|
| 95 | ## activation du calcul d equilibrage de charge |
|---|
| 96 | adjust=n |
|---|
| 97 | ## activation du filtre fft |
|---|
| 98 | use_filtre_fft=n |
|---|
| 99 | #use_filtre_fft=y |
|---|
| 100 | ## niveau d'impression de controle |
|---|
| 101 | prt_level=10 |
|---|
| 102 | prt_level=1 |
|---|
| 103 | |
|---|
| 104 | ## coordonnees type planeto |
|---|
| 105 | disvert_type = 2 |
|---|
| 106 | |
|---|
| 107 | ## the generic physics |
|---|
| 108 | planet_type = generic |
|---|
| 109 | |
|---|
| 110 | ## sorties instantanees dans la dynamique (fichiers dyn_hist.nc and co.) |
|---|
| 111 | ok_dyn_ins=n |
|---|
| 112 | #ok_dyn_ins=y |
|---|
| 113 | ## periode d'ecriture des sorties instantanees dans la dynamique |
|---|
| 114 | ## (en pas dynamiques) |
|---|
| 115 | iecri=30 |
|---|
| 116 | iecri=120 |
|---|
| 117 | iecri=480 |
|---|
| 118 | #iecri=4800 |
|---|
| 119 | #iecri=9600 |
|---|
| 120 | |
|---|
| 121 | ## sorties de valeurs moyennes dans la dynamique (fichiers dyn_hist_ave.nc and co.) |
|---|
| 122 | ok_dyn_ave=n |
|---|
| 123 | ## periode de stockage des moyennes dans la dynamique (en jours) |
|---|
| 124 | periodav=1 |
|---|
| 125 | |
|---|
| 126 | |
|---|
| 127 | ok_guide = y |
|---|
| 128 | facwind = 5.e7 |
|---|
| 129 | #facwind = 1.e6 |
|---|
| 130 | facwind = 2.e7 |
|---|
| 131 | |
|---|
| 132 | ################################# |
|---|
| 133 | #### AND NOW THE OLD run.def #### |
|---|
| 134 | ################################# |
|---|
| 135 | ## -----> REMOVED |
|---|
| 136 | ## was present and conflicting with LMDZ5 |
|---|
| 137 | ## (although useless in old OLD LMDz DYN CORE |
|---|
| 138 | ## - iecri [periode d'ecriture du fichier histoire (en jour)] |
|---|
| 139 | ## - periodav [periode de stockage fichier histmoy (en jour) dynzon?] |
|---|
| 140 | ################################# |
|---|
| 141 | ################################# |
|---|
| 142 | |
|---|
| 143 | #----------------------------------------------------------------------- |
|---|
| 144 | #Parametres de controle du run: |
|---|
| 145 | #------------------------------ |
|---|
| 146 | |
|---|
| 147 | # Nombre de jours d'integration |
|---|
| 148 | nday = 720 |
|---|
| 149 | nday = 20 |
|---|
| 150 | nday = 100 |
|---|
| 151 | nday = 1800 |
|---|
| 152 | nday = 250000 |
|---|
| 153 | nday = 65000 |
|---|
| 154 | |
|---|
| 155 | # nombre de pas par jour (multiple de iperiod) ( ici pour dt = 1 min ) |
|---|
| 156 | # day_step = 25 |
|---|
| 157 | # day_step = 50 |
|---|
| 158 | day_step = 70 |
|---|
| 159 | day_step = 140 |
|---|
| 160 | day_step = 280 |
|---|
| 161 | # day_step = 560 |
|---|
| 162 | # day_step = 240 |
|---|
| 163 | # day_step = 480 |
|---|
| 164 | # day_step = 960 |
|---|
| 165 | |
|---|
| 166 | # periode pour le pas Matsuno (en pas) |
|---|
| 167 | iperiod=5 |
|---|
| 168 | |
|---|
| 169 | # periode de sortie des variables de controle (en pas) |
|---|
| 170 | iconser=120 |
|---|
| 171 | |
|---|
| 172 | # periode de la dissipation (en pas) |
|---|
| 173 | # -- this is ignored by LMDZ5 which rely on dissip_period instead |
|---|
| 174 | idissip=5 |
|---|
| 175 | |
|---|
| 176 | # choix de l'operateur de dissipation (star ou non star ) |
|---|
| 177 | lstardis=.true. |
|---|
| 178 | |
|---|
| 179 | # avec ou sans coordonnee hybrides |
|---|
| 180 | hybrid=.false. # calculs 1D sandrine #.true. |
|---|
| 181 | |
|---|
| 182 | # nombre d'iterations de l'operateur de dissipation gradiv |
|---|
| 183 | nitergdiv=1 |
|---|
| 184 | |
|---|
| 185 | # nombre d'iterations de l'operateur de dissipation nxgradrot |
|---|
| 186 | nitergrot=2 |
|---|
| 187 | |
|---|
| 188 | # nombre d'iterations de l'operateur de dissipation divgrad |
|---|
| 189 | niterh=2 |
|---|
| 190 | |
|---|
| 191 | # temps de dissipation des plus petites long.d ondes pour u,v (gradiv) |
|---|
| 192 | tetagdiv=3600. 8000. 160000. 8000. 1000. 2500. 4000. |
|---|
| 193 | |
|---|
| 194 | # temps de dissipation des plus petites long.d ondes pour u,v(nxgradrot) |
|---|
| 195 | tetagrot=5400. 8000. 160000. 8000. 1000. 2500. 4000. |
|---|
| 196 | |
|---|
| 197 | # temps de dissipation des plus petites long.d ondes pour h ( divgrad) |
|---|
| 198 | tetatemp=5400. 8000. 160000. 8000. 1000. 2500. 4000. |
|---|
| 199 | |
|---|
| 200 | # coefficient pour gamdissip |
|---|
| 201 | coefdis=0. |
|---|
| 202 | |
|---|
| 203 | # choix du shema d'integration temporelle (Matsuno ou Matsuno-leapfrog) |
|---|
| 204 | purmats=.false. |
|---|
| 205 | |
|---|
| 206 | # avec ou sans physique |
|---|
| 207 | physic=.true. |
|---|
| 208 | |
|---|
| 209 | # periode de la physique (en pas) (en pas dynamiques, n'a de sens que si iflag_phys=1) |
|---|
| 210 | iphysiq=5 |
|---|
| 211 | iphysiq=7 |
|---|
| 212 | # iphysiq=1 |
|---|
| 213 | |
|---|
| 214 | # choix d'une grille reguliere |
|---|
| 215 | grireg=.true. |
|---|
| 216 | |
|---|
| 217 | # frequence (en pas) de l'ecriture du fichier diagfi |
|---|
| 218 | ecritphy=500 |
|---|
| 219 | ecritphy=7000 |
|---|
| 220 | |
|---|
| 221 | # longitude en degres du centre du zoom |
|---|
| 222 | clon=0. |
|---|
| 223 | |
|---|
| 224 | # latitude en degres du centre du zoom |
|---|
| 225 | clat=0. |
|---|
| 226 | |
|---|
| 227 | # facteur de grossissement du zoom,selon longitude |
|---|
| 228 | grossismx=1. |
|---|
| 229 | |
|---|
| 230 | # facteur de grossissement du zoom ,selon latitude |
|---|
| 231 | grossismy=1. |
|---|
| 232 | |
|---|
| 233 | # Fonction f(y) hyperbolique si = .true. , sinon sinusoidale |
|---|
| 234 | fxyhypb=.false. |
|---|
| 235 | |
|---|
| 236 | # extension en longitude de la zone du zoom ( fraction de la zone totale) |
|---|
| 237 | dzoomx= 0. |
|---|
| 238 | |
|---|
| 239 | # extension en latitude de la zone du zoom ( fraction de la zone totale) |
|---|
| 240 | dzoomy=0. |
|---|
| 241 | |
|---|
| 242 | # raideur du zoom en X |
|---|
| 243 | taux=2. |
|---|
| 244 | |
|---|
| 245 | # raideur du zoom en Y |
|---|
| 246 | tauy=2. |
|---|
| 247 | |
|---|
| 248 | # Fonction f(y) avec y = Sin(latit.) si = .TRUE. , Sinon y = latit. |
|---|
| 249 | ysinus= .false. |
|---|
| 250 | |
|---|
| 251 | ###################### SPONGE LAYER IN OLD DYNAMICS |
|---|
| 252 | ###################### -- this is ignored by LMDZ5 !!! NO LONGER !!! |
|---|
| 253 | ## Avec sponge layer |
|---|
| 254 | # callsponge = .true. |
|---|
| 255 | ## Sponge: mode0(u=v=0), mode1(u=umoy,v=0), mode2(u=umoy,v=vmoy) |
|---|
| 256 | # mode_sponge= 2 |
|---|
| 257 | ## Sponge: hauteur de sponge (km) |
|---|
| 258 | # hsponge= 50 |
|---|
| 259 | ## Sponge: tetasponge (secondes) |
|---|
| 260 | # tetasponge = 50000 |
|---|
| 261 | ####################### |
|---|
| 262 | |
|---|
| 263 | # some definitions for the physics, in file 'callphys.def' |
|---|
| 264 | INCLUDEDEF=callphys.def |
|---|
