Changeset 3989 for trunk/LMDZ.COMMON/libf/evolution/glaciers.F90

Timestamp:

Dec 11, 2025, 12:56:05 PM (6 weeks ago)

Author:

jbclement

Message:

PEM:
Massive structural refactor of the PEM codebase for improved readability, consistency and maintainability. The goal is to modernize, standardize and consolidate the code while removing legacy complexity. In detail, this change:

• Performs large-scale cleanup removing unused code, obsolete routines, duplicated functionality and deprecated initialization logic;
• Removes "*_mod" wrappers;
• Replaces mixed naming conventions with clearer variable names, domain-based file/module names and purpose-based routine names;
• Adds native reading/writing capabilities to the PEM removing the cumbersome dependency on Mars PCM subroutines;
• Centralizes soil/ice/orbital/PEM configuration logic into dedicated modules;
• Simplifies routines and updates calls/interfaces to match the new structure.

This refactor significantly clarifies the codebase and provides a cleaner foundation for forthcoming developments.
JBC

File:

• trunk/LMDZ.COMMON/libf/evolution/glaciers.F90 (moved) (moved from trunk/LMDZ.COMMON/libf/evolution/glaciers_mod.F90) (5 diffs)

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

@@ -1 @@
-MODULE glaciers_mod
-
-implicit none
-
-! Different types of ice
-type :: ice
-    real :: h2o
-    real :: co2
-end type ice
-
-! Perennial ice manage by the PEM
-type(ice), dimension(:,:), allocatable :: perice
+MODULE glaciers
+
+implicit none
 
 ! Flags for ice flow
-logical :: h2oice_flow  ! True by default, to compute H2O ice flow. Read in "run_PEM.def"
-logical :: co2ice_flow  ! True by default, to compute CO2 ice flow. Read in "run_PEM.def"
-
-! Threshold to consider H2O ice as watercap
-real :: inf_h2oice_threshold   ! To consider the amount of H2O ice as an infinite reservoir [kg.m-2]
-
-! Ice properties
-real, parameter :: rho_co2ice = 1650. ! Density of CO2 ice [kg.m-3]
-real, parameter :: rho_h2oice = 920.  ! Density of H2O ice [kg.m-3]
+logical :: h2oice_flow  ! True by default, flag to compute H2O ice flow
+logical :: co2ice_flow  ! True by default, flag to compute CO2 ice flow
 
 !=======================================================================
 contains
-!=======================================================================
-
-SUBROUTINE set_perice4PCM(ngrid,nslope,PCMfrost,is_perh2oice,PCMh2oice,PCMco2ice)
-
-use metamorphism,   only: iPCM_h2ofrost
-use comslope_mod,   only: subslope_dist, def_slope_mean
-use phys_constants, only: pi
-
-
-implicit none
-
-! Arguments
-!----------
-integer, intent(in) :: ngrid, nslope
-real, dimension(:,:,:), intent(inout) :: PCMfrost
-logical, dimension(ngrid),        intent(out) :: is_perh2oice
-real,    dimension(ngrid,nslope), intent(out) :: PCMco2ice, PCMh2oice
-
-! Local variables
-!----------------
-integer :: i
-
-! Code
-!-----
-write(*,*) '> Reconstructing perennial ic for the PCM'
-PCMco2ice(:,:) = perice(:,:)%co2
-PCMh2oice(:,:) = 0. ! Because in the Mars PCM, only the variation of perennial H2O ice is monitored, not the absolute quantity
-do i = 1,ngrid
-    ! Is H2O ice still considered as an infinite reservoir for the PCM?
-    if (sum(perice(i,:)%h2o*subslope_dist(i,:)/cos(pi*def_slope_mean(:)/180.)) > inf_h2oice_threshold) then
-        ! There is enough ice to be considered as an infinite reservoir
-        is_perh2oice(i) = .true.
-    else
-        ! Too little ice to be considered as an infinite reservoir so ice is transferred to the frost
-        is_perh2oice(i) = .false.
-        PCMfrost(i,iPCM_h2ofrost,:) = PCMfrost(i,iPCM_h2ofrost,:) + perice(i,:)%h2o
-        perice(i,:)%h2o = 0.
-    endif
-enddo
-
-END SUBROUTINE set_perice4PCM
-!=======================================================================
-
-SUBROUTINE ini_ice(ngrid,nslope)
-
-implicit none
-
-! Arguments
-!----------
-integer, intent(in) :: ngrid, nslope
-
-! Local variables
-!----------------
-
-! Code
-!-----
-if (.not. allocated(perice)) allocate(perice(ngrid,nslope))
-
-END SUBROUTINE ini_ice
-!=======================================================================
-
-SUBROUTINE end_ice()
-
-implicit none
-
-! Arguments
-!----------
-
-! Local variables
-!----------------
-
-! Code
-!-----
-if (allocated(perice)) deallocate(perice)
-
-END SUBROUTINE end_ice
-
 !=======================================================================
 
@@ -186 +93 @@
 SUBROUTINE compute_hmaxglaciers(ngrid,nslope,iflat,def_slope_mean,Tice,name_ice,hmax)
 
-use ice_table_mod,  only: rho_ice
-use abort_pem_mod,  only: abort_pem
+use ice_table,  only: rho_ice
+use stop_pem,  only: stop_clean
 use phys_constants, only: pi, g
 
@@ -222 +129 @@
         tau_d = 5.e3
     case default
-        call abort_pem("compute_hmaxglaciers","Type of ice not known!",1)
+        call stop_clean("compute_hmaxglaciers","Type of ice not known!",1)
 end select
 
@@ -249 +156 @@
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 
-use ice_table_mod,  only: rho_ice
-use abort_pem_mod,  only: abort_pem
+use ice_table,      only: rho_ice
+use stop_pem,       only: stop_clean
 use phys_constants, only: pi
 
@@ -345 +252 @@
 END SUBROUTINE computeTcondCO2
 
-END MODULE glaciers_mod
+END MODULE glaciers