[284] | 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 = 80 |
---|
| 39 | |
---|
| 40 | # Number of vertical layers : integer (default=19) |
---|
| 41 | llm = 10 |
---|
| 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 = 1000 |
---|
| 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 = 240 |
---|
| 61 | |
---|
| 62 | # Dissipation time for grad(div) : real (default=5000) |
---|
| 63 | tau_graddiv = 9000 |
---|
| 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 = 9000 |
---|
| 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 = 9000 |
---|
| 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 = 4.5e5 |
---|
| 85 | |
---|
| 86 | # Interval in s between two outputs : integer (default=??) |
---|
| 87 | write_period = 12500 |
---|
| 88 | |
---|
| 89 | |
---|
| 90 | #---------------- Physical parameters ---------------- |
---|
| 91 | |
---|
| 92 | # Scale factor for small planet experiment : real (default=1.) |
---|
| 93 | scale_factor = 125 |
---|
| 94 | |
---|
| 95 | # Initial state : |
---|
| 96 | # [jablonowsky06|academic|dcmip[1-4]|heldsz|dcmip2_schaer_noshear] (default=jablonowsky06) |
---|
| 97 | etat0 = dcmip3 |
---|
| 98 | |
---|