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

#---------- Output parameters ----------#
# Do you want to output onto the terminal? Default = .true.
# out2term=.true.

# Do you want to output into a log file? Default = .false.
# out2log=.false.

# Which level of verbosity do you want to? Default = 2
# verbosity_lvl=2

# Value of output rate for diagnostic files? Default = 1
# output_rate=1

# Value of output rate for backup files? Set 0 to disable. Default = 0
# backup_rate=0

#---------- Orbital parameters ----------#
# Do you want to follow an orbital forcing read in file set by "orbit_file_path"? Default = .false.
# evo_orbit=.false.

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

# Orbital forcing ASCII file path. Default = obl_ecc_lsp.asc
# orbit_file_path=obl_ecc_lsp.asc

# Do you want the obliquity to evolve when following orbital forcing file? Default = .true.
# evo_obl=.true.

# Do you want the eccentricity to evolve when following orbital forcing file? Default = .true.
# evo_ecc=.true.

# Do you want the ls perihelie to evolve when following orbital forcing file? Default = .true.
# evo_lsp=.true.

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

#---------- Stopping criteria parameters ----------#
# If evo_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 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.

#---------- Ice management parameters ----------#
# Amount of H2O ice to give to the PCM as a huge perennial reservoir. Default = 9200. kg.m-2 (= 10 m)
# h2oice_huge_resvr=9200.

# Amount of CO2 ice to give to the PCM as a huge perennial reservoir. Default = 16500. kg.m-2 (= 10 m)
# co2ice_huge_resvr=16500.

# Some definitions in file "run_pcm.def" and "callphys.def"
INCLUDEDEF=run_pcm.def
