Changeset 3989 for trunk/LMDZ.COMMON/libf/evolution/deftank

Timestamp:

Dec 11, 2025, 12:56:05 PM (8 days ago)

Author:

jbclement

Message:

PEM:
Massive structural refactor of the PEM codebase for improved readability, consistency and maintainability. The goal is to modernize, standardize and consolidate the code while removing legacy complexity. In detail, this change:

• Performs large-scale cleanup removing unused code, obsolete routines, duplicated functionality and deprecated initialization logic;
• Removes "*_mod" wrappers;
• Replaces mixed naming conventions with clearer variable names, domain-based file/module names and purpose-based routine names;
• Adds native reading/writing capabilities to the PEM removing the cumbersome dependency on Mars PCM subroutines;
• Centralizes soil/ice/orbital/PEM configuration logic into dedicated modules;
• Simplifies routines and updates calls/interfaces to match the new structure.

This refactor significantly clarifies the codebase and provides a cleaner foundation for forthcoming developments.
JBC

File:

• trunk/LMDZ.COMMON/libf/evolution/deftank/run_PEM.def (modified) (4 diffs)

trunk/LMDZ.COMMON/libf/evolution/deftank/run_PEM.def

@@ -7 +7 @@
 
 #---------- Output parameters ----------#
-# Frequency of outputs for the PEM? Default = 1
-ecritpem=1
+# Output rate? Default = 1
+output_rate=1
 
 #---------- Orbital parameters ----------#
 # Do you want to follow an orbital forcing read in "obl_ecc_lsp.asc"? Default = .false.
-# evol_orbit_pem=.false.
+# evol_orbit=.false.
 
-# If evol_orbit_pem=.true., number of Earth years before present to start the PEM run. Default = 0
+# If evol_orbit=.true., number of Earth years before present to start the PEM run. Default = 0
 # year_earth_bp_ini=0
 
@@ -30 +30 @@
 
 #---------- Stopping criteria parameters ----------#
-# If evol_orbit_pem=.false., maximal number of iterations if no stopping criterion is reached. Default=100000000
-# Max_iter_pem=100000000
+# If evol_orbit=.false., maximal number of iterations if no stopping criterion is reached. Default=100000000.
+# nyears_max=100000000
 
 # Acceptance rate of sublimating H2O ice surface change. Default = 0.2
-# h2o_ice_crit=0.2
+# h2oice_crit=0.2
 
 # Acceptance rate of sublimating CO2 ice surface change. Default = 0.2
-# co2_ice_crit=0.2
+# co2ice_crit=0.2
 
 # Acceptance rate of pressure surface change. Default = 0.15
-# ps_criterion = 0.15
+# ps_crit = 0.15
 
 # Acceptance of change for obliquity. Default = 1.
@@ -53 +53 @@
 #---------- Subsurface parameters ----------#
 # Do you want to run with subsurface physical processes in the PEM? Default = .true.
-# soil_pem=.true.
+# do_soil=.true.
 
-# Do you want to run with adsoprtion in the PEM? Default = .false.
-# adsorption_pem=.false.
+# Do you want to run with adsorption/desorption in the PEM? Default = .false.
+# sorption=.false.
 
 # Do you want to modify the soil thermal properties with the pressure? Default = .false.
@@ -77 +77 @@
 
 #---------- 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=9200.
+# Amount of H2O ice to initialize the huge reservoir if the variable is not present in "startpem.nc". Default = 9200. kg.m-2 (= 10 m)
+# h2oice_huge_ini=9200.
 
 # Threshold to consider the amount of H2O ice as an infinite reservoir. Default = 460. kg.m-2 (= 0.5 m)
-# inf_h2oice_threshold=460.
+# h2oice_cap_threshold=460.
 
 # Do you want the H2O ice to flow along subslope inside a cell? Default = .true.