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. |
---|