## Datagcm path
datadir=datadir
## Orbit / general options
## ~~~~~~~~~~~~~~~~~~~~~~~
# Run with or without tracer transport ?
tracer    = .true.
# Diurnal cycle ?  if diurnal=false, diurnally averaged solar heating
diurnal   = .false.
# Seasonal cycle ? if season=false, Ls stays constant, to value set in "start"
season    = .false.
noseason_day = 1550
# Fixed zenith angle in 1D
global1d = .true.
# Run 1D or 3D version (false) / fast version volatile transport (true)
fast    = .false.
# Write some more output on the screen ?
lwrite    = .false.
# Save statistics in file "stats.nc" ?
callstats = .false.
# Test energy conservation of model physics ?
enertest  = .false.
# run with triton orbit ?
triton = .false.
# conservation n2
conservn2= .false.
# check if cpp and mugaz from start.nc are consistent with values computed by comp_cpp_mugaz with gases.def
check_cpp_match = .true.
check_physics_inputs = .true.
check_physics_outputs = .true.
check_ps_min = 0.01

## Radiative transfer options
## ~~~~~~~~~~~~~~~~~~~~~~~~~~
# call radiative transfer?
callrad    = .true.
# the rad. transfer is computed every "iradia" physical timestep
iradia     = 1
# call multilayer correlated-k radiative transfer ?
corrk      = .true.
oldplutocorrk=.true.
ch4fix = .true.
vmrch4fix = 0.5
nlte=.true.
strobel=.false.
cooling=.true.
pref=.05
alpha_top = 1e-11
deltap = 0.7
# folder in which correlated-k data is stored ?
corrkdir   = pluton_ch4var_COfix
# strictboundcorrk
strictboundcorrk = .false.
tplanckmin = 20
dtplanck=0.1
# call visible gaseous absorption in radiative transfer ?
callgasvis = .true.
# call continuum in radiative transfer ?
continuum = .false.
# Include Rayleigh scattering in the visible ?
rayleigh   = .false.
# Test physics timescale in 1D ?
testradtimes = .false.
# Output spectral OLR in 1D/3D?
specOLR    = .true.
# Output global radiative balance in file 'rad_bal.out' - slow for 1D!!
meanOLR    = .false.
# Variable gas species: Radiatively active ?
varactive  = .true.
# Variable gas species: Fixed vertical distribution ?
#   (not to be used in time integration mode)
varfixed   = .false.
# Use fixed vertical profile, 1 step, no iteration ?
kastprof     = .false.

# Star type (inherited from generic model)
startype = 1
Fat1AU = 1366.0


## Other physics options
## ~~~~~~~~~~~~~~~~~~~~~
# call turbulent vertical diffusion ?
calldifv = .true.
kmixmin=0.001
vertdiff=.true.
# use turbdiff instead of vdifc ? (if oldplutovdifc = .false.)
# UseTurbDiff = .true.
oldplutovdifc=.true.
# call convective adjustment ?
calladj  = .true.
# call thermal conduction in the soil ?
callsoil = .true.
lay1_soil = 2e-4
# geothermal flux at the bottom layer
fluxgeo = 0.
# call thermal conduction
callconduct=.false.
# call molecular viscosity 3D only
callmolvis=.false.
# call molecular diffusion
callmoldiff=.false.
# flux qui arrive au-dessus de latmosphere
phitop=0.


## Tracer and aerosol options
## ~~~~~~~~~~~~~~~~~~~~~~~~~~
# atm mass update due to tracer evaporation/condensation?
mass_redistrib = .false.
# Gravitational sedimentation of tracers ?
sedimentation = .true.
oldsedimpluto= .true.
naerkind = 1
## Condensation
## ~~~~~~~~~~~
generic_condensation = .false.
# call N2 condensation ?
n2cond       = .false.
# N2 condensation subtimestep ?
nbcond = 20
# Set initial temperature profile to 1 K above N2 condensation everywhere?
nearn2cond   = .false.
# Compute methane cycle ?
methane = .true.
condmetsurf = .false.
# Compute CO cycle ?
carbox = .true.
condcosurf=.false.
# Compute methane clouds ?
metcloud = .false.
# Compute CO clouds ?
monoxcloud = .false.
# Haze production ?
haze = .true.
hazeconservch4    = .true.
# number mixing ratio of haze particles for ch4 clouds
Nmix_ch4      = 100000.
# number mixing ratio of haze particles for co clouds
Nmix_co      = 100000.
# Radiatively active haze ?
aerohaze=.true.
hazeprop_file=optprop_rannou_fractal_r010nm_N1_1e4_d2.dat
# fractal particle ?
fractal=.true.
nb_monomer=10
# haze particule mean radius
rad_haze = 10e-9
# Fixed haze radius
haze_radproffix = .false.
hazerad_file = rad_haze.txt
# Fixed haze distribution? # t==  haze_prof (hazemmr_file) else  hazecloud()
haze_proffix = .false.
# hazemmr_file = hazemmr.txt
hazedens_file = dens_haze.txt

## Surface properties
## ~~~~~~~~~~~~~~~~~~~~~~~~~~
mode_n2=0     ! 0=fixed, 1=change w thickness, 2=special for SP
deltab=0.    ! if mod_n2=1, difference of albedo depending of thickness
metlateq=35  ! latitude separating dark/bright methane
tholateq=25  ! latitude separating eq/pole tholins
mode_ch4=1
specalb=.false.  ! Specific treatment for tholins (diff values)

alb_n2b=0.7
alb_n2a=0.7
alb_ch4=0.5
alb_ch4_eq=0.5
alb_ch4_s=0.5
alb_co=0.4
alb_tho=0.2
alb_tho_eq=0.2

emis_n2b=0.8
emis_n2a=0.8
emis_ch4=0.8
emis_co=0.8
emis_tho=0.8
emis_tho_eq=0.8

# Change TI with time
changeti=.false.
changetid=.false.
ITCH4=800.
ITN2=800
ITH2O=2000.