1 | |
---|
2 | #---------------- Planet ---------------- |
---|
3 | |
---|
4 | # Planetary radius : real (default=6371220) |
---|
5 | radius = 6.37122e6 |
---|
6 | |
---|
7 | # Gravity : real (default=9.80616) |
---|
8 | g = 9.80616 |
---|
9 | |
---|
10 | # kappa=Rd/cpp : real (default=0.2857143) |
---|
11 | kappa = 0.2857143 |
---|
12 | |
---|
13 | # Heat capacity at constant pressure : real (default=1004.70885) |
---|
14 | cpp = 1004.5 |
---|
15 | |
---|
16 | # Reference pressure : real (default=101325) |
---|
17 | preff = 1e5 |
---|
18 | |
---|
19 | # Planetary rotation rate : real (default=7.292E-5) |
---|
20 | omega = 0 |
---|
21 | |
---|
22 | |
---|
23 | #---------------- Constants ---------------- |
---|
24 | |
---|
25 | # Reference temperature T0 for DCMIP test cases : real (default=300) |
---|
26 | ncar_T0 = 300 |
---|
27 | |
---|
28 | # Reference pressure p0 for DCMIP test cases, used by disvert if disver=ncar : real (default=1e5) |
---|
29 | ncar_p0 = 1e5 |
---|
30 | |
---|
31 | # Exponent for B(eta), used by disvert if disvert=ncar : integer (default=1) |
---|
32 | ncar_disvert_c = 1 |
---|
33 | |
---|
34 | |
---|
35 | #---------------- Mesh ---------------- |
---|
36 | |
---|
37 | # Number of subdivisions on a main triangle : integer (default=40) |
---|
38 | nbp = 64 |
---|
39 | |
---|
40 | # Number of vertical layers : integer (default=19) |
---|
41 | llm = 15 |
---|
42 | |
---|
43 | # Vertical grid : [std|ncar|ncarl30] (default=std) |
---|
44 | disvert = ncar |
---|
45 | |
---|
46 | # Model layer thickness in m, used if disvert=ncar : real (default=300) |
---|
47 | ncar_dz = 800 |
---|
48 | |
---|
49 | # Mesh optimisation : number of iterations : integer (default=0) |
---|
50 | optim_it = 100 |
---|
51 | |
---|
52 | # Sub splitting of main rhombus : integer (default=1) |
---|
53 | nsplit_i = 2 |
---|
54 | nsplit_j = 2 |
---|
55 | |
---|
56 | |
---|
57 | #---------------- Numerics ---------------- |
---|
58 | |
---|
59 | # Time step in s : real (default=480) |
---|
60 | dt = 300 |
---|
61 | |
---|
62 | # Dissipation time for grad(div) : real (default=5000) |
---|
63 | tau_graddiv = 10800 |
---|
64 | |
---|
65 | # Exponent of grad(div) disspation : integer (default=1) |
---|
66 | nitergdiv = 2 |
---|
67 | |
---|
68 | # Dissipation time for curl(curl) : real (default=5000) |
---|
69 | tau_gradrot = 10800 |
---|
70 | |
---|
71 | # Exponent of curl(curl) disspation : integer (default=1) |
---|
72 | nitergrot = 2 |
---|
73 | |
---|
74 | # Dissipation time for div(grad) : real (default=5000) |
---|
75 | tau_divgrad = 10800 |
---|
76 | |
---|
77 | # Exponent of div(grad) disspation : integer (default=1) |
---|
78 | niterdivgrad = 2 |
---|
79 | |
---|
80 | |
---|
81 | #---------------- Time ---------------- |
---|
82 | |
---|
83 | # Run length in s : real (default=??) |
---|
84 | run_length = 518400 |
---|
85 | |
---|
86 | # Interval in s between two outputs : integer (default=??) |
---|
87 | write_period = 21600 |
---|
88 | |
---|
89 | |
---|
90 | #---------------- Physical parameters ---------------- |
---|
91 | |
---|
92 | # Scale factor for small planet experiment : real (default=1.) |
---|
93 | scale_factor = 1 |
---|
94 | |
---|
95 | # Initial state : |
---|
96 | # [jablonowsky06|academic|dcmip[1-4]|heldsz|dcmip2_schaer_noshear] (default=jablonowsky06) |
---|
97 | etat0 = dcmip2_mountain |
---|
98 | |
---|