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

Timestamp:

Feb 16, 2026, 10:28:56 AM (12 days ago)

Author:

jbclement

Message:

PEM:

• Making the computation of ice tendencies more reliable by doing it after 'read_startpem' to know exactly where perennial ice is (no matter if there is a "startpem.nc" or not). Moreover, when there is no "startpem.nc", location of perennial ice depends now on 'watercaptag' and on huge amount of frost. This prevents negative ice tendency while there is no ice which can happen with weird PCM inputs (i.e. 'watercaptag' = F & 'watercap' /= 0 & no "stratpem.nc").
• Few small safeguards throughout the code.

JBC

File:

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

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

@@ -121 +121 @@
 ! DEPENDENCIES
 ! ------------
-use geometry,   only: nlayer, nslope
+use geometry,   only: nlayer
 use atmosphere, only: set_ap, set_bp, set_ps_PCM, set_teta_PCM, set_u_PCM, set_v_PCM, compute_alt_coord, set_ap, set_bp
 use tracers,    only: nq, set_q_PCM
@@ -134 +134 @@
 ! LOCAL VARIABLES
 ! ---------------
-integer(di)                     :: i, j, k, ierr, funit
+integer(di)                     :: i, k, ierr, funit
 character(30)                   :: header
 real(dp), dimension(1,1)        :: tmp
@@ -326 +326 @@
 use planet,             only: TI_breccia, TI_bedrock, alpha_clap_h2o, beta_clap_h2o
 use frost,              only: h2ofrost_PCM, h2o_frost4PCM, co2frost_PCM, co2_frost4PCM
-use surf_ice,           only: is_h2o_perice_PCM, h2oice_huge_ini, co2_perice_PCM
+use surf_ice,           only: is_h2o_perice_PCM, h2oice_huge_ini, co2_perice_PCM, threshold_h2oice_cap
 use soil,               only: do_soil, TI, depth_breccia, depth_bedrock, index_breccia, index_bedrock, inertiedat, layer, inertiedat_PCM
 use soil_temp,          only: ini_tsoil_pem, compute_tsoil
@@ -335 +335 @@
 use layered_deposits,   only: layering, do_layering, array2map, ini_layering, add_stratum
 use surf_ice,           only: rho_co2ice, rho_h2oice
+use slopes,             only: subslope_dist, def_slope_mean
+use maths,              only: pi
 use display,            only: print_msg
 use utility,            only: int2str
@@ -344 +346 @@
 ! ARGUMENTS
 ! ---------
-real(dp), dimension(:,:),       intent(in)    :: tsurf_avg_yr1, tsurf_avg_yr2 ! Surface temperature at the last but one PCM run and at the last PCM run
-real(dp),                       intent(in)    :: ps_avg_glob                  ! Mean average pressure on the planet
-real(dp), dimension(:,:),       intent(in)    :: ps_ts                        ! Surface pressure timeseries
-real(dp), dimension(:,:),       intent(in)    :: q_co2_ts, q_h2o_ts           ! MMR of CO2 and H2O
-real(dp), dimension(:,:),       intent(in)    :: h2o_surfdensity_avg          ! Average of surface water density
-real(dp), dimension(:,:),       intent(in)    :: d_h2oice, d_co2ice           ! Tendencies on ice
-real(dp), dimension(:,:),       intent(inout) :: h2o_ice, co2_ice             ! Surface ice
-real(dp), dimension(:,:,:),     intent(inout) :: tsoil_avg                    ! Soil temperature
-real(dp), dimension(:,:,:),     intent(inout) :: h2o_soildensity_avg          ! Average of soil water soil density
-real(dp), dimension(:,:),       intent(inout) :: icetable_depth               ! Ice table depth
-real(dp), dimension(:,:),       intent(inout) :: icetable_thickness           ! Ice table thickness
-real(dp), dimension(:,:,:),     intent(inout) :: ice_porefilling              ! Subsurface ice pore filling
-real(dp), dimension(:,:,:),     intent(inout) :: h2o_ads_reg, co2_ads_reg     ! Mass of H2O and CO2 adsorbed
-type(layering), dimension(:,:), intent(inout) :: layerings_map                ! Layerings
-real(dp), dimension(:),         intent(inout) :: delta_h2o_ads, delta_co2_ads ! Mass of H2O and CO2 exchanged due to adsorption/desorption
+real(dp),       dimension(:,:),   intent(in)    :: tsurf_avg_yr1, tsurf_avg_yr2 ! Surface temperature at the last but one PCM run and at the last PCM run
+real(dp),                         intent(in)    :: ps_avg_glob                  ! Mean average pressure on the planet
+real(dp),       dimension(:,:),   intent(in)    :: ps_ts                        ! Surface pressure timeseries
+real(dp),       dimension(:,:),   intent(in)    :: q_co2_ts, q_h2o_ts           ! MMR of CO2 and H2O
+real(dp),       dimension(:,:),   intent(in)    :: h2o_surfdensity_avg          ! Average of surface water density
+real(dp),       dimension(:,:),   intent(in)    :: d_h2oice, d_co2ice           ! Tendencies on ice
+real(dp),       dimension(:,:),   intent(inout) :: h2o_ice, co2_ice             ! Surface ice
+real(dp),       dimension(:,:,:), intent(inout) :: tsoil_avg                    ! Soil temperature
+real(dp),       dimension(:,:,:), intent(inout) :: h2o_soildensity_avg          ! Average of soil water soil density
+real(dp),       dimension(:,:),   intent(inout) :: icetable_depth               ! Ice table depth
+real(dp),       dimension(:,:),   intent(inout) :: icetable_thickness           ! Ice table thickness
+real(dp),       dimension(:,:,:), intent(inout) :: ice_porefilling              ! Subsurface ice pore filling
+real(dp),       dimension(:,:,:), intent(inout) :: h2o_ads_reg, co2_ads_reg     ! Mass of H2O and CO2 adsorbed
+type(layering), dimension(:,:),   intent(inout) :: layerings_map                ! Layerings
+real(dp),       dimension(:),     intent(inout) :: delta_h2o_ads, delta_co2_ads ! Mass of H2O and CO2 exchanged due to adsorption/desorption
 
 ! LOCAL VARIABLES
@@ -381 +383 @@
 
     ! H2O ice
-    call print_msg("'h2o_ice' is initialized with default value 'h2oice_huge_ini' where 'watercaptag' is true.")
+    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').")
     h2o_ice(:,:) = 0._dp
     do i = 1,ngrid
-        if (is_h2o_perice_PCM(i)) h2o_ice(i,:) = h2oice_huge_ini + h2ofrost_PCM(i,:) - h2o_frost4PCM(i,:)
+        if (is_h2o_perice_PCM(i)) then
+            h2o_ice(i,:) = h2oice_huge_ini
+        else if (sum((h2ofrost_PCM(i,:) - h2o_frost4PCM(i,:))*subslope_dist(i,:)/cos(pi*def_slope_mean(:)/180._dp)) > threshold_h2oice_cap) then
+            h2o_ice(i,:) = h2oice_huge_ini
+        end if
     end do
 
     ! CO2 ice
-    call print_msg("'co2_ice' is initialized with 'perennial_co2ice' found in the PCM.")
+    call print_msg("'co2_ice' is initialized with 'perennial_co2ice' and yearly minimum of frost found in the PCM.")
     co2_ice(:,:) = co2_perice_PCM(:,:) + co2frost_PCM(:,:) - co2_frost4PCM(:,:)