1 | datadir=/home1/tbertran/MODELES/KBO/datagcm/ |
---|
2 | ## Orbit / general options |
---|
3 | ## ~~~~~~~~~~~~~~~~~~~~~~~ |
---|
4 | # Run with or without tracer transport ? |
---|
5 | tracer = .true. |
---|
6 | # Diurnal cycle ? if diurnal=false, diurnally averaged solar heating |
---|
7 | diurnal = .false. |
---|
8 | # Seasonal cycle ? if season=false, Ls stays constant, to value set in "start" |
---|
9 | season = .true. |
---|
10 | # Run 3D version (false) / fast version (true) |
---|
11 | fast = .true. |
---|
12 | # Write some more output on the screen ? |
---|
13 | lwrite = .false. |
---|
14 | # Save statistics in file "stats.nc" ? |
---|
15 | callstats = .false. |
---|
16 | # Test energy conservation of model physics ? |
---|
17 | enertest = .false. |
---|
18 | # run with triton orbit ? |
---|
19 | triton = .false. |
---|
20 | # specific paleo run ? |
---|
21 | paleo = .false. |
---|
22 | # nb Earth years as paleoclimate step ? (if paleo = true) |
---|
23 | paleoyears = 50000. |
---|
24 | |
---|
25 | ## Radiative transfer options |
---|
26 | ## ~~~~~~~~~~~~~~~~~~~~~~~~~~ |
---|
27 | # call radiative transfer? |
---|
28 | callrad = .true. |
---|
29 | # the rad. transfer is computed every "iradia" physical timestep |
---|
30 | iradia = 1 |
---|
31 | # call multilayer correlated-k radiative transfer ? |
---|
32 | corrk = .false. |
---|
33 | # folder in which correlated-k data is stored ? |
---|
34 | corrkdir = pluton_ch4var_COfix |
---|
35 | # call visible gaseous absorption in radiative transfer ? |
---|
36 | callgasvis = .false. |
---|
37 | # Output spectral OLR in 3D? |
---|
38 | specOLR = .false. |
---|
39 | # Output global radiative balance in file 'olr.txt' - slow for 1D!! |
---|
40 | meanOLR = .false. |
---|
41 | |
---|
42 | ## Star type (inherited from generic model) |
---|
43 | ## ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ |
---|
44 | startype = 1 |
---|
45 | # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ |
---|
46 | # The choices are: |
---|
47 | # |
---|
48 | # startype = 1 Sol (G2V-class main sequence) |
---|
49 | # startype = 2 Ad Leo (M-class, synthetic) |
---|
50 | # startype = 3 GJ644 |
---|
51 | # startype = 4 HD128167 |
---|
52 | # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ |
---|
53 | # Stellar flux at 1 AU. Examples: |
---|
54 | # 1366.0 W m-2 Sol today |
---|
55 | # 1017.0 W m-2 Sol today x 0.75 = weak early Sun |
---|
56 | # 18.462 W m-2 The feeble Gl581 |
---|
57 | Fat1AU = 1366.0 |
---|
58 | |
---|
59 | ## Tracer and aerosol options |
---|
60 | ## ~~~~~~~~~~~~~~~~~~~~~~~~~~ |
---|
61 | # Call N2 condensation ? |
---|
62 | N2cond = .true. |
---|
63 | # N2 condensation subtimestep ? |
---|
64 | nbcond = 20 |
---|
65 | # Gravitational sedimentation of tracers ? |
---|
66 | sedimentation = .false. |
---|
67 | # Compute glacial flow (for fast option only)? |
---|
68 | glaflow = .true. |
---|
69 | # Compute methane cycle ? |
---|
70 | methane = .true. |
---|
71 | # Compute CO cycle ? |
---|
72 | carbox = .false. |
---|
73 | # Haze production ? |
---|
74 | haze = .false. |
---|
75 | # Fast haze production (with option fast=.true.) ? |
---|
76 | fasthaze = .false. |
---|
77 | # number mixing ratio of haze particles for ch4 clouds |
---|
78 | Nmix_ch4 = 100000. |
---|
79 | # number mixing ratio of haze particles for co clouds |
---|
80 | Nmix_co = 100000. |
---|
81 | # atmospheric redistribution : spring constant (s) |
---|
82 | tau_n2 = 1. |
---|
83 | tau_ch4 = 1.E6 |
---|
84 | tau_co = 1. !E5 |
---|
85 | # day fraction for blocking CH4/CO condensation in SP (0 by default) (with methane or carbox=.true.) |
---|
86 | dayfrac=1.0 |
---|
87 | # geothermal flux at the bottom layer |
---|
88 | fluxgeo = 0. |
---|
89 | |
---|
90 | ## Surface properties |
---|
91 | ## ~~~~~~~~~~~~~~~~~~~~~~~~~~ |
---|
92 | # constant albedo with time |
---|
93 | constalb=.false. |
---|
94 | |
---|
95 | alb_n2b=0.7 |
---|
96 | alb_n2a=0.7 |
---|
97 | alb_ch4_eq=0.5 |
---|
98 | alb_ch4=0.5 |
---|
99 | alb_co=0.4 |
---|
100 | alb_tho=0.1 |
---|
101 | alb_tho_eq=0.1 |
---|
102 | emis_n2b=0.8 |
---|
103 | emis_n2a=0.8 |
---|
104 | emis_ch4=0.8 |
---|
105 | emis_co=0.8 |
---|
106 | emis_tho=0.8 |
---|
107 | emis_tho_eq=0.8 |
---|
108 | |
---|
109 | # Change TI with time |
---|
110 | changeti=.false. |
---|
111 | ITCH4=800. |
---|
112 | ITN2=800 |
---|
113 | ITH2O=800. |
---|
114 | |
---|
115 | ## Other physics options |
---|
116 | ## ~~~~~~~~~~~~~~~~~~~~~ |
---|
117 | # call thermal conduction in the soil ? |
---|
118 | callsoil = .true. |
---|
119 | # flux qui arrive au-dessus de latmosphere |
---|
120 | phitop=0. |
---|
121 | |
---|