#------------------------------------#
# Run control parameters for the PEM #
#------------------------------------#

#---------- Output parameters ----------#
# Output rate? Default = 1
output_rate=1

#---------- Orbital parameters ----------#
# Do you want to follow an orbital forcing read in "obl_ecc_lsp.asc"? Default = .false.
# evol_orbit=.false.

# If evol_orbit=.true., number of Earth years before present to start the PEM run. Default = 0
# pem_ini_earth_date=0

# Do you want the obliquity to evolve when following "obl_ecc_lsp.asc"? Default = .true.
# evol_obl=.true.

# Do you want the eccentricity to evolve when following "obl_ecc_lsp.asc"? Default = .true.
# evol_ecc=.true.

# Do you want the ls perihelie to evolve when following "obl_ecc_lsp.asc"? Default = .true.
# evol_lsp=.true.

# Time step length of the PEM in Martian years? Default = 1.
# dt=1

#---------- Stopping criteria parameters ----------#
# If evol_orbit=.false., maximal number of iterations if no stopping criterion is reached. Default=100000000.
# nmax_yr_run=100000000

# Acceptance rate of sublimating H2O ice surface change. Default = 0.2
# h2oice_crit=0.2

# Acceptance rate of sublimating CO2 ice surface change. Default = 0.2
# co2ice_crit=0.2

# Acceptance rate of pressure surface change. Default = 0.15
# ps_crit = 0.15

# Acceptance of change for obliquity. Default = 1.
# max_change_obl=1.

# Acceptance of change for eccentricity? Default = 5.e-3
# max_change_ecc=5.e-3

# Acceptance of change for Lsp? Default = 20.
# max_change_lsp=20.

#---------- Subsurface parameters ----------#
# Do you want to run with subsurface physical processes in the PEM? Default = .true.
# do_soil=.true.

# Do you want to run with adsorption/desorption in the PEM? Default = .false.
# do_sorption=.false.

# Do you want to modify the soil thermal properties with the pressure? Default = .false.
# reg_thprop_dependp=.false.

# Value of the geothermal flux. Default = 0.
# flux_geo=0.

# Depth at which the breccia layer begins. Default = 10 m
# depth_breccia=10.

# Depth at which the bedrock layer begins. Default = 1000 m
# depth_bedrock=1000.

# Is the ice table computed at equilibrium? Default = .true.
# icetable_equilibrium=.true.

# Is the ice table computed with the dynamic method? Default = .false.
# icetable_dynamic=.false.

#---------- Ice management parameters ----------#
# Amount of H2O ice to initialize the huge reservoir if the variable is not present in "startevol.nc". Default = 9200. kg.m-2 (= 10 m)
# h2oice_huge_ini=9200.

# Threshold to consider the amount of H2O ice as an infinite reservoir. Default = 460. kg.m-2 (= 0.5 m)
# threshold_h2oice_cap=460.

# Do you want the H2O ice to flow along subslope inside a cell? Default = .true.
# h2oice_flow=.true.

# Do you want the CO2 ice to flow along subslope inside a cell? Default = .true.
# co2ice_flow=.true.

#---------- Layering parameters ----------#
# Do you want to run with the layering algorithm? Default = .false.
# do_layering=.false.
 
# Value of the dust tendency. Default = 5.78e-2 kg.m-2.y-1
# d_dust=5.78e-2
 
# Do you want to impose a dust-to-ice ratio instead of a dust tendency? Default = .false.
# impose_dust_ratio=.false.
 
# If impose_dust_ratio=.true., value of the dust-to-ice ratio. Default = 0.1
# dust2ice_ratio=0.1

# Some definitions for the physics in file "run_pcm.def"
INCLUDEDEF=run_pcm.def
