PEM (Planetary Evolution Model) - Revision history

Jbclement at 12:12, 21 July 2025

2025-07-21T12:12:14Z

2023-03-09T10:49:20Z

‎Choices of parameters (run_PEM.def)

2023-03-09T10:14:57Z

2023-02-10T15:44:15Z

‎Choices of parameters (run_PEM.def)

2023-02-06T11:16:45Z

‎Choices of parameters (run_PEM.def)

2023-02-06T11:11:07Z

‎Choices of parameters (run_PEM.def)

2023-02-02T11:36:11Z

‎Choices of parameters (run_PEM.def)

2023-02-02T11:17:06Z

‎Physical processes modelled

2023-02-02T11:07:42Z

‎Overview of a PEM run

2023-02-02T11:07:31Z

‎Overview of a PEM run

@@ Line 162: / Line 162: @@
 ''' Reservoirs '''
 * water_reservoir_nom: [Real]. Value of the available water in the subsurface. default =1e4 in kg/m^2
 [[Category:Planetary-Evolution-Model]]

@@ Line 141: / Line 141: @@
 * evol_orbit_pem: [Boolean] Do you want to follow an orbital forcing predefined (read in ob_ex_lsp.asc for example)? (default=false)
 ** If evol_orbit_pem is set to true, year_bp_ini needs to be specified
-* year_bp_ini: [Integer] Number of (Earth) year before present to start the pem run if evol_orbit_pem=.true. , default=0
+** year_bp_ini: [Integer] Number of (Earth) year before present to start the pem run if evol_orbit_pem=.true. , default=0
 ** var_obl: [Boolean] Do you want to vary the obliquity when following ob_ex_lsp.asc? Default =.true.
 ** var_ex: [Boolean] Do you want to vary the excenticity when following ob_ex_lsp.asc? Default =.true.

← Older revision		Revision as of 10:14, 9 March 2023
Line 177:		Line 177:

	An example of script to run the PEM in a cluster can be found in the following directory : trunk/???		An example of script to run the PEM in a cluster can be found in the following directory : trunk/???
		+
		+	[[Category:Planetary-Evolution-Model]]

@@ Line 140: / Line 140: @@
 * evol_orbit_pem: [Boolean] Do you want to follow an orbital forcing predefined (read in ob_ex_lsp.asc for example)? (default=false)
 * year_bp_ini: [Integer] Number of (Earth) year before present to start the pem run if evol_orbit_pem=.true. , default=0
 * Max_iter_pem: [Integer] Maximal number of iteration if none of the stopping criterion is reached and if evol_orbit_pem=.false., default=99999999
 * dt_pem: [Integer] Time step of the PEM in year, default=1

@@ Line 157: / Line 157: @@
 ''' Reservoirs '''
-* water_reservoir_nom: [Real]. Value of the avaible water in the subsurface. default =1e4 in kg/m^2
+* water_reservoir_nom: [Real]. Value of the available water in the subsurface. default =1e4 in kg/m^2
 == How to compile the model ==

@@ Line 150: / Line 150: @@
 * soil_pem: [Boolean]. Do you want to run with subsurface physical processes in the PEM? default=.true.
 * adsorption_pem: [Boolean]. Do you want to run with adsoprtion in the PEM? default=.true.
-* fluxgeo: [Real]. Value of the geothermal flux? default= 0.
+* reg_thprop_dependp: [Boolean]. Do you want to modify the soil thermal properties with the pressure ? default=.false.
 ''Layering:''
 * fluxgeo: [Real]. Value of the geothermal flux? default= 0.

@@ Line 135: / Line 135: @@
 Parameters of the PEM run can be specified in run_PEM.def. This file is included at the end of the run.def. For now, these are the possible parameters:
-* evol_orbit_pem: Boolean. Do you want to follow an orbital forcing predefined (read in ob_ex_lsp.asc for example)? (default=false)
-* year_bp_ini: Integer. Number of (Earth) year before present to start the pem run if evol_orbit_pem=.true. , default=0
-* Max_iter_pem: Integer. Maximal number of iteration if none of the stopping criterion is reached and if evol_orbit_pem=.false., default=99999999
+'''General parameters for the run '''
-* dt_pem: Integer. Time step of the PEM in year, default=1
-* alpha_criterion: Real. Acceptance rate of sublimating ice surface change, default=0.2
+* evol_orbit_pem: [Boolean] Do you want to follow an orbital forcing predefined (read in ob_ex_lsp.asc for example)? (default=false)
-* soil_pem: Boolean. Do you want to run with subsurface physical processes in the PEM? default=.true.
+* year_bp_ini: [Integer] Number of (Earth) year before present to start the pem run if evol_orbit_pem=.true. , default=0
 == How to compile the model ==

@@ Line 53: / Line 53: @@
 == Physical processes modelled ==
-The PEM is still in developement, so far, physical processes modelled are quite specific to Mars at low obliquity.
+The PEM is still in development, so far, the physical processes modeled are quite specific to Mars at low obliquity.
-✅ Ice accumulation (Improvement possibles like changes of topography, watercaptag etc)
+'''Current processes that have been validated:'''
-❌ Lag deposit, stratification (See Simon Neviere, Futur intern of Lucas)
+✅ Ice accumulation
-∼ Glacier flow (CO2) , ✅ subgrid-scale slope phenomena
+✅ CO2 / H_2O adsorption
-∼ Subsurface Ground Ice  (In validation)
+✅ Soil temperature from the surface to ten's of kms.
 ❌ Hydrology: lake, river, ocean, etc. (Alexandre Gauvain)
 == PEM inputs and outputs ==

@@ Line 42: / Line 42: @@
 In either mode, a run works like this:
-* Run 2 years of GCM: From this, we extract data such as daily average pressure, daily averaged vmr of atmospheric components, daily minimum of ice etc... You thuss need to use '''XIOS'''
+* Run 2 years of GCM: From this, we extract data such as daily average pressure, daily averaged vmr of atmospheric components, daily minimum of ice etc... You thus need to use '''XIOS'''
 * The GCM restart file are renamed to be read by the PEM. They will be modified by the PEM to be able to start a new GCM simulation.
 * If we are in a M2 configuration mode, a beginning year needs to be specified in the run_PEM.def, and the file ob_ex_lsp.asc needs to be present.