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

# 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 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=1.e4

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

# 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=5.e-2

# 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=16.e3

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

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