| 1 | |
|---|
| 2 | #---------------- Mesh ---------------- |
|---|
| 3 | |
|---|
| 4 | # Number of subdivisions on a main triangle : integer (default=40) |
|---|
| 5 | nbp = 60 |
|---|
| 6 | |
|---|
| 7 | # Number of vertical layers : integer (default=19) |
|---|
| 8 | llm = 54 |
|---|
| 9 | |
|---|
| 10 | # Vertical grid : [std|ncar|ncarl30] (default=std) |
|---|
| 11 | #disvert=read_apbp |
|---|
| 12 | disvert=plugin |
|---|
| 13 | |
|---|
| 14 | # Mesh optimisation : number of iterations : integer (default=0) |
|---|
| 15 | optim_it = 1000 |
|---|
| 16 | |
|---|
| 17 | # Sub splitting of main rhombus : integer (default=1) |
|---|
| 18 | nsplit_i = 1 |
|---|
| 19 | nsplit_j = 1 |
|---|
| 20 | |
|---|
| 21 | #number of openmp task on vertical level |
|---|
| 22 | omp_level_size=1 |
|---|
| 23 | |
|---|
| 24 | #---------------- Numerics ---------------- |
|---|
| 25 | |
|---|
| 26 | # Advection called every itau_adv time steps : integer (default=2) |
|---|
| 27 | itau_adv = 1 |
|---|
| 28 | |
|---|
| 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 |
|---|
| 33 | |
|---|
| 34 | # Number of tracers : integer (default=1) |
|---|
| 35 | nqtot = 11 |
|---|
| 36 | |
|---|
| 37 | |
|---|
| 38 | #---------------- Time and output ---------------- |
|---|
| 39 | |
|---|
| 40 | # Time style : [none|dcmip] (default=dcmip) |
|---|
| 41 | time_style = none |
|---|
| 42 | |
|---|
| 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 |
|---|
| 48 | |
|---|
| 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 |
|---|
| 53 | |
|---|
| 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 |
|---|
| 62 | |
|---|
| 63 | #etat0=start_file |
|---|
| 64 | #etat0_start_file_colocated=true |
|---|
| 65 | |
|---|
| 66 | # start file name (default is start.nc) |
|---|
| 67 | # restart file name (default is restart.nc) |
|---|
| 68 | start_file_name=start_icosa |
|---|
| 69 | restart_file_name=restart_icosa |
|---|
| 70 | |
|---|
| 71 | etat0_start_iteration_reset=true |
|---|
| 72 | etat0_start_iteration=0 |
|---|
| 73 | |
|---|
| 74 | # Dissipation time for grad(div) : real (default=5000) |
|---|
| 75 | tau_graddiv = 18000 |
|---|
| 76 | |
|---|
| 77 | # Exponent of grad(div) disspation : integer (default=1) |
|---|
| 78 | nitergdiv = 2 |
|---|
| 79 | |
|---|
| 80 | # Dissipation time for curl(curl) : real (default=5000) |
|---|
| 81 | tau_gradrot = 18000 |
|---|
| 82 | |
|---|
| 83 | # Exponent of curl(curl) disspation : integer (default=1) |
|---|
| 84 | nitergrot = 2 |
|---|
| 85 | |
|---|
| 86 | # Dissipation time for div(grad) : real (default=5000) |
|---|
| 87 | tau_divgrad = 18000 |
|---|
| 88 | |
|---|
| 89 | # Exponent of div(grad) disspation : integer (default=1) |
|---|
| 90 | niterdivgrad = 2 |
|---|
| 91 | |
|---|
| 92 | #---------------- Planet ---------------- |
|---|
| 93 | |
|---|
| 94 | INCLUDEDEF=mars_const.def |
|---|
| 95 | |
|---|
| 96 | # -------------------------------- Physics ------------------------------------- |
|---|
| 97 | |
|---|
| 98 | # Physics package : [none|held_suarez|dcmip] (default=none) |
|---|
| 99 | physics = phys_external |
|---|
| 100 | |
|---|
| 101 | # call physics every itau_physics (dynamical) time step |
|---|
| 102 | itau_physics=5 |
|---|
| 103 | #itau_physics=10 |
|---|
| 104 | |
|---|
| 105 | # physics initial state (startfi.nc) if startphy_file=.true. |
|---|
| 106 | startphy_file=.false. |
|---|