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

#---------- Planet type ----------#
planet_type=mars

#---------- Output parameters ----------#
# Frequency of outputs for the PEM? Default = 1
ecritpem=1

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

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

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

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

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

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

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

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

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

# Acceptance rate of pressure surface change. Default = 0.15
# ps_criterion = 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.
# soil_pem=.true.

# Do you want to run with adsoprtion in the PEM? Default = .false.
# adsorption_pem=.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.
# fluxgeo=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 "startfi_PEM.nc". Default = 9200. kg.m-2 (= 10 m)
# ini_huge_h2oice=9200.

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

# Do you want H2O frost to transform into perennial H2O ice? Default = .false.
# metam_h2oice=.false.

# Threshold to consider frost is becoming perennial H2O ice. Default = 460. kg.m-2 (= 0.5 m)
# metam_h2oice_threshold=460.

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

# Do you want CO2 frost to transform into perennial CO2 ice? Default = .false.
# metam_co2ice=.false.

# Threshold to consider frost is becoming perennial CO2 ice. Default = 16500. kg.m-2 (= 10 m)
# metam_co2ice_threshold=16500.

# 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.
# layering=.false.
 
# Value of the dust tendency. Default = 1.e-3 kg.m-2.y-1
# d_dust=1.e-3
 
# 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.01
# dust2ice_ratio=0.01

# Some definitions for the physics, in file 'callphys.def'
INCLUDEDEF=callphys.def