| [987] | 1 | {\footnotesize |
|---|
| 2 | \begin{verbatim} |
|---|
| 3 | #----------------------------------------------------------------------- |
|---|
| 4 | # Run parameters for the rcm1d.e model |
|---|
| 5 | #----------------------------------------------------------------------- |
|---|
| 6 | |
|---|
| 7 | #### Time integration parameters |
|---|
| 8 | # |
|---|
| 9 | # Initial date (in martian sols ; =0 at Ls=0) |
|---|
| 10 | day0=0 |
|---|
| 11 | # Initial local time (in hours, between 0 and 24) |
|---|
| 12 | time=0 |
|---|
| 13 | # Number of time steps per sol |
|---|
| 14 | day_step=48 |
|---|
| 15 | # Number of sols to run |
|---|
| 16 | ndt =400 |
|---|
| 17 | |
|---|
| 18 | #### Physical parameters |
|---|
| 19 | # |
|---|
| 20 | # Surface pressure (Pa) |
|---|
| 21 | psurf=7000. |
|---|
| 22 | # Gravity (ms^-2) |
|---|
| 23 | g=3.72 |
|---|
| 24 | # Molar mass of atmosphere (g) |
|---|
| 25 | mugaz=43.49 |
|---|
| 26 | # Specific heat capacity of atmosphere? |
|---|
| 27 | cpp=744.5 |
|---|
| 28 | # latitude (in degrees) |
|---|
| 29 | latitude=0.0 |
|---|
| 30 | # orbital distance at perihelion (AU) |
|---|
| 31 | periastr=1.558 |
|---|
| 32 | # orbital distance at aphelion (AU) |
|---|
| 33 | apoastr=1.558 |
|---|
| 34 | # obliquity (degrees) |
|---|
| 35 | obliquit=0.0 |
|---|
| 36 | # Solar zenith angle (degrees) |
|---|
| 37 | szangle=60.0 |
|---|
| 38 | |
|---|
| 39 | # Albedo of bare ground |
|---|
| 40 | albedo=0.2 |
|---|
| 41 | # Emissivity of bare ground |
|---|
| 42 | emis=1.0 |
|---|
| 43 | # Soil thermal inertia (SI) |
|---|
| 44 | inertia=400 |
|---|
| 45 | # zonal eastward component of the geostrophic wind (m/s) |
|---|
| 46 | u=10. |
|---|
| 47 | # meridional northward component of the geostrophic wind (m/s) |
|---|
| 48 | v=0. |
|---|
| 49 | # Initial CO2 ice on the surface (kg.m-2) |
|---|
| 50 | co2ice=0 |
|---|
| 51 | # hybrid vertical coordinate ? (.true. for hybrid and .false. for sigma levels) |
|---|
| 52 | hybrid=.false. |
|---|
| 53 | # autocompute vertical discretisation? (useful for exoplanet runs) |
|---|
| 54 | autozlevs=.false. |
|---|
| 55 | % pressure ceiling |
|---|
| 56 | pceil=40.0 |
|---|
| 57 | |
|---|
| 58 | ###### Initial atmospheric temperature profile |
|---|
| 59 | # |
|---|
| 60 | # Type of initial temperature profile |
|---|
| 61 | # ichoice=1 Constant Temperature: T=tref |
|---|
| 62 | # ichoice=2 Savidjari profile (as Seiff but with dT/dz=cte) |
|---|
| 63 | # ichoice=3 Lindner (polar profile) |
|---|
| 64 | # ichoice=4 inversion |
|---|
| 65 | # ichoice=5 Seiff (standard profile, based on Viking entry) |
|---|
| 66 | # ichoice=6 constant T + gaussian perturbation (levels) |
|---|
| 67 | # ichoice=7 constant T + gaussian perturbation (km) |
|---|
| 68 | # ichoice=8 Read in an ascii file "profile" |
|---|
| 69 | ichoice=5 |
|---|
| 70 | # Reference temperature tref (K) |
|---|
| 71 | tref=200 |
|---|
| 72 | # Add a perturbation to profile if isin=1 |
|---|
| 73 | isin=0 |
|---|
| 74 | # peak of gaussian perturbation (for ichoice=6 or 7) |
|---|
| 75 | pic=26.522 |
|---|
| 76 | # width of the gaussian perturbation (for ichoice=6 or 7) |
|---|
| 77 | largeur=10 |
|---|
| 78 | # height of the gaussian perturbation (for ichoice=6 or 7) |
|---|
| 79 | hauteur=30. |
|---|
| 80 | |
|---|
| 81 | # some definitions for the physics, in file 'callphys.def' |
|---|
| 82 | INCLUDEDEF=callphys.def |
|---|
| 83 | \end{verbatim} |
|---|
| 84 | } |
|---|
