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

Timestamp:

Mar 17, 2026, 4:22:44 PM (34 hours 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/ice_table.F90 (modified) (4 diffs)

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

@@ -28 +28 @@
 logical(k4), protected :: icetable_dynamic     ! Flag to compute the icetable depth according to the dynamic method
 
-! VARIABLES
-! ---------
-real(dp), allocatable, dimension(:,:)   :: icetable_depth     ! Depth of the ice table [m]
-real(dp), allocatable, dimension(:,:)   :: icetable_thickness ! Thickness of the ice table [m]
-real(dp), allocatable, dimension(:,:,:) :: ice_porefilling    ! Amount of porefilling in each layer in each grid [m^3/m^3]
-
 contains
 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
-
-!=======================================================================
-SUBROUTINE ini_ice_table()
-!-----------------------------------------------------------------------
-! NAME
-!     ini_ice_table
-!
-! DESCRIPTION
-!     Allocate module arrays for ice table fields.
-!
-! AUTHORS & DATE
-!     L. Lange
-!     JB Clement, 2023-2025
-!
-! NOTES
-!
-!-----------------------------------------------------------------------
-
-! DEPENDENCIES
-! ------------
-use geometry, only: ngrid, nslope, nsoil
-
-! DECLARATION
-! -----------
-implicit none
-
-! CODE
-! ----
-allocate(icetable_depth(ngrid,nslope))
-allocate(icetable_thickness(ngrid,nslope))
-allocate(ice_porefilling(ngrid,nsoil,nslope))
-icetable_depth(:,:) = 0._dp
-icetable_thickness(:,:) = 0._dp
-ice_porefilling(:,:,:) = 0._dp
-
-END SUBROUTINE ini_ice_table
-!=======================================================================
-
-!=======================================================================
-SUBROUTINE end_ice_table
-!-----------------------------------------------------------------------
-! NAME
-!     end_ice_table
-!
-! DESCRIPTION
-!     Deallocate ice table arrays.
-!
-! AUTHORS & DATE
-!     L. Lange
-!     JB Clement, 2023-2025
-!
-! NOTES
-!
-!-----------------------------------------------------------------------
-
-! DECLARATION
-! -----------
-implicit none
-
-! CODE
-! ----
-if (allocated(icetable_depth)) deallocate(icetable_depth)
-if (allocated(icetable_thickness)) deallocate(icetable_thickness)
-if (allocated(ice_porefilling)) deallocate(ice_porefilling)
-
-END SUBROUTINE end_ice_table
-!=======================================================================
 
 !=======================================================================
@@ -256 +183 @@
 
 !=======================================================================
-SUBROUTINE compute_deltam_icetable(icetable_thickness_old,ice_porefilling_old,tsurf,tsoil,delta_icetable)
+SUBROUTINE compute_deltam_icetable(icetable_depth,icetable_thickness,ice_porefilling,icetable_thickness_old,ice_porefilling_old,tsurf,tsoil,delta_icetable)
 !-----------------------------------------------------------------------
 ! NAME
@@ -286 +213 @@
 ! ARGUMENTS
 ! ---------
-real(dp),    dimension(:,:),   intent(in)  :: icetable_thickness_old   ! Ice table thickness at the previous iteration [m]
-real(dp),    dimension(:,:,:), intent(in)  :: ice_porefilling_old      ! Ice pore filling at the previous iteration [m]
-real(dp),    dimension(:,:),   intent(in)  :: tsurf                    ! Surface temperature [K]
-real(dp),    dimension(:,:,:), intent(in)  :: tsoil                    ! Soil temperature [K]
-real(dp),    dimension(:),     intent(out) :: delta_icetable ! Mass of H2O ice that has been condensed on the ice table / sublimates from the ice table [kg/m^2]
+real(dp),    dimension(:,:),   intent(in)  :: icetable_depth         ! Ice table depth [m]
+real(dp),    dimension(:,:),   intent(in)  :: icetable_thickness     ! Ice table thickness [m]
+real(dp),    dimension(:,:,:), intent(in)  :: ice_porefilling        ! Ice pore filling [m]
+real(dp),    dimension(:,:),   intent(in)  :: icetable_thickness_old ! Ice table thickness at the previous iteration [m]
+real(dp),    dimension(:,:,:), intent(in)  :: ice_porefilling_old    ! Ice pore filling at the previous iteration [m]
+real(dp),    dimension(:,:),   intent(in)  :: tsurf                  ! Surface temperature [K]
+real(dp),    dimension(:,:,:), intent(in)  :: tsoil                  ! Soil temperature [K]
+real(dp),    dimension(:),     intent(out) :: delta_icetable         ! Mass of H2O ice that has been condensed on the ice table / sublimates from the ice table [kg/m^2]
 
 ! LOCAL VARIABLES
@@ -674 +604 @@
 
 ! Mass of H2O exchange between the ssi and the atmosphere
-call compute_deltam_icetable(icetable_thickness_old,ice_porefilling_old,tsurf_avg,tsoil_avg,delta_icetable)
+call compute_deltam_icetable(icetable_depth,icetable_thickness,ice_porefilling,icetable_thickness_old,ice_porefilling_old,tsurf_avg,tsoil_avg,delta_icetable)
 
 END SUBROUTINE evolve_ice_table