1 | ## Orbit / general options |
---|
2 | ## ~~~~~~~~~~~~~~~~~~~~~~~ |
---|
3 | # Run with or without tracer transport ? |
---|
4 | tracer = .false. |
---|
5 | # Diurnal cycle ? if diurnal=false, diurnally averaged solar heating |
---|
6 | diurnal = .true. |
---|
7 | # Seasonal cycle ? if season=false, Ls stays constant, to value set in "start" |
---|
8 | season = .true. |
---|
9 | # Tidally resonant orbit ? must have diurnal=false, correct rotation rate in newstart |
---|
10 | tlocked = .false. |
---|
11 | # Tidal resonance ratio ? ratio T_orbit to T_rotation |
---|
12 | nres = 1 |
---|
13 | # Write some more output on the screen ? |
---|
14 | lwrite = .false. |
---|
15 | # Save statistics in file "stats.nc" ? |
---|
16 | callstats = .false. |
---|
17 | # Test energy conservation of model physics ? |
---|
18 | enertest = .true. |
---|
19 | # check if cpp matches gases.def (should be false in 1d) |
---|
20 | check_cpp_match = .false. |
---|
21 | # force cpp and mugaz values from gases.def (should be false in 1d) |
---|
22 | force_cpp = .false. |
---|
23 | |
---|
24 | ## Directory where external input files are |
---|
25 | ## ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ |
---|
26 | datadir = ../../datagcm |
---|
27 | |
---|
28 | ## Radiative transfer options |
---|
29 | ## ~~~~~~~~~~~~~~~~~~~~~~~~~~ |
---|
30 | # call radiative transfer? |
---|
31 | callrad = .true. |
---|
32 | # the rad. transfer is computed every "iradia" physical timestep |
---|
33 | iradia = 5 |
---|
34 | # Use blackbody for stellar spectrum ? |
---|
35 | stelbbody = .false. |
---|
36 | # Stellar blackbody temperature ? |
---|
37 | stelTbb = 2000.000 |
---|
38 | # call multilayer correlated-k radiative transfer ? |
---|
39 | corrk = .true. |
---|
40 | # folder in which correlated-k data is stored ? |
---|
41 | corrkdir = 23x23 |
---|
42 | # call visible gaseous absorption in radiative transfer ? |
---|
43 | callgasvis = .true. |
---|
44 | # call continuum in radiative transfer ? |
---|
45 | Continuum = .true. |
---|
46 | # Allow temperatures above kmatrix tmax |
---|
47 | strictboundcorrk = .false. |
---|
48 | # Include Rayleigh scattering in the visible ? |
---|
49 | rayleigh = .true. |
---|
50 | # Uniform absorption coefficient in radiative transfer? |
---|
51 | graybody = .false. |
---|
52 | # Constant absorption coefficient in visible |
---|
53 | # (in m^2/kg; only if graybody=true): |
---|
54 | # tau_surf= kappa*P/g |
---|
55 | kappa_VI = 0. |
---|
56 | # Constant absorption coefficient in IR |
---|
57 | # (in m^2/kg; only if graybody=true): |
---|
58 | kappa_IR = 0. |
---|
59 | # Use Newtonian cooling in place of radiative transfer ? |
---|
60 | newtonian = .false. |
---|
61 | # Radiative timescale for Newtonian cooling ? [only if newtonian = T] |
---|
62 | tau_relax = 30.00000 |
---|
63 | # Test physics timescale in 1D ? |
---|
64 | testradtimes = .false. |
---|
65 | # Characteristic planetary equilibrium (black body) temperature |
---|
66 | # This is used only in the aerosol radiative transfer setup. (see aerave.F) |
---|
67 | tplanet = 94. |
---|
68 | # Output spectral OLR in 1D/3D? |
---|
69 | specOLR = .false. |
---|
70 | # Output global radiative balance in file 'rad_bal.out' - slow for 1D!! |
---|
71 | meanOLR = .false. |
---|
72 | # Remove lower boundary (e.g. for gas giant sims) |
---|
73 | noradsurf = .false. |
---|
74 | |
---|
75 | |
---|
76 | # Compute effective variations of gravity field with altitude in physics ? |
---|
77 | ## Better be true if you run with chemistry |
---|
78 | eff_gz = .true. |
---|
79 | |
---|
80 | ## Chemistry |
---|
81 | ## ~~~~~~~~~ |
---|
82 | # call chemistry ? (must have tracer=true) |
---|
83 | callchim = .false. |
---|
84 | # Activate zonal mean ? |
---|
85 | moyzon_ch = .false. |
---|
86 | # chemistry is computed every "ichim" phys. timestep |
---|
87 | ichim = 200 |
---|
88 | |
---|
89 | ## Microphysics |
---|
90 | ## ~~~~~~~~~~~~ |
---|
91 | # call microphysics ? (must have tracer=true) |
---|
92 | callmufi = .false. |
---|
93 | ## If yes, compute clouds ? |
---|
94 | callclouds = .false. |
---|
95 | ### If yes, number of ices ? (must be compatible with traceur.def AND microphysical model ) |
---|
96 | nices = 3 |
---|
97 | # Desactivate radiative coupling of microphysics ( i.e. use mean opacity profile ) |
---|
98 | uncoupl_optic_haze = .true. |
---|
99 | # ~~~~~ Parameters for microphysics ~~~~~ |
---|
100 | # Path to microphys. config file ? |
---|
101 | config_mufi = ../../datagcm/microphysics/config.ne15.cfg |
---|
102 | # Pressure level of aer. production (Pa) ? |
---|
103 | p_prod = 1.0 |
---|
104 | # Aer. production rate (kg.m-2.s-1) ? |
---|
105 | tx_prod = 3.5e-13 |
---|
106 | # Equivalent radius production (m) ? |
---|
107 | rc_prod = 2.0e-8 |
---|
108 | # Radius of air (nitrogen) molecule (m) ? |
---|
109 | air_rad = 1.75e-10 |
---|
110 | |
---|
111 | |
---|
112 | ## Star type |
---|
113 | ## ~~~~~~~~~ |
---|
114 | startype = 1 |
---|
115 | # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ |
---|
116 | # The choices are: |
---|
117 | # |
---|
118 | # startype = 1 Sol (G2V-class main sequence) |
---|
119 | # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ |
---|
120 | # Stellar flux at 1 AU. Examples: |
---|
121 | # 1366.0 W m-2 Sol today |
---|
122 | # 1024.5 W m-2 Sol today x 0.75 = weak early Sun |
---|
123 | Fat1AU = 1366.0 |
---|
124 | |
---|
125 | ## Other physics options |
---|
126 | ## ~~~~~~~~~~~~~~~~~~~~~ |
---|
127 | # call turbulent vertical diffusion ? |
---|
128 | calldifv = .true. |
---|
129 | # use turbdiff instead of vdifc ? |
---|
130 | UseTurbDiff = .true. |
---|
131 | # call convective adjustment ? |
---|
132 | calladj = .true. |
---|
133 | # call thermal conduction in the soil ? |
---|
134 | callsoil = .true. |
---|