##General options
##~~~~~~~~~~~~~~~
#Run with or without tracer transport ?
tracer=.true.

#Diurnal cycle ?  if diurnal=False, diurnal averaged solar heating
diurnal=.true.

#Seasonal cycle ? if season=False, Ls stays constant, to value set in "start"
season = .true.   

#write some more output on the screen ?
lwrite = .false.

#Save statistics in file "stats.nc" ?
callstats =.false.

#Save EOF profiles in file "profiles" for Climate Database?
calleofdump = .false.

#Directory where external input files are:
datadir=datadir

## use (read/write) a startfi.nc file? (default=.true.)
startphy_file = .false.
surfemis      =0.5
surfalbedo    =0.5

## Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F)
## ~~~~~~~~~~~~~
#  =1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario,
#  =4 Mars Year 24 from TES assimilation (old version of MY24; dust_tes.nc file)
#  =6 "cold" (low dust) scenario ; =7 "warm" (high dust) scenario
#  =8 "clim" scenario
#  =24 Mars Year 24 from TES assimilation (ie: MCD reference case)
#  =25 Mars Year 25 from TES assimilation (ie: a year with a global dust storm)
#  =26 Mars Year 26 from TES assimilation
iaervar = 8
# Dust opacity at 700 Pa (when constant, i.e. for the iaervar=1 case)
tauvis=0.5
# Dust vertical distribution: 
# (=0: old distrib. (Pollack90), =1: top set by "topdustref",
#  =2: Viking scenario; =3 MGS scenario)
iddist  = 3
# Dust top altitude (km). (Matters only if iddist=1)
topdustref = 55.

## Physical Parameterizations :
## ~~~~~~~~~~~~~~~~~~~~~~~~~~
# call radiative transfer ?
callrad   = .false.
# call NLTE radiative schemes ?   matters only if callrad=T
callnlte = .false.
# NLTE 15um scheme to use.                                                                           
# 0-> Old scheme, static oxygen                                                                      
# 1-> Old scheme, dynamic oxygen                                                                     
# 2-> New scheme                                                                                     
nltemodel = 0 
# call CO2 NIR absorption ?   matters only if callrad=T
callnirco2 = .false.
# NIR NLTE correction ? matters only if callnirco2=T                                                 
nircorr=0
# call turbulent vertical diffusion ?
calldifv  = .false.
# call convective adjustment ?
calladj  = .false.
# Thermals
calltherm = .false.
callrichsl = .false.
# call CO2 condensation ?
callcond =.false.
# call thermal conduction in the soil ?
callsoil = .false.
# call Lott's gravity wave/subgrid topography scheme ?
calllott  = .false.
calllott_nonoro  = .false.
# Impose polar cap surface albedos as observed by TES?
TESicealbedo = .true.
## Coefficient for Northern cap albedoes
TESice_Ncoef=1.6
## Coefficient for Southern cap albedoes
TESice_Scoef=1.6


## Radiative transfer options :
## ~~~~~~~~~~~~~~~~~~~~~~~~~~
# the rad.transfer is computed every "iradia" physical timestep
iradia    = 1
# Output of the exchange coefficient mattrix ? for diagnostic only
callg2d  = .false.
# Rayleigh scattering : (should be .false. for now)
rayleigh  = .false.

## Tracer (dust water, ice and/or chemical species) options (used if tracer=T):
## ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# DUST: Transported dust ? (if >0, use 'dustbin' dust bins)
dustbin    = 0
# DUST: Radiatively active dust ? (matters if dustbin>0)
active  = .false.
# DUST: use mass and number mixing ratios to predict dust size ?
doubleq   = .false.
# DUST: lifted by GCM surface winds ?
lifting = .false.
# DUST: lifted by dust devils ?
callddevil = .false.
# DUST: Scavenging by H2O snowfall ?
scavenging = .false.
# DUST: Scavenging by instantaneous CO2 snowfall ?
scavco2cond = .false.

# DUST: freedust ?
freedust=.false.
# DUST: dust injection scheme (if >0, use 'dustinjection')
dustinjection=0
# DUST: dust injection scheme coefficient
coeff_injection=0.25
# DUST: Rocket Dust Storm scheme
rdstorm=.false.
# DUST: slpwind scheme
slpwind = .false.

# DUST/WATERICE: Gravitationnal sedimentation ?
sedimentation = .false.


# WATERICE: Radiatively active transported atmospheric water ice ?
activice   = .false.
# WATER: Compute water cycle
water = .false.
# WATER: Microphysical scheme for water-ice clouds?
microphys = .false.
# WATER: Microphysical time step for water-ice clouds?  (30 by default)
#imicro = 100
# WATER: parameter contact
mteta = 0.95
# WATER: current permanent caps at both poles. True IS RECOMMENDED
#        (with .true., North cap is a source of water and South pole
#         is a cold trap)
caps  = .true.
# WATER: Water ice albedo ?
albedo_h2o_ice = 0.35
# WATER: Water ice thermal inertia
inert_h2o_ice = 800
# WATER: Frost thickness threshold for albedo
frost_albedo_threshold = 0.005
nuice_sed=0.1
# PHOTOCHEMISTRY: include chemical species 
photochem  = .false.
#WATER: subgrid clouds
CLFvarying=.false.
#WATER: subgrid clouds
CLFfixval=1
#WATER: naerkind
naerkind=0

## Thermospheric options (relevant if tracer=T) :
##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# call thermosphere ?
callthermos = .false.
# WATER: included without cycle (only if water=.false.)
thermoswater = .false.
# call thermal conduction ? (only if callthermos=.true.)
callconduct = .false.
# call EUV heating ? (only if callthermos=.true.)
calleuv=.false.
# call molecular viscosity ? (only if callthermos=.true.)
callmolvis = .false.
# call molecular diffusion ? (only if callthermos=.true.)
callmoldiff = .false.
# call thermospheric photochemistry ? (only if callthermos=.true.)
thermochem = .false.
# date for solar flux calculation: (1985 < date < 2002)
## (Solar min=1996.4 ave=1993.4 max=1990.6)
solarcondate = 1993.4