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