Changeset 4174 for trunk/LMDZ.COMMON/libf/evolution/backup.F90

Timestamp:

Apr 8, 2026, 11:36:16 AM (3 days ago)

Author:

jbclement

Message:

PEM:
Refactor of H2O flux balancing:

• centralizes flux balancing logic;
• removes intermediate variables, especially related to of H2O mass bookkeeping;
• adds new stopping criterion when H2O flux balance fails8 (sanity check);
• introduces configurable weighting strategies for H2O flux balancing.

JBC

File:

• trunk/LMDZ.COMMON/libf/evolution/backup.F90 (modified) (3 diffs)

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

@@ -70 +70 @@
 !=======================================================================
 SUBROUTINE 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,h2oice_depth,flux_ssice_avg,h2o_ads_reg,co2_ads_reg,layerings_map,backup_idt)
+                           icetable_depth,icetable_thickness,ice_porefilling,h2oice_depth,h2o_ads_reg,co2_ads_reg,layerings_map,backup_idt)
 !-----------------------------------------------------------------------
 ! NAME
@@ -107 +107 @@
 real(dp),       dimension(:),     intent(in)           :: ps_avg, ps_dev
 real(dp),                         intent(in)           :: ps_avg_glob, ps_avg_glob_ini
-real(dp),       dimension(:,:),   intent(in)           :: tsurf_avg, tsurf_dev, icetable_depth, icetable_thickness, h2oice_depth, flux_ssice_avg
+real(dp),       dimension(:,:),   intent(in)           :: tsurf_avg, tsurf_dev, icetable_depth, icetable_thickness, h2oice_depth
 real(dp),       dimension(:,:,:), intent(in)           :: tsoil_avg, tsoil_dev, ice_porefilling, h2o_ads_reg, co2_ads_reg
 type(layering), dimension(:,:),   intent(in)           :: layerings_map
@@ -115 +115 @@
 ! LOCAL VARIABLES
 ! ---------------
-real(dp),    dimension(ngrid,nslope)           :: h2o_ice4PCM, co2_ice4PCM, tsurf4PCM, flux_geo4PCM, albedo4PCM, emissivity4PCM, flux_ssice4PCM, h2oice_depth4PCM
+real(dp),    dimension(ngrid,nslope)           :: h2o_ice4PCM, co2_ice4PCM, tsurf4PCM, flux_geo4PCM, albedo4PCM, emissivity4PCM, h2oice_depth4PCM
 real(dp),    dimension(ngrid,nlayer)           :: teta4PCM, air_mass4PCM
 real(dp),    dimension(ngrid,nsoil_PCM,nslope) :: tsoil4PCM, inertiesoil4PCM