- Timestamp:
- Jan 30, 2014, 12:16:55 PM (11 years ago)
- Location:
- trunk/LMDZ.UNIVERSAL
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.UNIVERSAL/README.aymeric
r1105 r1171 1 2 3 makelmdz_fcm -s 1 -d 64x48x64 -b 20x30 -arch GNOMEp -parallel mpi -p generic gcm 4 makelmdz_fcm -s 1 -d 128x96x64 -b 20x30 -arch GNOMEp -parallel mpi -p generic gcm 5 ##################################### 6 ## REF SANDRINE 7 makelmdz_fcm -cpp NODYN -s 1 -d 8x128x64 -b 20x30 -arch GNOMEp -parallel mpi -p generic gcm 8 ##################################### 9 10 11 12 13 14 makelmdz_fcm -s 1 -d 64x48x64 -b 17x23 -arch GNOMEp -parallel mpi -p generic gcm 15 16 test rings 17 makelmdz_fcm -cpp NODYN -s 1 -d 8x48x32 -b 17x23 -arch GNOMEp -parallel mpi -p generic gcm 18 makelmdz_fcm -cpp NODYN -s 1 -d 8x48x32 -b 20x30 -arch GNOMEp -parallel mpi -p generic gcm 19 20 ----------------------------------- 1 21 2 22 -
trunk/LMDZ.UNIVERSAL/deftank/PROTO_Saturn_Article/run.def
r1165 r1171 1 ########################################################################### 2 ### RUN.DEF ### RUN.DEF ### RUN.DEF ### RUN.DEF ### RUN.DEF ### RUN.DEF ### 3 ########################################################################### 4 ### A MERGED VERSION TO BE USED WITH OLD OR RECENT LMDZ CORE 5 ### i.e. compiles gcm.def+run.def in LMDZ5 with run.def in old planeto LMDZ 6 ### (*) indicates same parameters in old and recent LMDZ core 7 ########################################################################### 1 8 2 ######### 3 ######### THIS IS INTENDED TO WORK BOTH WITH LMDZ.GENERIC + old LMDZ core 4 ######### AND LMDZ.GENERIC + new LMDZ core 5 ######### 9 ########################################################################### 10 ### NUMERICAL INTEGRATIONS 11 ########################################################################### 12 ### nday --> number of days of integration (*) 13 ### iperiod --> period for the Matsuno step (in dynamical steps) (*) 14 ### day_step --> number of steps per day (multiple of iperiod) (*) 15 ### purmats --> temporal integration: Matsuno or Matsuno-leapfrog? (*) 16 ### grireg --> regular grid or not (*) 17 ### adjust --> flag to activate load optimization 18 ### use_filtre_fft --> flag to activate FFT filtering 19 ########################################################################### 20 nday = 65000 21 iperiod = 5 22 day_step = 70 23 purmats = .false. 24 grireg = .true. 25 adjust = n 26 use_filtre_fft = n 6 27 7 ### A INCLURE ? 8 ## calend=venus ## calendrier specifique 9 ## resetvarc=n ## reinit des variables de controle 10 ## iflag_trac=0 ## avec ou sans traceurs 28 ########################################################################### 29 ### WHAT RELATES TO PHYSICS 30 ########################################################################### 31 ### iphysiq --> period for calling physics (in dynamical steps) (*) 32 ### ecritphy --> frequency for outputs in diagfi (in dynamical steps) (*) 33 ### planet_type --> the planet 34 ### disvert_type --> type of vertical coordinates (2 for planeto) 35 ### hybrid --> use hybrid coordinates or not? (*) 36 ### iflag_phys=0 --> no physics (e.g. shallow water) 37 ### iflag_phys=1 --> physics (e.g. phylmd) 38 ### iflag_phys=2 --> simple newtonian relaxation 39 ### physic --> flag with or without physics (*) 40 ########################################################################### 41 iphysiq = 14 42 ecritphy = 7000 43 planet_type = generic 44 disvert_type = 2 45 hybrid = .false. 46 iflag_phys = 1 47 physic = .true. 11 48 12 ################################################################ 13 #### LMDZ5 14 #### IN GCM.DEF 15 #### -- those parameters are ignored when using old LMDZ core 16 ################################################################ 49 ########################################################################### 50 ### INITIALIZATION and NUDGING 51 ########################################################################### 52 ### read_start --> with or without start files ? 53 ### (without: initialization by iniacademic) 54 ### raz_date=0 --> use the date in start file 55 ### raz_date=1 --> erase the date and start from scratch 56 ### ok_guide --> flag for nudging towards winds (needs relax.nc) 57 ### facwind --> time constant for relaxation 58 ### 5.e8 medium 1.e7-5.e6 strong 1.e6 very fast 59 ########################################################################### 60 read_start = y 61 raz_date = 0 62 ok_guide = n 63 facwind = 1.e6 17 64 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 65 ########################################################################### 66 ### HORIZONTAL DISSIPATION 67 ########################################################################### 68 ### dissip_period --> period call to dissipation (in dynamical steps) 69 ### - must be a multiple of iperiod 70 ### - 0 means automatic computation 71 ### - use idissip instead if old LMDZ core (idissip ignored in LMDZ5) 72 ### dissip_fac_mid --> enhancement factor in middle atmosphere 73 ### - set to 1 to prevent dissipation enhancement 74 ### lstardis --> star or not star dissipation operator (*) 75 ### nitergdiv --> number of iterations for gradiv operator (*) 76 ### nitergrot --> number of iterations for nxgradrot operator (*) 77 ### niterh --> number of iterations for divgrad operator (*) 78 ### tetagdiv --> dissipation timescale of smallest wvl: u,v (gradiv) (*) 79 ### tetagrot --> dissipation timescale of smallest wvl: u,v (nxgradrot) (*) 80 ### tetatemp --> dissipation timescale of smallest wvl: h (divgrad) (*) 81 ### coefdis --> coefficient for gamdissip (*) 82 ########################################################################### 22 83 dissip_period = 5 84 dissip_fac_mid = 1. 85 lstardis = .true. 86 nitergdiv = 1 87 nitergrot = 2 88 niterh = 2 89 tetagdiv = 5400. 8000. 160000. 1000. 2500. 4000. 90 tetagrot = 7200. 8000. 160000. 1000. 2500. 4000. 91 tetatemp = 7200. 8000. 160000. 1000. 2500. 4000. 23 92 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 93 ########################################################################### 94 ### ENHANCEMENT OF HORIZONTAL DISSIPATION IN UPPERMOST LEVELS 95 ########################################################################### 96 ### - not activated if ok_strato=n 97 ### - dissip_fac_up --> enhancement factor 98 ### - dissip_deltaz --> altitude range of transition (km) 99 ### - dissip_hdelta --> scale height at transition (km) 100 ### - dissip_pupstart --> pressure at bottom of transition region (Pa) 101 ########################################################################### 102 ok_strato = n 103 dissip_fac_up = 5 104 dissip_deltaz = 150. 105 dissip_hdelta = 50. 106 dissip_pupstart = 50. 30 107 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 108 ########################################################################### 109 ### DAMPING ("SPONGE") LAYER AT MODEL TOP 110 ########################################################################### 111 ### iflag_top_bound=0 for no sponge 112 ### iflag_top_bound=1 for sponge over 4 topmost layers 113 ### iflag_top_bound=2 for sponge from top to ~1% of top layer pressure 114 ### mode_top_bound=0 for no sponge 115 ### mode_top_bound=1 for u,v --> 0 116 ### mode_top_bound=2 for u,v --> zonal mean 117 ### mode_top_bound=3 for u,v,h --> zonal mean 118 ### tau_top_bound --> damping frequency at last layer 119 ### e-5 medium / e-4 strong yet reasonable / e-3 very strong 120 ########################################################################### 121 iflag_top_bound = 1 122 mode_top_bound = 1 123 tau_top_bound = 1.e-4 36 124 37 ## Avec ou sans strato // i.e. Couche eponge et second palier pour dissip horizontale 38 ok_strato=n 125 ########################################################################### 126 ### OUTPUT FROM DYNAMICS 127 ########################################################################### 128 ### prt_level --> level of control prints (1-10) 129 ### ok_dyn_ins --> flag instantaneous outputs (dyn_hist.nc) 130 ### iecri --> frequency of instantaneous outputs (in dynamical steps) 131 ### ok_dyn_ave --> flag averaged outputs (dyn_hist_ave.nc) 132 ### periodav --> period of output of averages in dynamics (in days) 133 ### ok_dynzon --> flag outputs dynzon 134 ### NB: iecri and periodav used to be in old LMDZ core although useless 135 ### iconser --> frequency of outputs of control variables (in dyn steps)(*) 136 ########################################################################### 137 prt_level = 1 138 ok_dyn_ins = n 139 iecri = 480 140 ok_dyn_ave = n 141 periodav = 1 142 ok_dynzon = 0 143 iconser = 120 39 144 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 145 ########################################################################### 146 ### ZOOMED GRID (same parameters in old and recent LMDZ core) 147 ########################################################################### 148 ### clon --> longitude of center of zoom (deg) (*) 149 ### clat --> latitude of center of zoom (deg) (*) 150 ### grossismx --> zoom factor in longitude (*) 151 ### grossismy --> zoom factor in latitude (*) 152 ### dzoomx --> longitude extent of zoom zone (fraction of total area) (*) 153 ### dzoomy --> latitude extent of zoom zone (fraction of total area) (*) 154 ### taux --> "stifness" of zoom in longitude (*) 155 ### tauy --> "stifness" of zoom in latitude (*) 156 ### fxyhypb --> f(y) hyperbolic (T) or sinusoidal (F) (*) 157 ### ysinus --> f(y) = sin(lat) (T) or lat (F) (*) 158 ########################################################################### 159 clon = 0. 160 clat = 0. 161 grossismx = 1. 162 grossismy = 1. 163 dzoomx = 0. 164 dzoomy = 0. 165 taux = 2. 166 tauy = 2. 167 fxyhypb = .false. 168 ysinus = .false. 48 169 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 68 # Coefficient pour la couche eponge (valeur derniere couche) 69 tau_top_bound=1.e-5 70 #tau_top_bound=5.e-5 71 #tau_top_bound=1.e-4 72 73 ## guidage par vent dans relax.nc 74 ok_guide = y 75 ok_guide = n 76 facwind = 2.e7 77 facwind = 5.e8 78 79 ################################################################ 80 #### LMDZ5 81 #### IN RUN.DEF 82 #### -- those parameters are ignored when using old LMDZ core 83 ################################################################ 84 85 ## Type de calendrier utilise 86 ## valeur possible: earth_360d (defaut), earth_365d, earth_366d 87 ip_ebil_phy=1 88 ip_ebil_dyn=1 89 calend=earth_360d 90 91 ### Jour de l'etat initial ( = 350 si 20 Decembre ,par expl. ,comme ici ) 92 #dayref=10120 93 ### Annee de l'etat initial ( avec 4 chiffres ) 94 #anneeref=2000 95 96 ## Remise a zero de la date initiale 97 raz_date=0 98 ### for restart runs 99 # use 1 for 'from scratch' runs 100 101 ## flag de sortie dynzon 102 ok_dynzon=n 103 ## activation du calcul d equilibrage de charge 104 adjust=n 105 ## activation du filtre fft 106 use_filtre_fft=n 107 ## niveau d'impression de controle 108 prt_level=10 109 prt_level=1 110 111 ## coordonnees type planeto 112 disvert_type = 2 113 114 ## the generic physics 115 planet_type = generic 116 117 ## sorties instantanees dans la dynamique (fichiers dyn_hist.nc and co.) 118 ok_dyn_ins=n 119 #ok_dyn_ins=y 120 ## periode d'ecriture des sorties instantanees dans la dynamique 121 ## (en pas dynamiques) 122 #iecri=30 123 #iecri=120 124 iecri=480 125 #iecri=960 126 #iecri=4800 127 #iecri=9600 128 129 ## sorties de valeurs moyennes dans la dynamique (fichiers dyn_hist_ave.nc and co.) 130 ok_dyn_ave=n 131 ## periode de stockage des moyennes dans la dynamique (en jours) 132 periodav=1 133 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 = 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= 5400. 8000. 160000. 1000. 2500. 4000. 193 194 # temps de dissipation des plus petites long.d ondes pour u,v(nxgradrot) 195 tetagrot= 7200. 8000. 160000. 1000. 2500. 4000. 196 197 # temps de dissipation des plus petites long.d ondes pour h ( divgrad) 198 tetatemp= 7200. 8000. 160000. 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=14 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 !!! 170 ########################################################################### 171 ### LEFTOVERS 172 ### calend=earth_360d ## specific calendar 173 ### resetvarc=n ## reinit des variables de controle 174 ### iflag_trac=0 ## avec ou sans traceurs 175 ### dayref=10120 ## initial date 176 ### anneeref=2000 ## year of initial state 177 ### --- old version of spong layer 253 178 ## Avec sponge layer 254 179 # callsponge = .true. … … 259 184 ## Sponge: tetasponge (secondes) 260 185 # tetasponge = 50000 261 #######################262 186 263 # some definitions for the physics, in file 'callphys.def' 187 ########################################################################### 188 ### INCLUDE OTHER DEF FILES (physics, specific settings, etc...) 189 ########################################################################### 264 190 INCLUDEDEF=callphys.def 265 #INCLUDEDEF=modif.def191 INCLUDEDEF=modif.def
Note: See TracChangeset
for help on using the changeset viewer.