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

Timestamp:

Mar 9, 2026, 10:29:53 AM (4 days 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/surf_ice.F90 (modified) (8 diffs)

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

@@ -114 +114 @@
 ! DEPENDENCIES
 ! ------------
-use utility, only: real2str
-use display, only: print_msg
+use utility,  only: real2str
+use display,  only: print_msg, LVL_NFO
+use stoppage, only: stop_clean
 
 ! DECLARATION
@@ -129 +130 @@
 threshold_h2oice_cap = threshold_h2oice_cap_in
 h2oice_huge_ini = h2oice_huge_ini_in
-call print_msg('threshold_h2oice_cap = '//real2str(threshold_h2oice_cap))
-call print_msg('h2oice_huge_ini      = '//real2str(h2oice_huge_ini))
+call print_msg('threshold_h2oice_cap = '//real2str(threshold_h2oice_cap),LVL_NFO)
+call print_msg('h2oice_huge_ini      = '//real2str(h2oice_huge_ini),LVL_NFO)
+if (threshold_h2oice_cap < 0._dp) call stop_clean(__FILE__,__LINE__,'''threshold_h2oice_cap'' must be positive or null!',1)
+if (h2oice_huge_ini < 0._dp) call stop_clean(__FILE__,__LINE__,'''h2oice_huge_ini'' must be positive or null!',1)
 
 END SUBROUTINE set_surf_ice_config
@@ -223 +226 @@
 use slopes,   only: subslope_dist, def_slope_mean
 use maths,    only: pi
-use display,  only: print_msg
+use display,  only: print_msg, LVL_NFO
 
 ! DECLARATION
@@ -241 +244 @@
 ! CODE
 ! ----
-call print_msg('> Building surface ice/frost for the PCM')
+call print_msg('> Building surface ice/frost for the PCM',LVL_NFO)
 co2_ice4PCM(:,:) = co2_ice(:,:)
 h2o_ice4PCM(:,:) = 0._dp ! Because in the Mars PCM, only the variation of perennial H2O ice is monitored, not the absolute quantity
@@ -282 +285 @@
 use geometry,  only: ngrid, nslope
 use evolution, only: dt
-use display,   only: print_msg
+use display,   only: print_msg, LVL_NFO
 
 ! DECLARATION
@@ -301 +304 @@
 ! ----
 ! Evolution of CO2 ice for each physical point
-call print_msg('> Evolution of CO2 ice')
+call print_msg('> Evolution of CO2 ice',LVL_NFO)
 
 co2_ice_old = co2_ice
@@ -338 +341 @@
 use evolution,     only: dt
 use stopping_crit, only: stopping_crit_h2o, stopFlags
-use display,       only: print_msg
+use display,       only: print_msg, LVL_NFO
 
 ! DECLARATION
@@ -358 +361 @@
 ! CODE
 ! ----
-call print_msg('> Evolution of H2O ice')
+call print_msg('> Evolution of H2O ice',LVL_NFO)
 
 if (ngrid == 1) then ! In 1D