1 | #---------------------------------# |
---|
2 | # Run parameters for the 1D model # |
---|
3 | #---------------------------------# |
---|
4 | |
---|
5 | ## Planetary constants |
---|
6 | ## ~~~~~~~~~~~~~~~~~~~~~~~ |
---|
7 | ## NB: those are mandatory |
---|
8 | ## ~~~~~~~~~~~~~~~~~~~~~~~ |
---|
9 | ## LENGTH OF A DAY in s |
---|
10 | daysec = 887750. |
---|
11 | ## GRAVITY in m s-2 |
---|
12 | g = 26.0 |
---|
13 | ## LENGTH OF A YEAR in days |
---|
14 | year_day = 4332 |
---|
15 | ## MIN DIST STAR-PLANET in AU [periastron] |
---|
16 | periastr = 5.00 |
---|
17 | ## MAX DIST STAR-PLANET in AU [apoastron] |
---|
18 | apoastr = 5.00 |
---|
19 | ## DATE OF PERIASTRON in days |
---|
20 | peri_day = 170. |
---|
21 | ## OBLIQUITY in deg |
---|
22 | obliquit = 3.12 |
---|
23 | ## SURFACE PRESSURE in Pa |
---|
24 | psurf = 1000000 |
---|
25 | |
---|
26 | ## Time integration parameters |
---|
27 | ## ~~~~~~~~~~~~~~~~~~~~~~~~~~~ |
---|
28 | # Initial date (in solar days,=0 at Ls=0) |
---|
29 | day0 = 0 |
---|
30 | # Initial local time (in hours, between 0 and 24) |
---|
31 | time = 12 |
---|
32 | # Number of time steps per sol |
---|
33 | day_step = 1 |
---|
34 | # Number of sols to run |
---|
35 | ndt = 250000 |
---|
36 | |
---|
37 | ## Vertical levels |
---|
38 | ## ~~~~~~~~~~~~~~~ |
---|
39 | # hybrid vertical coordinate ? (.true. for hybrid and .false. for sigma levels) |
---|
40 | hybrid = .false. |
---|
41 | # autocompute vertical discretisation? (useful for exoplanet runs) |
---|
42 | autozlevs = .false. |
---|
43 | # Ceiling pressure (Pa) ? |
---|
44 | pceil = 10.0 |
---|
45 | |
---|
46 | ## Thermal properties |
---|
47 | ## ~~~~~~~~~~~~~~~~~~ |
---|
48 | # Simulate global averaged conditions ? |
---|
49 | global1d = .true. |
---|
50 | # Latitude (deg) [only used if global1d = F] |
---|
51 | latitude = 0.0 |
---|
52 | # Solar Zenith angle (deg) [only used if global1d = T] |
---|
53 | szangle = 45. |
---|
54 | # Force specific heat capacity and molecular mass values |
---|
55 | force_cpp = .false. |
---|
56 | # Specific heat capacity in J K-1 kg-1 [only used if force_cpp = T] |
---|
57 | cpp = 12000. |
---|
58 | # Molecular mass in g mol-1 [only used if force_cpp = T] |
---|
59 | mugaz = 3. |
---|
60 | # Albedo of bare ground |
---|
61 | albedo = 0. |
---|
62 | # Emissivity of bare ground |
---|
63 | emis = 1.0 |
---|
64 | # Soil thermal inertia (SI) |
---|
65 | inertia = 0.001 |
---|
66 | # Initial CO2 ice on the surface (kg.m-2) |
---|
67 | co2ice = 0 |
---|
68 | |
---|
69 | ## Wind profile |
---|
70 | ## ~~~~~~~~~~~~ |
---|
71 | ## zonal eastward component of the geostrophic wind (m/s) |
---|
72 | u = 0. |
---|
73 | # meridional northward component of the geostrophic wind (m/s) |
---|
74 | v = 0. |
---|
75 | |
---|
76 | ## Initial atmospheric temperature profile |
---|
77 | ## ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ |
---|
78 | # |
---|
79 | # Type of initial temperature profile |
---|
80 | # ichoice=1 Constant Temperature: T=tref |
---|
81 | # [Mars] ichoice=2 Savidjari profile (as Seiff but with dT/dz=cte) |
---|
82 | # [Mars] ichoice=3 Lindner (polar profile) |
---|
83 | # [Mars] ichoice=4 inversion |
---|
84 | # [Mars] ichoice=5 Seiff (standard profile, based on Viking entry) |
---|
85 | # ichoice=6 constant T + gaussian perturbation (levels) |
---|
86 | # ichoice=7 constant T + gaussian perturbation (km) |
---|
87 | # ichoice=8 Read in an ascii file "profile" |
---|
88 | ichoice = 1 |
---|
89 | # Reference temperature tref (K) |
---|
90 | tref = 110. ## cold start |
---|
91 | tref = 295. ## warm start |
---|
92 | # Add a perturbation to profile if isin=1 |
---|
93 | isin = 0 |
---|
94 | # peak of gaussian perturbation (for ichoice=6 or 7) |
---|
95 | pic = 26.522 |
---|
96 | # width of the gaussian perturbation (for ichoice=6 or 7) |
---|
97 | largeur = 10 |
---|
98 | # height of the gaussian perturbation (for ichoice=6 or 7) |
---|
99 | hauteur = 30. |
---|
100 | |
---|