Changeset 4110 for trunk/LMDZ.COMMON/libf/evolution/layered_deposits.F90

Timestamp:

Mar 9, 2026, 10:29:53 AM (2 weeks ago)

Author:

jbclement

Message:

PEM:

• Introduction of a configurable display/logging system with options 'out2term', 'out2log', 'verbosity_lvl'. All messages now use verbosity levels ('LVL_NFO', 'LVL_WRN', 'LVL_ERR' and 'LVL_DBG').
• Code encapsulation improvements with systematic privacy/protection of module variables.
• Addition of workflow safety checks for required executables, dependencies (e.g. 'ncdump'), input files and callphys keys.
• Renaming of PEM starting and diagnostic files ("startevol.nc" into "startevo.nc", "diagevol.nc" into "diagevo.nc").

JBC

File:

• trunk/LMDZ.COMMON/libf/evolution/layered_deposits.F90 (modified) (11 diffs)

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

@@ -28 +28 @@
 ! ---------
 ! Numerical parameter
-logical(k4),  protected :: do_layering
+logical(k4), protected :: do_layering
 
 ! Physical parameters
-logical(k4), protected :: impose_dust_ratio            ! Impose dust-to-ice ratio instead of dust tendency (see 'dust2ice_ratio')
-real(dp),    protected :: dust2ice_ratio               ! Dust-to-ice ratio when ice condenses (10% by default)
-real(dp),    protected :: d_dust                       ! Tendency of dust [kg.m-2.y-1]
-real(dp),    parameter :: dry_lag_porosity  = 0.2_dp   ! Porosity of dust lag
-real(dp),    parameter :: wet_lag_porosity  = 0.15_dp  ! Porosity of water ice lag
-real(dp),    parameter :: regolith_porosity = 0.4_dp   ! Porosity of regolith
-real(dp),    parameter :: breccia_porosity  = 0.1_dp   ! Porosity of breccia
-real(dp),    parameter :: co2ice_porosity   = 0.05_dp  ! Porosity of CO2 ice
-real(dp),    parameter :: h2oice_porosity   = 0.1_dp   ! Porosity of H2O ice
-real(dp),    parameter :: rho_dust          = 2500._dp ! Density of dust [kg.m-3]
-real(dp),    parameter :: rho_rock          = 3200._dp ! Density of rock [kg.m-3]
-real(dp),    parameter :: h_patchy_dust     = 5.e-4_dp ! Height under which a dust lag is considered patchy [m]
-real(dp),    parameter :: h_patchy_ice      = 5.e-4_dp ! Height under which a H2O ice lag is considered patchy [m]
+logical(k4), protected, private :: impose_dust_ratio            ! Impose dust-to-ice ratio instead of dust tendency (see 'dust2ice_ratio')
+real(dp),    protected, private :: dust2ice_ratio               ! Dust-to-ice ratio when ice condenses (10% by default)
+real(dp),    protected, private :: d_dust                       ! Tendency of dust [kg.m-2.y-1]
+real(dp),    parameter, private :: dry_lag_porosity  = 0.2_dp   ! Porosity of dust lag
+real(dp),    parameter, private :: wet_lag_porosity  = 0.15_dp  ! Porosity of water ice lag
+real(dp),    parameter, private :: regolith_porosity = 0.4_dp   ! Porosity of regolith
+real(dp),    parameter, private :: breccia_porosity  = 0.1_dp   ! Porosity of breccia
+real(dp),    parameter, private :: co2ice_porosity   = 0.05_dp  ! Porosity of CO2 ice
+real(dp),    parameter, private :: h2oice_porosity   = 0.1_dp   ! Porosity of H2O ice
+real(dp),    parameter, private :: rho_dust          = 2500._dp ! Density of dust [kg.m-3]
+real(dp),    parameter, private :: rho_rock          = 3200._dp ! Density of rock [kg.m-3]
+real(dp),    parameter, private :: h_patchy_dust     = 5.e-4_dp ! Height under which a dust lag is considered patchy [m]
+real(dp),    parameter, private :: h_patchy_ice      = 5.e-4_dp ! Height under which a H2O ice lag is considered patchy [m]
 
 ! Parameters for CO2 flux correction (inspired by Levrard et al. 2007)
-real(dp), parameter :: hmin_lag_co2 = 0.5_dp ! Minimal height of the lag deposit to reduce the sublimation rate [m]
-real(dp), parameter :: fred_sublco2 = 0.1_dp ! Reduction factor of sublimation rate
+real(dp), parameter, private :: hmin_lag_co2 = 0.5_dp ! Minimal height of the lag deposit to reduce the sublimation rate [m]
+real(dp), parameter, private :: fred_sublco2 = 0.1_dp ! Reduction factor of sublimation rate
 
 ! DATA STRUCTURES
@@ -99 +99 @@
 ! DEPENDENCIES
 ! ------------
-use utility, only: real2str, bool2str
-use display, only: print_msg
+use utility,  only: real2str, bool2str
+use display,  only: print_msg, LVL_NFO
+use stoppage, only: stop_clean
 
 ! DECLARATION
@@ -117 +118 @@
 d_dust = d_dust_in
 dust2ice_ratio = dust2ice_ratio_in
-call print_msg('do_layering       = '//bool2str(do_layering))
-call print_msg('impose_dust_ratio = '//bool2str(impose_dust_ratio))
-call print_msg('d_dust            = '//real2str(d_dust))
-call print_msg('dust2ice_ratio    = '//real2str(dust2ice_ratio))
+call print_msg('do_layering       = '//bool2str(do_layering),LVL_NFO)
+call print_msg('impose_dust_ratio = '//bool2str(impose_dust_ratio),LVL_NFO)
+call print_msg('d_dust            = '//real2str(d_dust),LVL_NFO)
+call print_msg('dust2ice_ratio    = '//real2str(dust2ice_ratio),LVL_NFO)
+if (dust2ice_ratio < 0._dp .or. dust2ice_ratio > 1._dp) call stop_clean(__FILE__,__LINE__,'''dust2ice_ratio'' outside admissible range [0.,1.]!',1)
 
 END SUBROUTINE set_layered_deposits_config
@@ -144 +146 @@
 ! ------------
 use utility, only: real2str
-use display, only: print_msg
+use display, only: print_msg, LVL_DBG
 
 ! DECLARATION
@@ -156 +158 @@
 ! CODE
 ! ----
-call print_msg('Top elevation       [m]: '//real2str(this%top_elevation))
-call print_msg('Height of CO2 ice   [m]: '//real2str(this%h_co2ice))
-call print_msg('Height of H2O ice   [m]: '//real2str(this%h_h2oice))
-call print_msg('Height of dust      [m]: '//real2str(this%h_dust))
-call print_msg('Height of pores     [m]: '//real2str(this%h_pore))
-call print_msg('Pore-filled ice [m3/m3]: '//real2str(this%poreice_volfrac))
+call print_msg('Top elevation       [m]: '//real2str(this%top_elevation),LVL_DBG)
+call print_msg('Height of CO2 ice   [m]: '//real2str(this%h_co2ice),LVL_DBG)
+call print_msg('Height of H2O ice   [m]: '//real2str(this%h_h2oice),LVL_DBG)
+call print_msg('Height of dust      [m]: '//real2str(this%h_dust),LVL_DBG)
+call print_msg('Height of pores     [m]: '//real2str(this%h_pore),LVL_DBG)
+call print_msg('Pore-filled ice [m3/m3]: '//real2str(this%poreice_volfrac),LVL_DBG)
 
 END SUBROUTINE print_stratum
@@ -529 +531 @@
 ! ------------
 use utility, only: int2str
-use display, only: print_msg
+use display, only: print_msg, LVL_DBG
 
 ! DECLARATION
@@ -550 +552 @@
 i = this%nb_str
 do while (associated(current))
-    call print_msg('Stratum '//int2str(i))
+    call print_msg('Stratum '//int2str(i),LVL_DBG)
     call print_stratum(current)
-    call print_msg('')
+    call print_msg('',LVL_DBG)
     current => current%down
     i = i - 1
@@ -822 +824 @@
 ! ------------
 use utility, only: int2str
-use display, only: print_msg
+use display, only: print_msg, LVL_DBG
 
 ! DEPENDENCIES
@@ -843 +845 @@
 ! ----
 do ig = 1,ngrid
-    call print_msg('Grid cell '//int2str(ig))
+    call print_msg('Grid cell '//int2str(ig),LVL_DBG)
     do islope = 1,nslope
-        call print_msg('Slope number '//int2str(islope))
+        call print_msg('Slope number '//int2str(islope),LVL_DBG)
         call print_layering(this(ig,islope))
-        call print_msg('~~~~~~~~~~')
+        call print_msg('~~~~~~~~~~',LVL_DBG)
     end do
-    call print_msg('--------------------')
+    call print_msg('--------------------',LVL_DBG)
 end do
 
@@ -1247 +1249 @@
 use stoppage, only: stop_clean
 use utility,  only: real2str
-use display,  only: print_err
+use display,  only: print_msg, LVL_ERR
 
 ! DECLARATION
@@ -1561 +1563 @@
 
 if (unexpected) then
-    call print_err('    dh_co2ice [m/y] = '//real2str(dh_co2ice))
-    call print_err('    dh_h2oice [m/y] = '//real2str(dh_h2oice))
-    call print_err('    dh_dust   [m/y] = '//real2str(dh_dust))
+    call print_msg('    dh_co2ice [m/y] = '//real2str(dh_co2ice),LVL_ERR)
+    call print_msg('    dh_h2oice [m/y] = '//real2str(dh_h2oice),LVL_ERR)
+    call print_msg('    dh_dust   [m/y] = '//real2str(dh_dust),LVL_ERR)
     call stop_clean(__FILE__,__LINE__,'situation for the layering construction not managed!',1)
 end if