1 | # |
---|
2 | #----------------------------------------------------------------------- |
---|
3 | # Run parameters for the 1D 'testphys1d.e' model |
---|
4 | #----------------------------------------------------------------------- |
---|
5 | |
---|
6 | #### In/Output parameters |
---|
7 | # |
---|
8 | # Initialization with start files if they are present and writing of restart files |
---|
9 | startfiles_1D=.false. |
---|
10 | # Writing subsurface outputs |
---|
11 | diagsoil=.false. |
---|
12 | |
---|
13 | #### Time integration parameters |
---|
14 | # |
---|
15 | # Initial date (in martian sols ; =0 at Ls=0) |
---|
16 | day0=0 |
---|
17 | # Initial local time (in hours, between 0 and 24) |
---|
18 | time=0. |
---|
19 | # Number of time steps per sol |
---|
20 | day_step=48 |
---|
21 | # Number of sols to run |
---|
22 | ndt = 669 |
---|
23 | |
---|
24 | #### Physical parameters |
---|
25 | # |
---|
26 | # Surface pressure (Pa) |
---|
27 | psurf=610. |
---|
28 | # Reference dust opacity at 610 Pa, in the visible (true tau~tauref*psurf/610) |
---|
29 | tauvis=0.2 |
---|
30 | # latitude (in degrees) |
---|
31 | latitude=0. |
---|
32 | # Albedo of bare ground |
---|
33 | albedo=0.2 |
---|
34 | # Soil thermal inertia (SI) |
---|
35 | inertia=400 |
---|
36 | # zonal eastward component of the geostrophic wind (m/s) |
---|
37 | u=10. |
---|
38 | # meridional northward component of the geostrophic wind (m/s) |
---|
39 | v=0. |
---|
40 | # Initial CO2 ice on the surface (kg.m-2) |
---|
41 | co2ice=0. |
---|
42 | # Water ice cap on ground ? |
---|
43 | watercaptag=.false. |
---|
44 | ## Slopes parameters: |
---|
45 | # (for slope insolation scheme, when callslope=.true. in callphys.def) |
---|
46 | # slope_inclination angle (deg) 0: horizontal, 90: vertical |
---|
47 | slope_inclination=0. |
---|
48 | # slope orientation (deg) |
---|
49 | # 0 == Northward, 90 == Eastward, 180 == Southward, 270 == Westward |
---|
50 | slope_orientation=0. |
---|
51 | |
---|
52 | # hybrid vertical coordinate ? (.true. for hybrid and .false. for sigma levels) |
---|
53 | hybrid=.true. |
---|
54 | |
---|
55 | ###### Initial atmospheric temperature profile |
---|
56 | # |
---|
57 | # Type of initial temperature profile |
---|
58 | # ichoice=1 Constant Temperature: T=tref |
---|
59 | # ichoice=2 Savidjari profile (as Seiff but with dT/dz=cte) |
---|
60 | # ichoice=3 Lindner (polar profile) |
---|
61 | # ichoice=4 inversion |
---|
62 | # ichoice=5 Seiff (standard profile, based on Viking entry) |
---|
63 | # ichoice=6 constant T + gaussian perturbation (levels) |
---|
64 | # ichoice=7 constant T + gaussian perturbation (km) |
---|
65 | # ichoice=8 Read in an ascii file "profile" |
---|
66 | ichoice=5 |
---|
67 | # Reference temperature tref (K) |
---|
68 | tref=200 |
---|
69 | # Add a perturbation to profile if isin=1 |
---|
70 | isin=0 |
---|
71 | # peak of gaussian perturbation (for ichoice=6 or 7) |
---|
72 | pic=26.522 |
---|
73 | # width of the gaussian perturbation (for ichoice=6 or 7) |
---|
74 | largeur=10 |
---|
75 | # height of the gaussian perturbation (for ichoice=6 or 7) |
---|
76 | hauteur=30. |
---|
77 | |
---|
78 | |
---|
79 | # some definitions for the physics, in file 'callphys.def' |
---|
80 | INCLUDEDEF=callphys.def |
---|