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

Timestamp:

Mar 16, 2026, 3:24:12 PM (13 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/pem.F90 (modified) (10 diffs)

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

@@ -23 +23 @@
 ! PEM modules
 use allocation,         only: ini_allocation, end_allocation
-use atmosphere,         only: ps_PCM, evolve_pressure, CO2cond_ps_PCM, build4PCM_atmosphere
+use atmosphere,         only: ps_PCM, evolve_pressure, CO2cond_ps_PCM
+use backup,             only: save_clim_state, backup_rate
 use clim_state_init,    only: read_start, read_startfi, read_startpem
-use clim_state_rec,     only: write_restart, write_restartfi, write_restartpem
 use config,             only: read_rundef, read_display_config
 use display,            only: print_ini, print_end, print_msg, LVL_NFO, LVL_WRN
@@ -39 +39 @@
 use physics,            only: g
 use slopes,             only: subslope_dist, def_slope_mean
-use soil,               only: do_soil, set_soil, TI, build4PCM_soil
+use soil,               only: do_soil, set_soil, TI
 use soil_temp,          only: tsoil_PCM, shift_tsoil2surf, evolve_soil_temp
 use soil_therm_inertia, only: update_soil_TI
 use sorption,           only: do_sorption, compute_totmass_adsorbed, evolve_regolith_adsorption
 use stopping_crit,      only: stopFlags, stopping_crit_pressure, stopping_crit_h2o, stopping_crit_h2oice, stopping_crit_co2ice
-use surface,            only: emissivity_PCM, build4PCM_surf_rad_prop
-use surf_ice,           only: evolve_co2ice, evolve_h2oice, balance_h2oice_reservoirs, build4PCM_perice
-use surf_temp,          only: tsurf_PCM, adapt_tsurf2disappearedice, build4PCM_tsurf
+use surface,            only: emissivity_PCM
+use surf_ice,           only: evolve_co2ice, evolve_h2oice, balance_h2oice_reservoirs
+use surf_temp,          only: tsurf_PCM, adapt_tsurf2disappearedice
 use tendencies,         only: compute_tendice, evolve_tend_co2ice, evolve_tend_h2oice
-use tracers,            only: adapt_tracers2pressure, build4PCM_tracers, nq
+use tracers,            only: adapt_tracers2pressure
 use utility,            only: real2str
 use workflow_status,    only: i_pem_run, read_workflow_status, update_workflow_status
@@ -71 +71 @@
 real(dp)                              :: ps_avg_glob_old ! Global average pressure of previous time step [Pa]
 real(dp)                              :: ps_avg_glob     ! Global average pressure of current time step [Pa]
-real(dp), dimension(:),   allocatable :: ps4PCM          ! Surface pressure reconstruction to feed back into PCM [Pa]
-real(dp), dimension(:,:), allocatable :: teta4PCM        ! Potential temperature reconstruction to feed back into PCM [K]
-real(dp), dimension(:,:), allocatable :: air_mass4PCM    ! Air mass reconstruction to feed back into PCM [kg]
-real(dp)                              :: pa4PCM          ! Transition pressure (for hybrid coord.) [Pa]
-real(dp)                              :: preff4PCM       ! Reference pressure [Pa]
 ! Ice-related:
 real(dp),    dimension(:,:),   allocatable :: h2o_ice                    ! H2O ice [kg.m-2]
@@ -89 +84 @@
 real(dp),    dimension(:,:,:), allocatable :: minPCM_h2ofrost               ! Minimum of H2O frost over the last PCM year [kg.m-2]
 real(dp),    dimension(:,:,:), allocatable :: minPCM_co2frost               ! Minimum of CO2 frost over the last PCM year [kg.m-2]
-real(dp),    dimension(:,:),   allocatable :: h2o_ice4PCM                ! H2O ice reconstruction to feed back into PCM [kg.m-2]
-real(dp),    dimension(:,:),   allocatable :: co2_ice4PCM                ! CO2 ice reconstruction to feed back into PCM [kg.m-2]
-logical(k4), dimension(:),     allocatable :: is_h2o_perice              ! Location of H2O infinite reservoir, called 'watercaptag' in PCM
 logical(k4), dimension(:,:),   allocatable :: is_co2ice_disappeared      ! Flag to check if CO2 ice disappeared at the previous timestep
 ! Surface-related:
@@ -100 +92 @@
 real(dp), dimension(:,:), allocatable :: zshift_surf         ! Elevation shift for the surface [m]
 real(dp), dimension(:,:), allocatable :: zlag                ! Newly built lag thickness [m]
-real(dp), dimension(:,:), allocatable :: albedo4PCM          ! Albedo reconstruction to feed back into PCM
-real(dp), dimension(:,:), allocatable :: emissivity4PCM      ! Emissivity reconstruction to feed back into PCM
-real(dp), dimension(:,:), allocatable :: tsurf4PCM           ! Surface temperature reconstruction to feed back into PCM [K]
 ! Soil-related:
 real(dp), dimension(:,:,:),   allocatable :: tsoil_avg                ! Average soil temperature [K]
@@ -112 +101 @@
 real(dp), dimension(:),       allocatable :: delta_h2o_ads            ! Quantity of H2O exchanged due to adsorption/desorption [kg/m^2]
 real(qp)                                  :: totmass_adsh2o           ! Total mass of H2O exchanged because of adsorption/desorption [kg]
-real(dp), dimension(:,:,:),   allocatable :: inertiesoil4PCM          ! Soil thermal inertia reconstruction to feed back into PCM [SI]
-real(dp), dimension(:,:,:),   allocatable :: tsoil4PCM                ! Average soil temperature reconstruction to feed back into PCM [K]
-real(dp), dimension(:,:),     allocatable :: flux_geo4PCM             ! Geothermal flux reconstruction to feed back into PCM [W/m2]
 real(dp), dimension(:,:,:),   allocatable :: h2o_ads_reg              ! H2O adsorbed in the regolith [kg/m^2]
 real(dp), dimension(:,:,:),   allocatable :: co2_ads_reg              ! CO2 adsorbed in the regolith [kg/m^2]
@@ -123 +109 @@
 real(dp), dimension(:,:),   allocatable :: q_co2_ts_ini ! Initial CO2 mass mixing ratio in the first layer [kg/kg]
 real(dp), dimension(:,:),   allocatable :: q_h2o_ts     ! H2O mass mixing ratio in the first layer [kg/kg]
-real(dp), dimension(:,:,:), allocatable :: q4PCM        ! Tracers reconstruction to feed back into PCM [kg/kg]
 ! Tendency-related:
 real(dp), dimension(:,:), allocatable :: d_co2ice     ! Tendency of perennial CO2 ice [kg/m2/y]
@@ -440 +425 @@
     idt = idt + 1
 
+    ! Save periodic backups of restart files
+    if (backup_rate > 0 .and. mod(idt,backup_rate) == 0) then
+        call save_clim_state(h2o_ice,co2_ice,tsurf_avg,tsurf_dev,tsoil_avg,tsoil_dev,ps_avg,ps_dev,ps_avg_glob,ps_avg_glob_ini, &
+                             icetable_depth,icetable_thickness,ice_porefilling,h2o_ads_reg,co2_ads_reg,layerings_map,idt)
+    end if
+
     ! Check the stopping criteria
     call print_msg("> Checking the stopping criteria",LVL_NFO)
@@ -469 +460 @@
 if (evo_orbit) call update_orbit(n_yr_sim,n_yr_run)
 
-! Build ice for the PCM
-allocate(h2o_ice4PCM(ngrid,nslope),co2_ice4PCM(ngrid,nslope),is_h2o_perice(ngrid))
-call build4PCM_perice(h2o_ice,co2_ice,is_h2o_perice,h2o_ice4PCM,co2_ice4PCM)
-
-! Build surface temperature for the PCM
-allocate(tsurf4PCM(ngrid,nslope))
-call build4PCM_tsurf(tsurf_avg,tsurf_dev,tsurf4PCM)
-
-! Build soil for the PCM
-if (do_soil) then
-    allocate(tsoil4PCM(ngrid,nsoil_PCM,nslope),inertiesoil4PCM(ngrid,nsoil_PCM,nslope),flux_geo4PCM(ngrid,nslope))
-    call build4PCM_soil(tsoil_avg,tsoil_dev,inertiesoil4PCM,tsoil4PCM,flux_geo4PCM)
-end if
-
-! Build atmosphere for the PCM
-allocate(ps4PCM(ngrid),teta4PCM(ngrid,nlayer),air_mass4PCM(ngrid,nlayer))
-call build4PCM_atmosphere(ps_avg,ps_dev,ps_avg_glob,ps_avg_glob_ini,ps4PCM,pa4PCM,preff4PCM,teta4PCM,air_mass4PCM)
-
-! Build tracers for the PCM
-allocate(q4PCM(ngrid,nlayer,nq))
-call build4PCM_tracers(ps4PCM,q4PCM)
-
-! Build surface radiative properties state for the PCM
-allocate(albedo4PCM(ngrid,nslope),emissivity4PCM(ngrid,nslope))
-call build4PCM_surf_rad_prop(h2o_ice,co2_ice,albedo4PCM,emissivity4PCM)
-
-! Write the "startevo.nc"
-call write_restartpem(h2o_ice,co2_ice,tsoil_avg,TI,icetable_depth,icetable_thickness,ice_porefilling,h2o_ads_reg,co2_ads_reg,layerings_map)
-
-! Write the "startfi.nc"
-call write_restartfi(is_h2o_perice,h2o_ice4PCM,co2_ice4PCM,tsurf4PCM,tsoil4PCM,inertiesoil4PCM,albedo4PCM,emissivity4PCM,flux_geo4PCM)
-
-! Write the "start.nc"
-call write_restart(ps4PCM,pa4PCM,preff4PCM,q4PCM,teta4PCM,air_mass4PCM)
+call save_clim_state(h2o_ice,co2_ice,tsurf_avg,tsurf_dev,tsoil_avg,tsoil_dev,ps_avg,ps_dev,ps_avg_glob,ps_avg_glob_ini, &
+                     icetable_depth,icetable_thickness,ice_porefilling,h2o_ads_reg,co2_ads_reg,layerings_map)
 
 ! Update the duration information of the workflow
@@ -508 +467 @@
 
 ! Deallocation
-deallocate(emissivity4PCM,albedo4PCM)
-deallocate(q4PCM)
-deallocate(ps4PCM,teta4PCM,air_mass4PCM)
-if (do_soil) deallocate(tsoil4PCM,inertiesoil4PCM,flux_geo4PCM)
-deallocate(tsurf4PCM)
-deallocate(co2_ice4PCM,h2o_ice4PCM,is_h2o_perice)
 deallocate(num)
 if (do_layering) then