1 | {\footnotesize |
---|
2 | \begin{verbatim} |
---|
3 | ## Orbit / general options |
---|
4 | ## ~~~~~~~~~~~~~~~~~~~~~~~ |
---|
5 | # Run with or without tracer transport ? |
---|
6 | tracer = .true. |
---|
7 | # Diurnal cycle ? if diurnal=false, diurnally averaged solar heating |
---|
8 | diurnal = .true. |
---|
9 | # Seasonal cycle ? if season=false, Ls stays constant, to value set in "start" |
---|
10 | season = .true. |
---|
11 | # Tidally resonant orbit ? must have diurnal=false, correct rotation rate in newstart |
---|
12 | tlocked = .false. |
---|
13 | # Tidal resonance ratio ? ratio T_orbit to T_rotation |
---|
14 | nres = 10 |
---|
15 | # Write some more output on the screen ? |
---|
16 | lwrite = .false. |
---|
17 | # Save statistics in file "stats.nc" ? |
---|
18 | callstats = .true. |
---|
19 | # Test energy conservation of model physics ? |
---|
20 | enertest = .true. |
---|
21 | |
---|
22 | ## Radiative transfer options |
---|
23 | ## ~~~~~~~~~~~~~~~~~~~~~~~~~~ |
---|
24 | # call radiative transfer? |
---|
25 | callrad = .true. |
---|
26 | # the rad. transfer is computed every "iradia" physical timestep |
---|
27 | iradia = 4 |
---|
28 | # call multilayer correlated-k radiative transfer ? |
---|
29 | corrk = .true. |
---|
30 | # folder in which correlated-k data is stored ? |
---|
31 | corrkdir = CO2_H2Ovar |
---|
32 | # call visible gaseous absorption in radiative transfer ? |
---|
33 | callgasvis = .true. |
---|
34 | # Include Rayleigh scattering in the visible ? |
---|
35 | rayleigh = .true. |
---|
36 | # Characteristic planetary equilibrium (black body) temperature |
---|
37 | # This is used only in the aerosol radiative transfer setup. (see aerave.F) |
---|
38 | tplanet = 215. |
---|
39 | # Output spectral OLR in 1D/3D? |
---|
40 | specOLR = .false. |
---|
41 | # Output global radiative balance in file 'rad_bal.out' - slow for 1D!! |
---|
42 | meanOLR = .true. |
---|
43 | # Variable gas species: Radiatively active ? |
---|
44 | varactive = .true. |
---|
45 | # Variable gas species: Fixed vertical distribution ? |
---|
46 | varfixed = .false. |
---|
47 | # Variable gas species: Saturation percentage value at ground ? |
---|
48 | satval = 0.0 |
---|
49 | |
---|
50 | ## Star type |
---|
51 | ## ~~~~~~~~~ |
---|
52 | startype = 1 |
---|
53 | # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ |
---|
54 | # The choices are: |
---|
55 | # |
---|
56 | # startype = 1 Sol (G2V-class main sequence) |
---|
57 | # startype = 2 Ad Leo (M-class, synthetic) |
---|
58 | # startype = 3 GJ644 |
---|
59 | # startype = 4 HD128167 |
---|
60 | # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ |
---|
61 | # Stellar flux at 1 AU. Examples: |
---|
62 | # 1366.0 W m-2 Sol today |
---|
63 | # 1024.5 W m-2 Sol today x 0.75 = weak early Sun |
---|
64 | # 18.462 W m-2 The feeble Gl581 |
---|
65 | # 19.960 W m-2 Gl581 with e=0.38 orbital average |
---|
66 | Fat1AU = 1024.5 |
---|
67 | |
---|
68 | ## Tracer and aerosol options |
---|
69 | ## ~~~~~~~~~~~~~~~~~~~~~~~~~~ |
---|
70 | # Gravitational sedimentation of tracers (just H2O ice for now) ? |
---|
71 | sedimentation = .false. |
---|
72 | |
---|
73 | ## Other physics options |
---|
74 | ## ~~~~~~~~~~~~~~~~~~~~~ |
---|
75 | # call turbulent vertical diffusion ? |
---|
76 | calldifv = .true. |
---|
77 | # call convective adjustment ? |
---|
78 | calladj = .true. |
---|
79 | # call thermal conduction in the soil ? |
---|
80 | callsoil = .true. |
---|
81 | |
---|
82 | ######################################################################### |
---|
83 | ## extra non-standard definitions for Early Mars |
---|
84 | ######################################################################### |
---|
85 | |
---|
86 | ## Tracer and aerosol options |
---|
87 | ## ~~~~~~~~~~~~~~~~~~~~~~~~~~ |
---|
88 | # Fixed aerosol distributions? |
---|
89 | aerofixed = .false. |
---|
90 | # Varying H2O cloud fraction? |
---|
91 | CLFvarying = .false. |
---|
92 | # H2O cloud fraction? |
---|
93 | CLFfixval = 1.0 |
---|
94 | # number mixing ratio of CO2 ice particles |
---|
95 | Nmix_co2 = 100000. |
---|
96 | # number mixing ratio of water ice particles |
---|
97 | Nmix_h2o = 100000. |
---|
98 | |
---|
99 | ## Water options |
---|
100 | ## ~~~~~~~~~~~~~ |
---|
101 | # Model water cycle |
---|
102 | water = .true. |
---|
103 | # Model water cloud formation |
---|
104 | watercond = .true. |
---|
105 | # Model water precipitation (including coagulation etc.) |
---|
106 | waterrain = .true. |
---|
107 | # WATER: Precipitation threshold (simple scheme only) ? |
---|
108 | rainthreshold = 0.0011 |
---|
109 | # Include hydrology ? |
---|
110 | hydrology = .true. |
---|
111 | # H2O snow (and ice) albedo ? |
---|
112 | albedosnow = 0.5 |
---|
113 | # Maximum sea ice thickness ? |
---|
114 | maxicethick = 0.05 |
---|
115 | # Freezing point of seawater (degrees C) ? |
---|
116 | Tsaldiff = 0.0 |
---|
117 | # Evolve surface water sources ? |
---|
118 | sourceevol = .true. |
---|
119 | |
---|
120 | ## CO2 options |
---|
121 | ## ~~~~~~~~~~~ |
---|
122 | # gas is non-ideal CO2 ? |
---|
123 | nonideal = .false. |
---|
124 | # call CO2 condensation ? |
---|
125 | co2cond = .true. |
---|
126 | # Set initial temperature profile to 1 K above CO2 condensation everywhere? |
---|
127 | nearco2cond = .false. |
---|
128 | \end{verbatim} |
---|
129 | } |
---|