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

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

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

# 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.

#---------- 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 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 length of the PEM in Martian years? 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 the glaciers 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 = 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.

#---------- Water reservoirs parameters --------------#
# Amont of water reservoir if the variable is not present in startfi_PEM.nc? Default = 1.e4
# water_reservoir_nom=1e4