- Timestamp:
- 10/11/14 23:27:39 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/trunk/param_sets/climate/Held_Suarez/run.def
r170 r278 1 INCLUDEDEF=earth_const.def2 1 3 # Scale factor for small planet experiment : real (default=1.) 4 scale_factor=1. 2 #---------------- Mesh ---------------- 5 3 6 # ------------------------------- Mesh --------------------------------- 7 # nbp 20 40 80 160 8 # T-edge length (km) 500 250 120 60 4 # Number of subdivisions on a main triangle : integer (default=40) 5 nbp = 40 9 6 10 # Number of subdivision on a main triangle (nbp) : integer (default=40)11 nbp=40 7 # Number of vertical layers : integer (default=19) 8 llm = 19 12 9 13 # sub splitting of main rhombus : integer (default=1) 14 nsplit_i=2 15 nsplit_j=3 10 # Vertical grid : [std|ncar|ncarl30] (default=std) 11 disvert = std 16 12 17 # Number of vertical layer (llm) : integer (default=19)18 llm=19 13 # Mesh optimisation : number of iterations : integer (default=0) 14 optim_it = 1000 19 15 20 # disvert : vertical discretisation : string (default='std') : std, ncar, ncar30l 21 disvert=std 16 # Sub splitting of main rhombus : integer (default=1) 17 nsplit_i = 2 18 nsplit_j = 2 22 19 23 # optim_it : mesh optimisation : number of iteration : integer (default=0)24 optim_it=100025 20 26 # ---------------------------------- Time --------------------------------21 #---------------- Numerics ---------------- 27 22 28 # scheme type : string ( default='runge_kutta') euler, leapfrog_matsuno, runge_kutta)29 scheme = runge_kutta 23 # Advection called every itau_adv time steps : integer (default=2) 24 itau_adv = 1 30 25 31 # matsuno period : integer ( default=5)32 matsuno_period = 5 26 # Time step in s : real (default=480) 27 dt = 480 33 28 34 # advection called every itau_adv time steps : integer (default 2) 35 # standard : umax=100m/s vs c=340m/s (ratio 1:3) 36 # in JW06 umax=35m/s vs c=340m/s (ratio 1:10) 37 itau_adv= 1 29 # Number of tracers : integer (default=1) 30 nqtot = 0 38 31 39 # time_style40 time_style=lmd41 32 42 # Number of days 43 ndays=1200 33 #---------------- Time and output ---------------- 44 34 45 # NUBER of day_step46 day_step= 180 35 # Time style : [none|dcmip] (default=dcmip) 36 time_style = none 47 37 48 # output field period : integer (default none)49 write_period=8640038 # Run length in s : real (default=??) 39 run_length = 103680000 50 40 51 # ---------------------------------- Misc -------------------------------- 41 # Interval in s between two outputs : integer (default=??) 42 write_period = 86400 52 43 53 # number of tracer (nqtot) : integer (default 1)54 nqtot=155 44 56 # etat0 : initial state : string (default=jablonowsky06) : 57 # jablonowsky06, academic, dcmip[1-4] ,heldsz 58 etat0=held_suarez 45 #---------------- Planet ---------------- 59 46 60 # ------------------------------ Dynamics -------------------------------- 47 # Planetary radius : real (default=6371220) 48 radius = 6.371e6 61 49 62 # caldyn : computation type for gcm equation : string (default=gcm) : gcm, adv63 caldyn=gcm 50 # Gravity : real (default=9.80616) 51 g = 9.8 64 52 65 caldyn_conserv=energy 66 caldyn_exner=direct 67 caldyn_hydrostat=direct 53 # Planetary rotation rate : real (default=7.292E-5) 54 omega = 7.292e-5 68 55 69 # guided_type : string (default=none) : none, ncar70 guided_type=none 56 # kappa=Rd/cpp : real (default=0.2857143) 57 kappa = 0.2857143 71 58 72 # ------------------------------ Dissipation --------------------------- 59 # Heat capacity at constant pressure : real (default=1004.70885) 60 cpp = 1004 73 61 74 # dissipation time graddiv : real (default=5000) 62 # Reference pressure : real (default=101325) 63 preff = 1e5 64 65 66 #---------------- Physical parameters ---------------- 67 68 # Initial state : 69 # [jablonowsky06|academic|dcmip[1-4]|heldsz|dcmip2_schaer_noshear|dcmip3] (default=jablonowsky06) 70 etat0 = held_suarez 71 72 # Physics package : [none|held_suarez|dcmip] (default=none) 73 physics = held_suarez 74 75 # Dissipation time for grad(div) : real (default=5000) 75 76 tau_graddiv = 50000 76 77 77 # number of iteration for graddiv: integer (default=1)78 # Exponent of grad(div) disspation : integer (default=1) 78 79 nitergdiv = 2 79 80 80 # dissipation time nxgradrot(default=5000)81 # Dissipation time for curl(curl) : real (default=5000) 81 82 tau_gradrot = 50000 82 83 83 # number of iteration for nxgradrot: integer (default=1)84 nitergrot =284 # Exponent of curl(curl) disspation : integer (default=1) 85 nitergrot = 2 85 86 86 # dissipation time divgrad (theta)(default=5000)87 tau_divgrad = 5000087 # Dissipation time for div(grad) : real (default=5000) 88 tau_divgrad = 50000 88 89 89 # number of iteration for divgrad: integer (default=1)90 niterdivgrad =290 # Exponent of div(grad) disspation : integer (default=1) 91 niterdivgrad = 2 91 92 92 # Rayleigh friction : string (default=none) : none, dcmip2_schaer_noshear, dcmip2_schaer_shear93 rayleigh_friction_type=none94 95 #----------------------------physics-------------------------------------96 physics=automatic
Note: See TracChangeset
for help on using the changeset viewer.