Changeset 4136 for trunk/LMDZ.COMMON/libf/evolution/soil.F90

Timestamp:

Mar 17, 2026, 11:12:37 AM (3 weeks ago)

Author:

jbclement

Message:

PEM:

• Temporary memory load reduced in "xios_data".
• Rename local variables in procedures to avoid masking variables in parent scope.

JBC

File:

• trunk/LMDZ.COMMON/libf/evolution/soil.F90 (modified) (4 diffs)

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

@@ -25 +25 @@
 ! PARAMETERS
 ! ----------
-! Density and porosity of regolith and breccia
-real(dp), parameter :: rho_regolith      = 2000._dp ! [kg/m^3], Zent et al. 1995; Buhler and Piqueux 2021
+real(dp), parameter :: rho_regolith      = 2000._dp ! Density of the martian regolith [kg/m^3], Zent et al. 1995; Buhler and Piqueux 2021
 real(dp), parameter :: regolith_porosity = 0.4_dp   ! Porosity of the martian regolith for a random loose packing of monodisperse sphere, Scott 1960
 real(dp), parameter :: breccia_porosity  = 0.1_dp   ! Porosity of breccia
@@ -255 +254 @@
 
 !=======================================================================
-SUBROUTINE set_soil(TI)
+SUBROUTINE set_soil(therm_initia)
 !-----------------------------------------------------------------------
 ! NAME
@@ -291 +290 @@
 ! ARGUMENTS
 ! ---------
-real(dp), dimension(:,:,:), intent(inout) :: TI ! Thermal inertia in the PEM [SI]
+real(dp), dimension(:,:,:), intent(inout) :: therm_initia ! Thermal inertia in the PEM [SI]
 
 ! LOCAL VARIABLES
@@ -334 +333 @@
     do islope = 1,nslope
         do iloop = 1,nsoil_PCM
-            TI(ig,iloop,islope) = TI_PCM(ig,iloop,islope)
+            therm_initia(ig,iloop,islope) = TI_PCM(ig,iloop,islope)
         end do
         if (nsoil > nsoil_PCM) then
             do iloop = nsoil_PCM + 1,nsoil
-                TI(ig,iloop,islope) = TI_PCM(ig,nsoil_PCM,islope)
+                therm_initia(ig,iloop,islope) = TI_PCM(ig,nsoil_PCM,islope)
             end do
         end if