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   = .false.
# Seasonal cycle ? if season=false, Ls stays constant, to value set in "start"
season    = .true. 
# Run 3D version (false) / fast version (true)
fast    = .true.
# 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.
# specific paleo run ?
paleo = .false.
# nb Earth years as paleoclimate step ? (if paleo = true)
paleoyears = 50000.

## 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      = .false.
# folder in which correlated-k data is stored ?
corrkdir   = pluton_ch4var_COfix
# call visible gaseous absorption in radiative transfer ?
callgasvis = .false.
# 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.
# N2 condensation subtimestep ?
nbcond = 20
# Gravitational sedimentation of tracers ?
sedimentation = .false.
# Compute glacial flow (for fast option only)?
glaflow = .true.
# Compute methane cycle ?
methane = .true.
# Compute CO cycle ?
carbox = .false.
# Haze production ? 
haze          = .false.
# Fast haze production (with option fast=.true.) ? 
fasthaze          = .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.
# atmospheric redistribution : spring constant (s)
tau_n2 = 1.
tau_ch4 = 1.E6
tau_co = 1. !E5
# day fraction for blocking CH4/CO condensation in SP (0 by default) (with methane or carbox=.true.)
dayfrac=1.0
# geothermal flux at the bottom layer
fluxgeo = 0.

## Surface properties 
## ~~~~~~~~~~~~~~~~~~~~~~~~~~
# constant albedo with time
constalb=.false.

alb_n2b=0.7
alb_n2a=0.7
alb_ch4_eq=0.5
alb_ch4=0.5
alb_co=0.4
alb_tho=0.1
alb_tho_eq=0.1
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.
ITCH4=800.
ITN2=800
ITH2O=800.

## Other physics options
## ~~~~~~~~~~~~~~~~~~~~~
# call thermal conduction in the soil ?
callsoil = .true.
# flux qui arrive au-dessus de latmosphere
phitop=0.