Changeset 4134 for trunk/LMDZ.COMMON/libf/evolution/clim_state_init.F90

Timestamp:

Mar 16, 2026, 3:24:12 PM (6 days ago)

Author:

jbclement

Message:

PEM:

• Add a periodic backup based on a single wrapper which builds PCM-compatible climate state and writes "restart" files.
• Remove condition "if (.not. allocated(*))" for more strict safeguard.

JBC

File:

• trunk/LMDZ.COMMON/libf/evolution/clim_state_init.F90 (modified) (6 diffs)

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

@@ -17 +17 @@
 ! ------------
 use numerics,  only: dp, di, k4
-use io_netcdf, only: open_nc, close_nc, get_var_nc, get_dim_nc, put_var_nc, start_name, start1D_name, startfi_name, startpem_name
+use io_netcdf, only: open_nc, close_nc, get_var_nc, get_dim_nc, put_var_nc, start_name, start1D_name, startfi_name, startevo_name
 
 ! DECLARATION
@@ -369 +369 @@
 ! ----
 ! Check if the file exists
-inquire(file = startpem_name,exist = here)
+inquire(file = startevo_name,exist = here)
 
 ! If the file is not here
@@ -375 +375 @@
 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)',LVL_WRN)
+    call print_msg('> The file "'//startevo_name//'" was not found (possibly because this is the first PEM run)',LVL_WRN)
 
     ! H2O ice
@@ -509 +509 @@
 ! If the file is present
 ! ~~~~~~~~~~~~~~~~~~~~~~
-call print_msg('> Reading "'//startpem_name//'"',LVL_NFO)
-call open_nc(startpem_name,'read')
+call print_msg('> Reading "'//startevo_name//'"',LVL_NFO)
+call open_nc(startevo_name,'read')
 
 ! H2O ice
@@ -524 +524 @@
     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),LVL_ERR)
-        call stop_clean(__FILE__,__LINE__,'nsoil defined in the PEM is different from the one read in "'//startpem_name//'"!',1)
+        call print_msg('nsoil (PEM) = '//int2str(nsoil)//' | nsoil ("'//startevo_name//'") = '//int2str(nsoil_startpem),LVL_ERR)
+        call stop_clean(__FILE__,__LINE__,'nsoil defined in the PEM is different from the one read in "'//startevo_name//'"!',1)
     end if
 
@@ -578 +578 @@
 
 ! Close
-call close_nc(startpem_name)
+call close_nc(startevo_name)
 
 END SUBROUTINE read_startpem