## Datagcm path
datadir=/home/tanguy/SIMUS/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    = .false.
# 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.
# specific paleo run ?
paleo = .true.
# nb Earth years as paleoclimate step ? (if paleo = true)
paleoyears = 1.
# conservation of nitrogen
conservn2 = .false.

## 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.
ch4fix = .false.
vmrch4fix = 0.5
nlte=.false.
cooling=.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
Fat1AU = 1366.0

## Tracer and aerosol options 
## ~~~~~~~~~~~~~~~~~~~~~~~~~~
# Call N2 condensation ?
N2cond       = .false.
# N2 condensation subtimestep ?
nbcond = 20
# Gravitational sedimentation of tracers ?
sedimentation = .false.
# Compute methane cycle ?
methane = .true.
condmetsurf=.true.
# Compute CO cycle ?
carbox = .true.
condcosurf=.true.
# Compute methane clouds ?
metcloud = .false.
# Compute CO clouds ?
monoxcloud = .false.
# Haze production ? 
haze          = .true.
hazeconservch4    = .true.
# Fast haze production (with option fast=.true.) ? 
fasthaze          = .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=.false.
proffix = .true.
# fractal particle ?
fractal=.false.
# geothermal flux at the bottom layer
fluxgeo = 0.
patchflux=0  ! Warm patch of flux in soil.F

## Surface properties 
## ~~~~~~~~~~~~~~~~~~~~~~~~~~
# options to change albedo with time
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=.true.

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.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=.true.
ITCH4=800.
ITN2=800
ITH2O=2000.

## 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=.false.
# call molecular diffusion
callmoldiff=.false.
# flux qui arrive au-dessus de latmosphere
phitop=0.