Difference between revisions of "Advanced Topics Mars PEM"

From Planets
Jump to: navigation, search
(Created page with "== Choices of parameters (run_PEM.def) == 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 p...")
 
(→‎Running the PEM in 1D)
 
(4 intermediate revisions by 2 users not shown)
Line 1: Line 1:
−
== Choices of parameters (run_PEM.def) ==
+
== Running the PEM with the launching script ==
  
−
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:
+
The script ''launchPEM.sh'' can take an argument:
 +
      # If there is no argument, then the script initiates a PEM simulation from scratch.
 +
      # If the argument is 're', then the script relaunches an existing PEM simulation. It will ask for parameters to know the starting point that you want to.  
  
 +
== Running the PEM in 1D ==
  
 +
One can run the PEM in a single column configuration, i.e. as a 1D model. Compilation is done using the [[The makelmdz fcm GCM Compilation Script|makelmdz_fcm]] script except that only a number of vertical layers needs be specified. For example, to compile a version of 32 vertical layers, one would run the following command:
 +
<pre>
 +
./makelmdz_fcm -arch somearch -p mars -d 32 -j 8 pem
 +
</pre>
  
−
'''General parameters for the run '''
+
To work with the 1D PEM, you must use the related PCM in 1D. For example, for the Mars PCM, see [[Mars 1D testphys1d program]].
  
−
* evol_orbit_pem: [Boolean] Do you want to follow an orbital forcing predefined (read in ob_ex_lsp.asc for example)? (default=false)
+
[[Category:Planetary-Evolution-Model]]
−
** 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
 
−
** 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.
 
−
** var_lsp: [Boolean] Do you want to vary the ls perihelie when following ob_ex_lsp.asc? Default =.true.
 
−
* 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
 
−
* ice_criterion: [Real] Acceptance rate of sublimating ice surface change, default=0.2
 
−
* ps_criterion: [Real] Acceptance rate of pressure surface change, default=0.15
 
−
 
 
−
'''Subsurface'''
 
−
 
 
−
* 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.
 
−
* 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.
 
−
* depth_breccia: [Real]. Depth at which the breccia layer begins? default= 50 m
 
−
* depth_bedrock: [Real]. Depth at which the bedrock layer begins? default= 1000 m
 
−
 
 
−
''' Reservoirs '''
 
−
* water_reservoir_nom: [Real]. Value of the available water in the subsurface. default =1e4 in kg/m^2
 

Latest revision as of 15:49, 23 October 2025

Running the PEM with the launching script

The script launchPEM.sh can take an argument: 

     # If there is no argument, then the script initiates a PEM simulation from scratch.
     # If the argument is 're', then the script relaunches an existing PEM simulation. It will ask for parameters to know the starting point that you want to.

Running the PEM in 1D

One can run the PEM in a single column configuration, i.e. as a 1D model. Compilation is done using the makelmdz_fcm script except that only a number of vertical layers needs be specified. For example, to compile a version of 32 vertical layers, one would run the following command: 

./makelmdz_fcm -arch somearch -p mars -d 32 -j 8 pem

To work with the 1D PEM, you must use the related PCM in 1D. For example, for the Mars PCM, see Mars 1D testphys1d program.

Retrieved from "http://lmdz-forge.lmd.jussieu.fr/mediawiki/Planets/index.php?title=Advanced_Topics_Mars_PEM&oldid=2908"