Changeset 4110 for trunk/LMDZ.COMMON/libf/evolution/clim_state_init.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/clim_state_init.F90 (modified) (15 diffs)

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

@@ -48 +48 @@
 use tracers,    only: nq, qnames, set_q_PCM
 use stoppage,   only: stop_clean
-use display,    only: print_msg
+use display,    only: print_msg, LVL_NFO
 
 ! DECLARATION
@@ -76 +76 @@
 inquire(file = start_name,exist = here)
 if (.not. here) call stop_clean(__FILE__,__LINE__,'cannot find required file "'//start_name//'"!',1)
-call print_msg('> Reading "'//start_name//'"')
+call print_msg('> Reading "'//start_name//'"',LVL_NFO)
 call open_nc(start_name,'read')
 
@@ -126 +126 @@
 use stoppage,   only: stop_clean
 use config,     only: read_callphys
-use display,    only: print_msg
+use display,    only: print_msg, LVL_NFO
 
 ! DECLARATION
@@ -147 +147 @@
 inquire(file = start1D_name,exist = here)
 if (.not. here) call stop_clean(__FILE__,__LINE__,'cannot find required file "'//start1D_name//'"!',1)
-call print_msg('> Reading "'//start1D_name//'"')
+call print_msg('> Reading "'//start1D_name//'"',LVL_NFO)
 open(newunit = funit,file = start1D_name,status = "old",action = "read",iostat = ierr)
 if (ierr /= 0) call stop_clean(__FILE__,__LINE__,'error opening file "'//start1D_name//'"!',ierr)
@@ -201 +201 @@
 use geometry,  only: ngrid, nslope, nsoil_PCM
 use stoppage,  only: stop_clean
-use config,    only: read_controldata
+use config,    only: read_control_data
 use slopes,    only: set_def_slope_mean, set_subslope_dist, set_iflat
 use surface,   only: set_albedodat_PCM, set_albedo_PCM, set_emissivity_PCM
@@ -209 +209 @@
 use surf_ice,  only: set_is_h2o_perice_PCM, set_co2_perice_PCM
 use soil,      only: set_TI_PCM, set_inertiedat_PCM
-use display,   only: print_msg
+use display,   only: print_msg, LVL_NFO
 
 ! DECLARATION
@@ -228 +228 @@
 
 ! Allocate the array to store the variables
-call print_msg('> Reading "'//startfi_name//'"')
+call print_msg('> Reading "'//startfi_name//'"',LVL_NFO)
 allocate(tmp1d(nslope + 1),tmp2d(ngrid,nslope))
 
 ! Get control data
-call read_controldata()
+call read_control_data()
 
 ! Open
@@ -308 +308 @@
 !
 ! DESCRIPTION
-!     Read the file "startevol.nc" which stores the PEM state.
+!     Read the file "startevo.nc" which stores the PEM state.
 !
 ! AUTHORS & DATE
@@ -337 +337 @@
 use slopes,             only: subslope_dist, def_slope_mean
 use maths,              only: pi
-use display,            only: print_msg
+use display,            only: print_msg, LVL_WRN, LVL_NFO, LVL_ERR
 use utility,            only: int2str
 
@@ -367 +367 @@
 integer(di)                               :: i, islope, k, nb_str_max, nsoil_startpem
 real(dp)                                  :: delta           ! Depth of the interface regolith/breccia, breccia/bedrock [m]
-real(dp), dimension(ngrid,nsoil,nslope)   :: TI_startpem     ! Soil thermal inertia saved in the startevol [SI]
-real(dp), dimension(ngrid,nsoil,nslope)   :: tsoil_startpem  ! Soil temperature saved in the startevol [K]
+real(dp), dimension(ngrid,nsoil,nslope)   :: TI_startpem     ! Soil thermal inertia saved in the startevo [SI]
+real(dp), dimension(ngrid,nsoil,nslope)   :: tsoil_startpem  ! Soil temperature saved in the startevo [K]
 real(dp), dimension(:,:,:,:), allocatable :: layerings_array ! Array for layerings
 
@@ -380 +380 @@
 if (.not. here) then
     ! It is possibly because it is the very first PEM run so everything needs to be initalized
-    call print_msg('> The file "'//startpem_name//'" was not found (possibly because this is the first PEM run)')
+    call print_msg('> The file "'//startpem_name//'" was not found (possibly because this is the first PEM run)',LVL_WRN)
 
     ! H2O ice
-    call print_msg("'h2o_ice' is initialized with default value 'h2oice_huge_ini' where 'watercaptag' is true and where yearly minimum of frost can be considered as a huge reservoir ('threshold_h2oice_cap').")
+    call print_msg("'h2o_ice' is initialized with default value 'h2oice_huge_ini' where 'watercaptag' is true and where yearly minimum of frost can be considered as a huge reservoir ('threshold_h2oice_cap').",LVL_NFO)
     h2o_ice(:,:) = 0._dp
     do i = 1,ngrid
@@ -394 +394 @@
 
     ! CO2 ice
-    call print_msg("'co2_ice' is initialized with 'perennial_co2ice' and yearly minimum of frost found in the PCM.")
+    call print_msg("'co2_ice' is initialized with 'perennial_co2ice' and yearly minimum of frost found in the PCM.",LVL_NFO)
     co2_ice(:,:) = co2_perice_PCM(:,:) + co2frost_PCM(:,:) - co2_frost4PCM(:,:)
 
@@ -494 +494 @@
     ! Layering
     if (do_layering) then
-        call print_msg('layerings_map is initialized with sub-surface strata.')
-        call print_msg("Ice is added with 'h2oice_huge_ini' where 'watercaptag' is true and otherwise with 'perennial_co2ice' found in the PCM.")
+        call print_msg('layerings_map is initialized with sub-surface strata.',LVL_NFO)
+        call print_msg("Ice is added with 'h2oice_huge_ini' where 'watercaptag' is true and otherwise with 'perennial_co2ice' found in the PCM.",LVL_NFO)
         do i = 1,ngrid
             if (is_h2o_perice_PCM(i)) then
@@ -516 +516 @@
 ! If the file is present
 ! ~~~~~~~~~~~~~~~~~~~~~~
-call print_msg('> Reading "'//startpem_name//'"')
+call print_msg('> Reading "'//startpem_name//'"',LVL_NFO)
 call open_nc(startpem_name,'read')
 
@@ -531 +531 @@
     call get_dim_nc('subsurface_layers',nsoil_startpem)
     if (nsoil_startpem /= nsoil) then
-        call print_msg('nsoil (PEM) = '//int2str(nsoil)//' | nsoil ("'//startpem_name//'") = '//int2str(nsoil_startpem))
+        call print_msg('nsoil (PEM) = '//int2str(nsoil)//' | nsoil ("'//startpem_name//'") = '//int2str(nsoil_startpem),LVL_ERR)
         call stop_clean(__FILE__,__LINE__,'nsoil defined in the PEM is different from the one read in "'//startpem_name//'"!',1)
     end if