[3550] | 1 | # |
---|
| 2 | #----------------------------------------------------------------------- |
---|
| 3 | # Run parameters for the 1D 'rcm1d.e' model |
---|
| 4 | #----------------------------------------------------------------------- |
---|
| 5 | |
---|
| 6 | #### Time integration parameters |
---|
| 7 | # |
---|
| 8 | # Initial date (in martian sols ; =0 at Ls=0) |
---|
| 9 | day0=1500 |
---|
| 10 | # Initial local time (in hours, between 0 and 24) |
---|
| 11 | time=0 |
---|
| 12 | # Number of time steps per sol |
---|
| 13 | day_step=10 |
---|
| 14 | # Number of sols to run |
---|
| 15 | ndt = 2000 |
---|
| 16 | ecritphy=20 |
---|
| 17 | |
---|
| 18 | cpp_mugaz_mode = 2 |
---|
| 19 | |
---|
| 20 | #### Physical parameters |
---|
| 21 | # |
---|
| 22 | # Day length, in s |
---|
| 23 | daysec = 551854.0625 |
---|
| 24 | # length of a year, in days |
---|
| 25 | year_day = 14178.3037109375 |
---|
| 26 | # Radius of the planet, in m |
---|
| 27 | rad = 1187000 |
---|
| 28 | # Surface pressure (Pa) |
---|
| 29 | psurf=1.1 |
---|
| 30 | # Gravity (ms^-2) |
---|
| 31 | g=0.617188990116119 |
---|
| 32 | # Molar mass of atmosphere (g) |
---|
| 33 | mugaz=28 |
---|
| 34 | # Specific heat capacity of atmosphere (J kg-1 K-1) # ignored if cpp_mugaz_mode = 1 |
---|
| 35 | cpp=1040.0 |
---|
| 36 | # latitude (in degrees) |
---|
| 37 | latitude=0.0 |
---|
| 38 | # orbital distance at perihelion (AU) |
---|
| 39 | periastr=29.65831 |
---|
| 40 | # orbital distance at aphelion (AU) |
---|
| 41 | apoastr=49.3050467 |
---|
| 42 | # date of periastron, in days |
---|
| 43 | peri_day = 87.3600006103516 |
---|
| 44 | # obliquity (degrees) |
---|
| 45 | obliquit=119.591 |
---|
| 46 | |
---|
| 47 | # Albedo of bare ground |
---|
| 48 | albedo=0.5 |
---|
| 49 | # Emissivity of bare ground |
---|
| 50 | emis=0.85 |
---|
| 51 | # Soil thermal inertia (SI) |
---|
| 52 | inertia=20 |
---|
| 53 | # zonal eastward component of the geostrophic wind (m/s) |
---|
| 54 | u=0 |
---|
| 55 | # meridional northward component of the geostrophic wind (m/s) |
---|
| 56 | v=0. |
---|
| 57 | # Initial N2 ice on the surface (kg.m-2) |
---|
| 58 | n2ice=0 |
---|
| 59 | # hybrid vertical coordinate ? (.true. for hybrid and .false. for sigma levels) |
---|
| 60 | hybrid=.false. |
---|
| 61 | |
---|
| 62 | # CH4 VMR (%) |
---|
| 63 | ch4ref=0.5 |
---|
| 64 | # CO VMR (%) |
---|
| 65 | coref=0.046 |
---|
| 66 | # haze ref (kg/kg) |
---|
| 67 | hazeref=0.000000385 |
---|
| 68 | # prec haze ref (kg/kg) |
---|
| 69 | prechazeref=0.0000004 |
---|
| 70 | |
---|
| 71 | # read haze from profhaze |
---|
| 72 | lecthaze = 1 |
---|
| 73 | |
---|
| 74 | ###### Initial atmospheric temperature profile |
---|
| 75 | # |
---|
| 76 | # Type of initial temperature profile |
---|
| 77 | # ichoice=1 Constant Temperature: T=tref |
---|
| 78 | # ichoice=2 Savidjari profile (as Seiff but with dT/dz=cte) |
---|
| 79 | # ichoice=3 Lindner (polar profile) |
---|
| 80 | # ichoice=4 inversion |
---|
| 81 | # ichoice=5 Seiff (standard profile, based on Viking entry) |
---|
| 82 | # ichoice=6 constant T + gaussian perturbation (levels) |
---|
| 83 | # ichoice=7 constant T + gaussian perturbation (km) |
---|
| 84 | # ichoice=8 Read in an ascii file "profile" |
---|
| 85 | # ichoice=1 |
---|
| 86 | ichoice=8 |
---|
| 87 | # Reference temperature tref (K) |
---|
| 88 | tref=38 |
---|
| 89 | # Add a perturbation to profile if isin=1 |
---|
| 90 | isin=0 |
---|
| 91 | # peak of gaussian perturbation (for ichoice=6 or 7) |
---|
| 92 | pic=26.522 |
---|
| 93 | # width of the gaussian perturbation (for ichoice=6 or 7) |
---|
| 94 | largeur=10 |
---|
| 95 | # height of the gaussian perturbation (for ichoice=6 or 7) |
---|
| 96 | hauteur=30. |
---|
| 97 | |
---|
| 98 | |
---|
| 99 | # some definitions for the physics, in file 'callphys.def' |
---|
| 100 | INCLUDEDEF=callphys.def |
---|