source: trunk/LMDZ.PLUTO/deftank/1d/callphys.def.nohaze @ 3552

Last change on this file since 3552 was 3420, checked in by afalco, 4 months ago

Pluto PCM: callphys.def 1D examples. AF

File size: 4.8 KB
Line 
1## Datagcm path
2datadir=datadir
3## Orbit / general options
4## ~~~~~~~~~~~~~~~~~~~~~~~
5# Run with or without tracer transport ?
6tracer    = .true.
7# Diurnal cycle ?  if diurnal=false, diurnally averaged solar heating
8diurnal   = .false.
9# Seasonal cycle ? if season=false, Ls stays constant, to value set in "start"
10season    = .false.
11noseason_day = 1550
12# Fixed zenith angle in 1D
13global1d = .true.
14# Run 1D or 3D version (false) / fast version volatile transport (true)
15fast    = .false.
16# Write some more output on the screen ?
17lwrite    = .false.
18# Save statistics in file "stats.nc" ?
19callstats = .false.
20# Test energy conservation of model physics ?
21enertest  = .false.
22# run with triton orbit ?
23triton = .false.
24# conservation n2
25conservn2= .false.
26# check if cpp and mugaz from start.nc are consistent with values computed by comp_cpp_mugaz with gases.def
27check_cpp_match = .true.
28check_physics_inputs = .true.
29check_physics_outputs = .true.
30check_ps_min = 0.01
31
32## Radiative transfer options
33## ~~~~~~~~~~~~~~~~~~~~~~~~~~
34# call radiative transfer?
35callrad    = .true.
36# the rad. transfer is computed every "iradia" physical timestep
37iradia     = 1
38# call multilayer correlated-k radiative transfer ?
39corrk      = .true.
40oldplutocorrk=.true.
41ch4fix = .true.
42vmrch4fix = 0.5
43nlte=.true.
44strobel=.false.
45cooling=.true.
46pref=.05
47alpha_top = 1e-11
48deltap = 0.7
49# folder in which correlated-k data is stored ?
50corrkdir   = pluton_ch4var_COfix
51# strictboundcorrk
52strictboundcorrk = .false.
53tplanckmin = 30
54dtplanck=0.1
55# call visible gaseous absorption in radiative transfer ?
56callgasvis = .true.
57# call continuum in radiative transfer ?
58continuum = .false.
59# Include Rayleigh scattering in the visible ?
60rayleigh   = .false.
61# Test physics timescale in 1D ?
62testradtimes = .false.
63# Output spectral OLR in 1D/3D?
64specOLR    = .true.
65# Output global radiative balance in file 'rad_bal.out' - slow for 1D!!
66meanOLR    = .false.
67# Variable gas species: Radiatively active ?
68varactive  = .true.
69# Variable gas species: Fixed vertical distribution ?
70#   (not to be used in time integration mode)
71varfixed   = .false.
72# Use fixed vertical profile, 1 step, no iteration ?
73kastprof     = .false.
74
75# Star type (inherited from generic model)
76startype = 1
77Fat1AU = 1366.0
78
79
80## Other physics options
81## ~~~~~~~~~~~~~~~~~~~~~
82# call turbulent vertical diffusion ?
83calldifv = .true.
84kmixmin=0.001
85vertdiff=.true.
86# use turbdiff instead of vdifc ? (if oldplutovdifc = .false.)
87UseTurbDiff = .true.
88oldplutovdifc=.true.
89# call convective adjustment ?
90calladj  = .true.
91# call thermal conduction in the soil ?
92callsoil = .true.
93lay1_soil = 2e-4
94# geothermal flux at the bottom layer
95fluxgeo = 0.
96# call thermal conduction
97callconduct=.false.
98# call molecular viscosity 3D only
99callmolvis=.false.
100# call molecular diffusion
101callmoldiff=.false.
102# flux qui arrive au-dessus de latmosphere
103phitop=0.
104
105
106## Tracer and aerosol options
107## ~~~~~~~~~~~~~~~~~~~~~~~~~~
108# atm mass update due to tracer evaporation/condensation?
109mass_redistrib = .false.
110# Gravitational sedimentation of tracers ?
111sedimentation = .true.
112oldsedimpluto= .true.
113naerkind = 1
114## Condensation
115## ~~~~~~~~~~~
116generic_condensation = .false.
117# call N2 condensation ?
118n2cond       = .false.
119# N2 condensation subtimestep ?
120nbcond = 20
121# Set initial temperature profile to 1 K above N2 condensation everywhere?
122nearn2cond   = .false.
123# Compute methane cycle ?
124methane = .true.
125condmetsurf = .false.
126# Compute CO cycle ?
127carbox = .true.
128condcosurf=.false.
129# Compute methane clouds ?
130metcloud = .false.
131# Compute CO clouds ?
132monoxcloud = .false.
133# Haze production ?
134haze = .true.
135hazeconservch4    = .true.
136# haze particule mean radius
137rad_haze = 10e-9
138# number mixing ratio of haze particles for ch4 clouds
139Nmix_ch4      = 100000.
140# number mixing ratio of haze particles for co clouds
141Nmix_co      = 100000.
142# Radiatively active haze ?
143aerohaze=.false.
144hazeprop_file=optprop_rannou_r2-200nm_nu003.dat
145# fractal particle ?
146fractal=.false.
147nb_monomer=10
148# Fixed haze radius
149haze_radproffix = .false.
150hazerad_file = rad_haze.txt
151# Fixed haze distribution? # t==  haze_prof (hazemmr_file) else  hazecloud()
152haze_proffix = .false.
153hazemmr_file = hazemmr.txt
154hazedens_file = dens_haze.txt
155
156## Surface properties
157## ~~~~~~~~~~~~~~~~~~~~~~~~~~
158mode_n2=0     ! 0=fixed, 1=change w thickness, 2=special for SP
159deltab=0.    ! if mod_n2=1, difference of albedo depending of thickness
160metlateq=35  ! latitude separating dark/bright methane
161tholateq=25  ! latitude separating eq/pole tholins
162mode_ch4=1
163specalb=.false.  ! Specific treatment for tholins (diff values)
164
165alb_n2b=0.7
166alb_n2a=0.7
167alb_ch4=0.5
168alb_ch4_eq=0.5
169alb_ch4_s=0.5
170alb_co=0.4
171alb_tho=0.2
172alb_tho_eq=0.2
173
174emis_n2b=0.8
175emis_n2a=0.8
176emis_ch4=0.8
177emis_co=0.8
178emis_tho=0.8
179emis_tho_eq=0.8
180
181# Change TI with time
182changeti=.false.
183changetid=.false.
184ITCH4=800.
185ITN2=800
186ITH2O=2000.
187
Note: See TracBrowser for help on using the repository browser.