Changeset 3991 for trunk/LMDZ.COMMON/libf/evolution/pem.F90

Timestamp:

Dec 16, 2025, 4:39:24 PM (5 weeks ago)

Author:

jbclement

Message:

PEM:
Apply documentation template everywhere: standardized headers format with short description, separators between functions/subroutines, normalized code sections, aligned dependencies/arguments/variables declaration.
JBC

File:

• trunk/LMDZ.COMMON/libf/evolution/pem.F90 (modified) (5 diffs)

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

@@ -26 +26 @@
 
 PROGRAM pem
-
+!-----------------------------------------------------------------------
+! NAME
+!     pem
+!
+! DESCRIPTION
+!     Main entry point for the Planetary Evolution Model (PEM).
+!
+! AUTHORS & DATE
+!     R. Vandemeulebrouck, 2022
+!     L. Lange, 2022
+!     JB Clement, 2023-2025
+!
+! NOTES
+!     Drives initialization, run loop, and outputs for PEM/PCM coupling.
+!-----------------------------------------------------------------------
+
+! DEPENDENCIES
+! ------------
 use phyetat0_mod,          only: phyetat0
 use phyredem,              only: physdem0, physdem1
@@ -99 +116 @@
 #endif
 
+! DECLARATION
+! -----------
 implicit none
 
@@ -106 +125 @@
 include "iniprint.h"
 
+! LOCAL VARIABLES
+! ---------------
 ! Same variable names as in the PCM
 integer, parameter :: ngridmx = 2 + (jjm - 1)*iim - 1/jjm
@@ -256 +277 @@
 #endif
 
-! Loop variables
 integer :: i, l, ig, nnq, t, islope, ig_loop, islope_loop, isoil, icap
 real    :: totmass_ini
 logical :: num_str
 
-! Code
-!-----
+! CODE
+! ----
 ! Elapsed time with system clock
 call system_clock(count_rate = cr)
@@ -682 +702 @@
         call stopping_crit_h2o(ngrid,nslope,cell_area,delta_h2o_adsorbed,delta_h2o_icetablesublim,h2o_ice,d_h2oice,S_atm_2_h2o,S_h2o_2_atm,S_atm_2_h2oice,S_h2oice_2_atm,stopCrit)
         call balance_h2oice_reservoirs(ngrid,nslope,S_atm_2_h2o,S_h2o_2_atm,S_atm_2_h2oice,S_h2oice_2_atm,h2o_ice,d_h2oice,d_h2oice_new)
-    
+
         do islope = 1,nslope
             do ig = 1,ngrid