datadir=/home1/tbertran/MODELES/KBO/datagcm/
## 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. 
# 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.

## 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.
ch4fix = .true.
vmrch4fix = 0.5
nlte=.false.
# folder in which correlated-k data is stored ?
corrkdir   = pluton_ch4var_COfix
# call visible gaseous absorption in radiative transfer ?
callgasvis = .true.
# Output spectral OLR in 3D?
specOLR    = .false.
# Output global radiative balance in file 'olr.txt' - slow for 1D!!
meanOLR    = .false.

## Star type (inherited from generic model)
## ~~~~~~~~~
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
# 1017.0 W m-2		Sol today x 0.75 = weak early Sun
# 18.462 W m-2		The feeble Gl581
Fat1AU = 1366.0

## Tracer and aerosol options 
## ~~~~~~~~~~~~~~~~~~~~~~~~~~
# Call N2 condensation ?
N2cond       = .true.
# Gravitational sedimentation of tracers ?
sedimentation = .true.
# Compute glacial flow (for fast option only)?
glaflow = .false.
# Compute methane cycle ?
methane = .true.
# Compute CO cycle ?
carbox = .true.
# Compute methane clouds ?
metcloud = .true.
# Compute CO clouds ?
monoxcloud = .true.
# Haze production ? 
haze          = .true.
# source additional haze ?
source_haze= .false.
# 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=.false.
# fractal particle ?
fractal=.false.
# Day fraction for limited cold trap in SP?
dayfrac=0.
# geothermal flux at the bottom layer
fluxgeo = 0.

## Surface properties 
## ~~~~~~~~~~~~~~~~~~~~~~~~~~
constalb=.true.
alb_n2b=0.7
alb_n2a=0.7
alb_ch4_eq=0.5
alb_ch4=0.5
alb_co=0.4
alb_tho=0.2
alb_tho_eq=0.1
emis_n2b=0.8
emis_n2a=0.8
emis_ch4=0.8
emis_co=0.8
emis_tho=1.0
emis_tho_eq=1.0

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

## Other physics options
## ~~~~~~~~~~~~~~~~~~~~~
# call turbulent vertical diffusion ?
calldifv = .true.
# call convective adjustment ?
calladj  = .true.
# call thermal conduction in the soil ?
callsoil = .true.
# call thermal conduction
callconduct=.true.
# call molecular viscosity 3D only
callmolvis=.true.
# call molecular diffusion
callmoldiff=.true.
# flux qui arrive au-dessus de latmosphere
phitop=0.