## 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.
noseason_day = 0
# 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

# Frequency (in physical steps) of diagfi file writing ?
diagfi_output_rate = 32
# Save statistics in file "stats.nc" ?
callstats = .false.
# Test energy conservation of model physics ?
enertest  = .false.

# Check the validity of input fields coming from the dynamics ?
check_physics_inputs  = .false.
# Check the validity of updated fields at the end of the physics step ?
check_physics_outputs = .false.
# Check minimum reasonable surface pressure (Pa) ?
check_ps_min          = 1e-8

# Conservation tracers diagnostics
conservn2  = .false.
conservch4 = .false.

# run with triton orbit ?
triton = .false.

## Directory where external input files are
## ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
datadir = /scratchu/bdebatz/Pluto/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.
oldplutocorrk    = .false.
# Compute radiative methane?
ch4fix           = .true.
vmrch4fix        = 0.5
vmrch4_proffix   = .false.
nlte             = .false.
strobel          = .false.
cooling          = .false.
pref             = 5.0e-2
alpha_top        = 1e-11
deltap           = 0.7
# folder in which correlated-k data is stored ?
corrkdir         = pluton_ch4var_COfix_17x27
# strictboundcorrk
strictboundcorrk = .false.
strictboundcia   = .false.
tplanckmin       = 5
# call visible gaseous absorption in radiative transfer ?
callgasvis       = .true.
# call continuum in radiative transfer ?
continuum        = .true.
# 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.
# Use fixed vertical profile, 1 step, no iteration ?
kastprof         = .false.

## Star type (startype = 1 <=> Sol)
## ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
startype = 1
Fat1AU = 1366.0


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


## Tracer options
## ~~~~~~~~~~~~~~
# atm mass update due to tracer evaporation/condensation?
mass_redistrib = .false.

### Volatile cycles
### ~~~~~~~~~~~~~~~
# call N2 condensation ?
n2cond      = .true.
# N2 condensation subtimestep ?
nbsub       = 20
# No formation of N2 frost ?
no_n2frost  = .true.
# Compute methane cycle ?
methane     = .true.
condmetsurf = .true.
# Compute CO cycle ?
carbox      = .true.
condcosurf  = .false.
# Compute methane clouds ?
metcloud    = .false.
# Compute CO clouds ?
monoxcloud  = .false.

## Microphysics options
## ~~~~~~~~~~~~~~~~~~~~
# call microphysics moment model?
callmufi    = .true.
# Path to microphysical config file.
config_mufi = /scratchu/bdebatz/Pluto/datadir/microphysics/config.cfg

### Radiatively active aerosols
### ~~~~~~~~~~~~~~~~~~~~~~~~~~~
# Number of radiatively active aerosols (if callmufi = true, must be = 2)?
naerkind      = 2
# Radiatively active haze aerosol?
optichaze     = .true.
##aersprop_file = optprop_rannou_r2-200nm_nu003.dat
aersprop_file = optprop_lora_pluto_r4-1000nm_nu003.dat
aerfprop_file = optprop_rannou_fractal_r020nm_N1_1e4_d2.dat
##aerfprop_file = optprop_lora_fractal_PH_r010nm_N1_5e5_d2.dat

### Aerosol options (if callmufi = .true.)
### ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# Equivalent radius production (m).
haze_rc_prod        = 8.0e-9
# Monomer radius (m).
haze_rm             = 20e-9
# Aerosol's fractal dimension.
haze_df             = 2.0
# Aerosol density (kg.m-3).
haze_rho            = 800.
# Radius of air molecule (m) ? (Pluto : N2)
air_rad             = 1.75e-10
# Use production from CH4 photolysis (if not, use production rate)?
call_haze_prod_pCH4 = .true.
hazeconservch4      = .true.
# (call_haze_prod_pCH4 == false): 
# Pressure level of aerosol production (Pa).
haze_p_prod         = 1.0e-2
# Aerosol production rate (kg.m-2.s-1).
haze_tx_prod        = 9.8e-14

### Old aerosol options (if callmufi = .false.)
### ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# haze production?
haze            = .false.
# Gravitational sedimentation of tracers ?
sedimentation   = .false.
# Fixed haze distribution? (false == hazecloud, true == haze_prof)
haze_proffix    = .false.
hazemmr_file    = hazemmr.txt
hazedens_file   = dens_haze.txt
# Haze particule mean radius
rad_haze        = 10e-9
haze_radproffix = .false.
hazerad_file    = rad_haze.txt
fractal         = .false.
nb_monomer      = 10
hazeprop_file   = optprop_rannou_fractal_r010nm_N1_1e4_d2.dat
# Number mixing ratio of haze particles for ch4 clouds
Nmix_ch4        = 100000.
# Number mixing ratio of haze particles for co clouds
Nmix_co         = 100000.


## Surface properties
## ~~~~~~~~~~~~~~~~~~~~~~~~~~
# options to change albedo with time
#!  1) ALBEDOS and EMISSIVITY
#!       A. N2
#!                 CASE (0) ! fixed albedo
#!                 CASE (1) ! Albedo decreases with thickness
#!                 CASE (2) ! Special Sputnik differences of albedo
#!                 CASE (3) ! Albedo increases (delta neg) or decreases (delta pos) with sublimationi rates
#!                 CASE (4) ! Albedo Difference in N/S (e.g. used for Triton)
#!                 CASE (5) ! Special Sputnik differences of albedo in small (1 pixel) patches (e.g. simulating dark patches / plumes)
#!            --> EMISSIVITY N2: based on the alpha/beta transition
#!       B. CO
#!       C. CH4
#!                 CASE (0) ! 2 albedos, one for the tropics, one for the poles
#!                 CASE (1) ! 3 albedos, one for the tropics, 2 for the poles (north and south)
#!                 CASE (2) ! 2 albedos + albedo feedback
#!                          SELECT CASE (feedback_met)
#!                            CASE (0) ! Default (linear from alb_ch4_eq)
#!                            CASE (1) ! Hyperbolic tangent old
#!                            CASE (2) ! hyperbolic tangent old
#!                            CASE (3) ! hyperbolic tangent equation with parameters
#!                 CASE (3) ! Eq, poles N, pole S + depending on Ls
#!       D. Tholins
#!                 CASE (0) ! Default, 2 albedos, one for the tropics, one for the poles
#!                 CASE (1) ! Special mode one region with a different albedo
#!       E. Tholins read from file
#                 ! specalb

# N2
mode_n2     = 0
thres_n2ice = 0.001
deltab      = 0. ! if mod_n2 = 1, difference of albedo depending of thickness

# CH4
mode_ch4     = 1
thres_ch4ice = 0.01
metlateq     = 35 ! latitude separating dark/bright methane

# CO
thres_coice = 1.

# Tholins
mode_tholins = 0
tholateq     = 25      ! latitude separating eq/pole tholins
specalb      = .false. ! Specific treatment for tholins (diff values)

alb_n2b    = 0.74
alb_n2a    = 0.74
alb_ch4    = 0.7
alb_ch4_eq = 0.55
alb_ch4_s  = 0.6
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    = 1.
emis_tho_eq = 1.

# Change TI with time
changeti=.false.
ITCH4  = 800.
ITN2   = 800.
ITH2O  = 800.
changetid=.false.
ITCH4d = 20.
ITN2d  = 20.
ITH2Od = 20.