#
#----------------------------------------------------------------------- 
#PEM run control parameters:                                          
#---------------------------  

#---------- ORBITAL parameters --------------#

# Do you want to follow an orbital forcing predefined (read in ob_ex_lsp.asc for example)? (default=false)
#  evol_orbit_pem=.false.

# Number of (Martian) year before present to start the pem run if evol_orbit_pem=.true. , default=0
#  year_bp_ini=-507217

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

# Do you want to vary the excenticity when following ob_ex_lsp.asc? Default =.true.
#  var_ex=.true.

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

#---------- Stopping criterion parameters --------------#

# Maximal number of iteration if none of the stopping criterion is reached and if evol_orbit_pem=.false., default=99999999
#  Max_iter_pem=99999999

# Acceptance rate of sublimating water ice surface change, default=0.2
# water_ice_criterion=0.2

# Acceptance rate of sublimating co2 ice surface change, default=0.2
# co2_ice_criterion=0.2

# Acceptance rate of pressure surface change, default=0.15
# ps_criterion = 0.15

# Time step of the PEM in year, default=1
# dt_pem=1

#---------- 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=.true.
# adsorption_pem = .true.

# Do you want to have the glacier to flow along subslope inside a mesh? default=.true.
# co2glaciersflow = .true.

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

#---------- Layering parameters --------------#

# Value of the geothermal flux? default= 0.
# fluxgeo = 0.

#  Depth at which the breccia layer begins? default= 50 m
# depth_breccia   = 10.

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

# Do we compute the ice table at equilibrium?
# icetable_equilibrium = .true.

# Do we compute the ice table with the dynamic method?
# icetable_dynamic = .false.

#---------- Water reservoirs parameters --------------#

# Default amont of water reservoir if the variable is not present in startfi_PEM.nc
# water_reservoir_nom = 1e4