## Orbit / general options
## ~~~~~~~~~~~~~~~~~~~~~~~
# Run with or without tracer transport ?
tracer    = .true.
# rate of outputs for diagfi files (in physics steps, default: 24)
diagfi_output_rate = 24
# 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 = .false.
# Check if physics inputs and outputs are ok
#check_physics_inputs = .true.
#check_physics_outputs = .true.

## 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     = 6
# Use blackbody for stellar spectrum ?
stelbbody    = .false.
# Stellar blackbody temperature ?
stelTbb      = 2000.000
# call multilayer correlated-k radiative transfer ?
corrk      = .true.
# folder in which correlated-k data is stored ?
corrkdir    = N2-0.000376CO2-H2Ovar_2026
#corrkdir    = Earth_110-710K
# corrkdir   = N2CO2poor_H2Ovar
# corrkdir   = megaCO2
# corrkdir   = null
# call visible gaseous absorption in radiative transfer ?
callgasvis = .true.
# Include collision-induced absorption (CIA) continuum absorption in radiative transfer ? (make sure it is not already accounted for in your correlated-k tables)
continuum_CIA  = .true.
# Include far wing continuum absorption in radiative transfer ? (make sure it is not already accounted for in your correlated-k tables)
continuum_far_wing  = .true.
# Include Rayleigh scattering in the visible ?
rayleigh   = .true.
# Uniform absorption coefficient in radiative transfer?
graybody     = .false.
# Constant absorption coefficient in visible 
#      (in m^2/kg; only if graybody=true):
#      tau_surf= kappa*P/g
kappa_VI = 1.e-4
# Constant absorption coefficient in IR
#      (in m^2/kg; only if graybody=true):
kappa_IR = 5.e-1 
# Use Newtonian cooling in place of radiative transfer ?
newtonian    = .false.
# Radiative timescale for Newtonian cooling ? [only if newtonian = T]
tau_relax    = 30.00000
# Test physics timescale in 1D ?
testradtimes = .false.
# Characteristic planetary equilibrium (black body) temperature
# This is used only in the aerosol radiative transfer setup. (see aerave.F)
tplanet    = 215.
# 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.
# Remove lower boundary (e.g. for gas giant sims)
noradsurf    = .false. 

## Host Star
## ~~~~~~~~~
stelspec_file = spectrum_Sun_2025.dat
tstellar = 5780.

# 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.

## 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.

#########################################################################
## extra non-standard definitions for Earth
#########################################################################
 
## Tracer and aerosol options 
## ~~~~~~~~~~~~~~~~~~~~~~~~~~
# Number of radiatively active aerosols
naerkind = 1
# Varying H2O cloud fraction?
CLFvarying    = .true.
# H2O cloud fraction?
CLFfixval     = 1.
# number mixing ratio of CO2 ice particles
Nmix_co2      = 1.e5
# basic dust opacity
dusttau       = 0.0
# water cloud pressure level (norm. by psurf)
cloudlvl      = 0.0
# atm mass update due to tracer evaporation/condensation?
mass_redistrib = .true.
# Radiatively active CO2 aerosol?
aeroco2       = .false.
# Fixed CO2 aerosol distribution?
aerofixco2    = .false.
# Radiatively active water aerosol?
aeroh2o       = .true.
# Fixed water aerosol distribution?
aerofixh2o  = .false.
# Radiatively active sulfur aersol?
aeroh2so4     = .false.
# fixed radii for h2o cloud particles?
radfixed=.true.
# number mixing ratio of water particles (for rafixed=.false.)
Nmix_h2o      = 4e6
# number mixing ratio of water ice particles (for rafixed=.false.)
Nmix_h2o_ice      = 2.e4
# radius of H2O water particles (for rafixed=.true.): (CHANGED FROM 10 TO 12 AFTER BENJAMIN)
rad_h2o=12.e-6
# radius of H2O ice particles (for rafixed=.true.):
rad_h2o_ice=35.e-6

## Water options 
## ~~~~~~~~~~~~~
# Model water cycle
water         = .true.
# Model water cloud formation
watercond     = .true.
# Model water precipitation (including coagulation etc.)
waterrain     = .true.
# Moist adjustment
moistadjustment = .true.
# Use simple precipitation scheme?
precip_scheme=4
# multiplicative constant in Boucher 95 precip scheme
Cboucher=0.6
# Include hydrology ?
hydrology     = .true.
# active runoff ?
activerunoff  = .true.
# H2O snow (and ice) albedo ?
albedosnow    = 0.65
# Maximum sea ice thickness ?
maxicethick   = 10.
# Freezing point of seawater (degrees C) ?
Tsaldiff      = 0.0
# Evolve surface water sources ?
sourceevol    = .false.
# compute lightning rate ?
compute_lightning     = .true.

## CO2 options 
## ~~~~~~~~~~~
# call CO2 condensation ?
co2cond       = .false.
# Set initial temperature profile to 1 K above CO2 condensation everywhere?
nearco2cond   = .false.

## Subsurface options
## ~~~~~~~~~~~~~~~~~~
# Number of subsurface layers (For Earth 14 is OK; for Mars 18)
nsoilmx=14
# Thickness of topmost soil layer (m) (For Earth 0.03 is OK; for Mars 0.0002)
lay1_soil=0.03
#Coefficient for soil layer thickness distribution (default is 2)
alpha_soil=2

## Ocean options
## ~~~~~~~~~~~~~
# Model slab-ocean (Main flag for slab ocean)
ok_slab_ocean = .true.
# The following flags can only be set to true if ok_slab_ocean is true
# Ekman transport
slab_ekman = .true.
# Ekman zonal advection
slab_ekman_zonadv = .true.
# Horizontal diffusion (if slab_gm=.true., coef_hdiff=8000, if .false., coef_hdiff=25000, can be changed)
slab_hdiff = .true.
# Slab-ocean timestep (in physics timesteps)
cpl_pas = 1
# Gent-McWilliams Scheme (can only be true if slab_ekman is true)
slab_gm   = .true.
# Slab convective adjustment? 0 - no, 1 - yes
slab_cadj = 1

