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

Timestamp:

Mar 17, 2026, 4:22:44 PM (3 weeks ago)

Author:

jbclement

Message:

PEM:

• Move ice table variables from "ice_table" to the main program.
• Merge "job_id_mod" and "job_timelimit_mod" into "job" which is relocated to the PEM folder.
• Rename local variables in procedures to avoid masking variables in parent scope.
• Few cleanings to delete remaining PEM-external "include" and "use".

JBC

File:

• trunk/LMDZ.COMMON/libf/evolution/pem.F90 (modified) (5 diffs)

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

@@ -19 +19 @@
 ! ------------
 ! Common modules
-use job_timelimit_mod, only: timelimit, antetime, timewall
 use parse_args_mod,    only: parse_args
-! PEM modules
+
 use allocation,         only: ini_allocation, end_allocation
 use atmosphere,         only: ps_PCM, evolve_pressure, CO2cond_ps_PCM
@@ -31 +30 @@
 use geometry,           only: ngrid, nslope, nday, nsoil_PCM, nsoil, cell_area, total_surface, nlayer
 use glaciers,           only: h2oice_flow, co2ice_flow, flow_co2glaciers, flow_h2oglaciers
-use ice_table,          only: icetable_equilibrium, icetable_dynamic, icetable_depth, icetable_thickness, ice_porefilling, evolve_ice_table
+use ice_table,          only: icetable_equilibrium, icetable_dynamic, evolve_ice_table
+use job,                only: timelimit, antetime, timewall
 use layered_deposits,   only: layering, do_layering, del_layering, evolve_layering, ptrarray, layering2surfice, surfice2layering, print_layering
 use maths,              only: pi
@@ -103 +103 @@
 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]
+real(dp), dimension(:,:),     allocatable :: icetable_depth           ! Depth of the ice table [m]
+real(dp), dimension(:,:),     allocatable :: icetable_thickness       ! Thickness of the ice table [m]
+real(dp), dimension(:,:,:),   allocatable :: ice_porefilling          ! Amount of porefilling [m^3/m^3]
 real(dp), dimension(:,:),     allocatable :: icetable_depth_old       ! Old depth of the ice table [m]
 real(dp), dimension(:),       allocatable :: delta_icetable ! Total mass of the H2O that has sublimated / condenses from the ice table [kg]
@@ -191 +194 @@
 ! Read the "startevo.nc"
 allocate(h2o_ice(ngrid,nslope),co2_ice(ngrid,nslope))
+allocate(icetable_depth(ngrid,nslope),icetable_thickness(ngrid,nslope),ice_porefilling(ngrid,nsoil,nslope))
 allocate(h2o_ads_reg(ngrid,nsoil,nslope),co2_ads_reg(ngrid,nsoil,nslope),delta_h2o_ads(ngrid),delta_co2_ads(ngrid))
 allocate(layerings_map(ngrid,nslope))
+icetable_depth(:,:) = 0._dp
+icetable_thickness(:,:) = 0._dp
+ice_porefilling(:,:,:) = 0._dp
 delta_h2o_ads(:) = 0._dp
 delta_co2_ads(:) = 0._dp
@@ -485 +492 @@
 deallocate(d_h2oice,d_co2ice,d_co2ice_ini)
 deallocate(delta_h2o_ads,delta_co2_ads,delta_icetable,icetable_depth_old)
+deallocate(icetable_depth,icetable_thickness,ice_porefilling)
 call end_allocation()