#_# Orbit / general options
## ~~~~~~~~~~~~~~~~~~~~~~~
# Run with or without tracer transport ?
tracer    = .true.
# Diurnal cycle ?  if diurnal=false, diurnally averaged solar heating
diurnal   = .true.
# Seasonal cycle ? if season=false, Ls stays constant, to value set in "start"
season    = .true.
# Tidally resonant orbit ? must have diurnal=false, correct rotation rate in newstart
tlocked   = .false.
# Tidal resonance ratio ? ratio T_orbit to T_rotation
nres      = 1
# 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.
# check if cpp and mugaz from start.nc are consistent with values computed by comp_cpp_mugaz with gases.def
check_cpp_match = .true.


### run with triton orbit ?
##triton = .false.
### conservation n2
##conservn2=.false.

## Directory where external input files are
## ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
datadir     = datadir

## Radiative transfer options
## ~~~~~~~~~~~~~~~~~~~~~~~~~~
# call radiative transfer?
callrad    = .true.
# the rad. transfer is computed every "iradia" physical timestep
iradia     = 5
# call multilayer correlated-k radiative transfer ?
corrk      = .true.
# folder in which correlated-k data is stored ?
 corrkdir    = pluton_ch4var_COfix
# strictboundcorrk
strictboundcorrk = .true.
# call visible gaseous absorption in radiative transfer ?
callgasvis = .true.
# call continuum in radiative transfer ?
continuum = .false.
# Include Rayleigh scattering in the visible ?
rayleigh   = .true.
# Test physics timescale in 1D ?
testradtimes = .false.
# Output spectral OLR in 1D/3D?
specOLR    = .false.
# 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.
# Variable gas species: Saturation percentage value at ground ?
satval     = .8
# Use fixed vertical profile, 1 step, no iteration ?
kastprof     = .false.

## Star type
## ~~~~~~~~~
startype = 1
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# The choices are:
#
#	startype = 1		Sol        (G2V-class main sequence)
#	startype = 2		Ad Leo     (M-class, synthetic)
#       startype = 3            GJ644
#       startype = 4            HD128167
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# Stellar flux at 1 AU. Examples:
# 1366.0 W m-2		Sol today
# 1024.5 W m-2		Sol today x 0.75 = weak early Sun
# 18.462 W m-2		The feeble Gl581
# 19.960 W m-2		Gl581 with e=0.38 orbital average
Fat1AU = 1366.0

## Tracer and aerosol options
## ~~~~~~~~~~~~~~~~~~~~~~~~~~
# Gravitational sedimentation of tracers (just H2O ice for now) ?
sedimentation = .true.

naerkind = 1
hazeprop_file = optprop_tholins_fractal100nm.dat

## Other physics options
## ~~~~~~~~~~~~~~~~~~~~~
# call turbulent vertical diffusion ?
calldifv = .true.
# use turbdiff instead of vdifc ?
UseTurbDiff = .true.
# call convective adjustment ?
calladj  = .true.
# call thermal conduction in the soil ?
callsoil = .true.
lay1_soil = 1.4e-2

#########################################################################
## extra non-standard definitions for Earth
#########################################################################

## Tracer and aerosol options
## ~~~~~~~~~~~~~~~~~~~~~~~~~~
# atm mass update due to tracer evaporation/condensation?
mass_redistrib = .false.
## Radiatively active CH4 aerosol?
#aeroch4       = .false.

# haze production?
haze = .false.
# Radiatively active haze aerosol?
aerohcaer       = .false.
# Fixed haze distribution?
aerofixhcaer  = .false.
# haze particule mean radius
#rad_hcaer = 10e-9


## Condensation
## ~~~~~~~~~~~
# call N2 condensation ?
n2cond       = .false.
# Set initial temperature profile to 1 K above N2 condensation everywhere?
#nearn2cond   = .false.