Changeset 3989 for trunk/LMDZ.COMMON/libf/evolution/sorption.F90

Timestamp:

Dec 11, 2025, 12:56:05 PM (5 weeks 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/sorption.F90 (moved) (moved from trunk/LMDZ.COMMON/libf/evolution/adsorption_mod.F90) (4 diffs)

trunk/LMDZ.COMMON/libf/evolution/sorption.F90

@@ -1 @@
-MODULE adsorption_mod
-
-implicit none
-
-logical                             :: adsorption_pem     ! True by default, to compute adsorption/desorption. Read in pem.def
+MODULE sorption
+
+implicit none
+
+logical                             :: do_sorption       ! Flag to compute adsorption/desorption
 real, allocatable, dimension(:,:,:) :: co2_adsorbed_phys ! co2 that is in the regolith [kg/m^2]
 real, allocatable, dimension(:,:,:) :: h2o_adsorbed_phys ! h2o that is in the regolith [kg/m^2]
@@ -85 +85 @@
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 
-use comsoil_h_PEM,         only: layer_PEM, index_breccia
+use soil,         only: layer_PEM, index_breccia
 use comslope_mod,          only: subslope_dist, def_slope_mean
 use vertical_layers_mod,   only: ap, bp
@@ -218 +218 @@
 SUBROUTINE regolith_co2adsorption(ngrid,nslope,nsoil_PEM,timelen,d_h2oglaciers,d_co2glaciers,waterice,co2ice,ps,q_co2,q_h2o,tsoil_PEM,TI_PEM,m_co2_completesoil,delta_mreg)
 
-use comsoil_h_PEM,         only: layer_PEM, index_breccia, index_breccia
+use soil,                  only: layer_PEM, index_breccia, index_breccia
 use comslope_mod,          only: subslope_dist, def_slope_mean
 use vertical_layers_mod,   only: ap, bp
@@ -349 +349 @@
 END SUBROUTINE regolith_co2adsorption
 
-END MODULE adsorption_mod
-
+END MODULE sorption