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