- Timestamp:
- Apr 14, 2023, 5:49:58 PM (21 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.MARS/util/startarchive2icosa/run_icosa.def_no_start
r2813 r2941 1 #INCLUDEDEF=saturn_const.def2 1 3 ##mpi_threading_mode=serialized 4 #mpi_threading_mode=funneled 2 #---------------- Mesh ---------------- 5 3 6 #halo_i=0 7 #halo_j=1 8 # ------------------------------- Mesh --------------------------------- 4 # Number of subdivisions on a main triangle : integer (default=40) 5 nbp = 60 9 6 10 # Number of subdivision on a main triangle (nbp) : integer (default=40) 11 #nbp=160 12 nbp=20 7 # Number of vertical layers : integer (default=19) 8 llm = 54 13 9 14 # Earth as reference ? 15 # nbp 20 40 80 160 16 # T-edge length (km) 500 250 120 60 17 18 ## sub splitting of main rhombus : integer (default=1) 19 #nsplit_i=1 20 #nsplit_j=1 21 #omp_level_size=1 22 ########################################### 23 ## There must be less MPIxOpenMP processes than the 10 x nsplit_i x nsplit_j tiles 24 ## typically for pure MPI runs, let nproc = 10 x nsplit_i x nsplit_j 25 ## it is better to have nbp/split >~ 10 26 ########################################### 27 ## 10 procs 28 nsplit_i=1 29 nsplit_j=2 30 #### 50 noeuds de 24 processeurs = 1200 procs 31 #nsplit_i=10 32 #nsplit_j=12 33 34 # Number of vertical layer (llm) : integer (default=19) 35 #llm=64 36 llm=73 37 38 # disvert : vertical discretisation : string (default='std') : std, ncar, ncar30l 10 # Vertical grid : [std|ncar|ncarl30] (default=std) 39 11 #disvert=read_apbp 40 12 disvert=plugin 41 13 42 # optim_it : mesh optimisation : number of iteration : integer (default=0) 43 #optim_it=1000 44 optim_it=100 14 # Mesh optimisation : number of iterations : integer (default=0) 15 optim_it = 1000 45 16 46 # ---------------------------------- Time--------------------------------- 17 # Sub splitting of main rhombus : integer (default=1) 18 nsplit_i = 1 19 nsplit_j = 1 47 20 48 ## scheme type : string : euler, leapfrog_matsuno, runge_kutta ) 49 # time_scheme = ARK2.3 (RK 2nd order at 3 levels) | leapfrog_matsuno | ... 50 time_scheme = ARK2.3 21 #number of openmp task on vertical level 22 omp_level_size=1 51 23 52 ## matsuno period : integer ( default=5) Only used for leapfrog_matsuno scheme 53 matsuno_period = 10 24 #---------------- Numerics ---------------- 54 25 55 # timestep : real (default=480s) 56 # dt = 92.473958333334 57 ### attention 3 chiffres significatifs dans XIOS 58 # number of dynamical steps per day (instead of dt) 59 day_step=960 26 # Advection called every itau_adv time steps : integer (default=2) 27 itau_adv = 1 60 28 61 # advection called every itau_adv time steps : integer (default 2)62 # standard : umax=100m/s vs c=340m/s (ratio 1:3)63 # in JW06 umax=35m/s vs c=340m/s (ratio 1:10)64 itau_adv=3 29 # Time step in s : real (default=480) 30 # dt = 720 31 # Alternative to specifying "dt", specify number of steps per day : day_step 32 day_step = 960 65 33 66 # number of timestep (default 100) (ignored if run_length is set) 67 #itaumax = 10 68 69 # run length : real (default=dt*itaumax) 70 71 #################################### 72 73 ## 1 martian day 74 #run_length = 88775 75 #run_length = 887750 76 77 ## 1 martian year : 59212925s 78 ## first martian month (61sols):5415275s 79 80 # number of days to run 81 ndays=1 82 83 #################################### 84 85 ##activate IO (default = true) 86 #enable_io = false 87 88 ## output with XIOS (only if compiled with XIOS): true/false (default true) 89 #xios_output=false 90 91 # output field period (valid XIOS or not) : integer (default none) 92 # [should be equal to dt if XIOS dynamical outputs are wanted] 93 # [because write_period sets synchronization frequency with XIOS] 94 write_period = 92.473958333 34 # Number of tracers : integer (default=1) 35 nqtot = 11 95 36 96 37 97 # itau_write_etat0=38052038 #---------------- Time and output ---------------- 98 39 99 # ---------------------------------- Misc -------------------------------- 40 # Time style : [none|dcmip] (default=dcmip) 41 time_style = none 100 42 101 # number of tracer (nqtot) : integer (default 1) 102 nqtot=1 43 # Run length in s : real (default=??) 44 # run_length = 103680000 45 run_length = 462.369791666 46 # Alternative to specifying "run_length", specify number of days to run : ndays 47 # ndays=1 103 48 104 # pression value where output is interpolated : real (default=0, no output) 105 ###out_pression_level=85000 49 # Interval in s between two outputs : integer (default=0) 50 #write_period = 10800 51 # output every dyn time step => dt = day_length/day_step ; 88775/960=92.4739583333 52 write_period= 92.4739583333333 106 53 107 # etat0 : initial state : string (default=jablonowsky06) : 108 # jablonowsky06, academic, ncar 54 55 # Initial state : 56 # [jablonowsky06|academic|dcmip[1-4]|heldsz|dcmip2_schaer_noshear] (default=jablonowsky06) 57 etat0=isothermal 58 etat0_isothermal_temp=200 59 # pertub initial state by a factor of ... 60 etat0_ps_white_noise=0.01 61 etat0_theta_rhodz_white_noise=0.01 109 62 110 63 #etat0=start_file 111 64 #etat0_start_file_colocated=true 112 65 113 ## -- to cross the 2y limit114 #etat0_start_iteration=0115 #run_length=38052116 117 etat0=isothermal118 etat0_isothermal_temp=200119 120 66 # start file name (default is start.nc) 121 67 # restart file name (default is restart.nc) 122 #start_file_name=start_icosa68 start_file_name=start_icosa 123 69 restart_file_name=restart_icosa 124 70 125 #etat0_start_iteration_reset=true126 #etat0_start_iteration=071 etat0_start_iteration_reset=true 72 etat0_start_iteration=0 127 73 128 ##### for one-day run from a profile to create a start 129 #etat0=temperature_profile 130 #temperature_profile_file=temp_profile.txt 131 #run_length=38052 74 # Dissipation time for grad(div) : real (default=5000) 75 tau_graddiv = 18000 132 76 133 # ------------------------------ Dynamics -------------------------------- 134 135 # caldyn : computation type for gcm equation : string (default=gcm) : gcm, adv 136 caldyn=gcm 137 138 # caldyn_conserv : string (default=energy) : energy,enstrophy 139 caldyn_conserv=energy 140 141 # caldyn_exner : scheme for computing Exner function : string (default=direct) : direct,lmdz 142 caldyn_exner=direct 143 144 # caldyn_hydrostat : scheme for computing geopotential : string (default=direct) : direct,lmdz 145 caldyn_hydrostat=direct 146 147 # guided_type : string (default=none) : none, ncar 148 guided_type=none 149 150 #Sponge layer 151 ### iflag_sponge=0 for no sponge 152 ### iflag_sponge=1 for sponge over 4 topmost layers 153 ### iflag_sponge=2 for sponge from top to ~1% of top layer pressure 154 ### tau_sponge --> damping frequency at last layer 155 ### e-5 medium / e-4 strong yet reasonable / e-3 very strong 156 ### mode_sponge=1 for u,v --> 0 157 ### mode_sponge=2 for u,v --> zonal mean (NOT IMPLEMENTED) 158 ### mode_sponge=3 for u,v,h --> zonal mean (NOT IMPLEMENTED) 159 #iflag_sponge = 1 160 iflag_sponge = 1 161 tau_sponge = 1.e-4 162 #mode_sponge = 1 163 164 165 # ---------------------------------- Dissipation --------------------------------- 166 167 # dissipation time graddiv : real (default=5000) 168 tau_graddiv = 2000 169 170 # number of iteration for graddiv : integer (default=1) 77 # Exponent of grad(div) disspation : integer (default=1) 171 78 nitergdiv = 2 172 79 173 # dissipation time nxgradrot(default=5000)174 tau_gradrot = 400080 # Dissipation time for curl(curl) : real (default=5000) 81 tau_gradrot = 18000 175 82 176 # number of iteration for nxgradrot: integer (default=1)177 nitergrot =283 # Exponent of curl(curl) disspation : integer (default=1) 84 nitergrot = 2 178 85 179 # dissipation time divgrad (theta) (default=5000) 180 tau_divgrad= 4000 181 # number of iteration for divgrad : integer (default=1) 182 niterdivgrad=2 86 # Dissipation time for div(grad) : real (default=5000) 87 tau_divgrad = 18000 183 88 184 # Vertical evolution of the dissipation coefficient 185 # vert_prof_dissip_type=legacy 186 vert_prof_dissip_type=plugin 89 # Exponent of div(grad) disspation : integer (default=1) 90 niterdivgrad = 2 187 91 188 # vertical dissipation setup (none|altitude|pressure) 189 vert_prof_dissip_setup=altitude 190 dissip_fac_mid=3. 191 dissip_fac_up=30. 192 dissip_startalt=70. 193 dissip_deltaz=30. 92 #---------------- Planet ---------------- 194 93 195 #Vertical coefficient for dissipation: 1 is for Earth, 2 Mars, 3 Venus, default is basic scheme 196 #vert_prof_dissip=2 197 #fac_mid=3. 198 #fac_up=30. 199 #dissip_deltaz=30. 200 #startalt=70. 201 202 # Rayleigh friction : string (default=none) : none, dcmip2_schaer_noshear, dcmip2_schaer_shear 203 rayleigh_friction_type=none 204 #rayleigh_friction_type=giant_liu_schneider 205 #rayleigh_limlat = 33. 206 #rayleigh_friction_tau = 8640000. 94 INCLUDEDEF=mars_const.def 207 95 208 96 # -------------------------------- Physics ------------------------------------- 209 97 210 # (itau_physics=160)*(dt=118.9125) --> half a saturn day as physical timestep 211 # ---- for some reason (day consistency, physical timestep must be less a day) 212 # ---- change also in run.def 98 # Physics package : [none|held_suarez|dcmip] (default=none) 99 physics = phys_external 100 101 # call physics every itau_physics (dynamical) time step 213 102 itau_physics=5 103 #itau_physics=10 214 104 215 # kind of physics : string : none, dcmip (default=none) 216 physics=phys_external 217 218 ## testcase physics for dcmip : INTEGER : 1, 2 (default=1) 219 #dcmip_physics=1 220 221 # -------------------- parameters for NCAR test cases ------------------------ 222 223 # NCAR advection test, initial tracer : string ( default='cos_bell') 224 # const, slotted_cyl, cos_bell, dbl_cos_bell_q1, dbl_cos_bell_q2, complement, hadley 225 ncar_adv_shape=cos_bell 226 227 # NCAR advection test, wind field : string (default='deform') : solid, deform, hadley 228 ncar_adv_wind=solid 229 230 # ncar_T0 : reference temperature for NCAR test cases : real (default=300) 231 # also used by disvert if disvert=ncar 232 ncar_T0=300 233 234 # ncar_p0 : reference pressure for NCAR test cases : real (default=1e5) 235 # also used by disvert if disvert=ncar 236 ncar_p0=1e5 237 238 # ncar_disvert_c : exponent for B(eta) : integer (default=1) 239 # used by disvert if disvert=ncar 240 241 ncar_disvert_c=1 242 243 # dcmip 4 testcase : integer (default=1) : 1, 2 244 dcmip4_testcase=1 245 #----------------------------------------------------------------------- 246 247 248 check_conservation = detailed 249 itau_check_conserv = 320 105 # physics initial state (startfi.nc) if startphy_file=.true. 106 startphy_file=.false.
Note: See TracChangeset
for help on using the changeset viewer.
