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

Timestamp:

Feb 5, 2025, 3:06:50 PM (9 hours ago)

Author:

jbclement

Message:

PEM:

• The sublimation flux can now be modified by the new function 'recomp_tend_h2o' to account for the growth of a dust lag layer (see Eran Vos's note for the formula).
• Addition of a function 'itp_tsoil' to compute the soil temperature at any point in the soil profile.

JBC

File:

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

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

@@ -65 +65 @@
 use pemetat0_mod,               only: pemetat0
 use read_data_PCM_mod,          only: read_data_PCM
-use recomp_tend_co2_mod,        only: recomp_tend_co2
+use recomp_tend_mod,            only: recomp_tend_co2, recomp_tend_h2o
 use compute_soiltemp_mod,       only: compute_tsoil_pem, shift_tsoil2surf
 use writediagpem_mod,           only: writediagpem, writediagsoilpem
@@ -211 +211 @@
 real, dimension(:,:,:,:), allocatable :: tsoil_timeseries                 ! Grid points x Soil x Slope x Times field: Soil temperature timeseries [K]
 real, dimension(:,:,:,:), allocatable :: tsoil_PEM_timeseries             ! Grid points x Soil x Slope x Times field: Soil temperature timeseries for PEM [K]
+real, dimension(:,:,:,:), allocatable :: tsoil_PEM_timeseries_old         ! Grid points x Soil x Slope x Times field: Soil temperature timeseries for PEM at the previous time step [K]
 real, dimension(:,:,:,:), allocatable :: watersoil_density_timeseries     ! Grid points x Soil x Slope x Times Water soil density timeseries [kg /m^3]
 real, dimension(:,:),     allocatable :: watersurf_density_avg            ! Grid points x Slope: Averaged  water surface density [kg/m^3]
@@ -227 +228 @@
 real, dimension(:,:),     allocatable :: zshift_surf                      ! Elevation shift for the surface [m]
 real, dimension(:,:),     allocatable :: zlag                             ! Newly built lag thickness [m]
+real, dimension(:,:),     allocatable :: icetable_depth_old               ! Old depth of the ice table
 
 ! Some variables for the PEM run
@@ -857 +859 @@
 ! II_d.3 Update soil temperature
         write(*,*)"> Updating soil temperature profile"
+        allocate(tsoil_avg_old(ngrid,nsoilmx_PEM),tsoil_PEM_timeseries_old(ngrid,nsoilmx_PEM,nslope,timelen))
+        tsoil_PEM_timeseries_old = tsoil_PEM_timeseries
         do islope = 1,nslope
             tsoil_avg_old = tsoil_PEM(:,:,islope)
@@ -878 +882 @@
 
 ! II_d.4 Update the ice table
-        allocate(icetable_thickness_old(ngrid,nslope),ice_porefilling_old(ngrid,nsoilmx_PEM,nslope))
+        allocate(icetable_thickness_old(ngrid,nslope),ice_porefilling_old(ngrid,nsoilmx_PEM,nslope),icetable_depth_old(ngrid,nslope))
         if (icetable_equilibrium) then
             write(*,*) "> Updating ice table (equilibrium method)"
@@ -887 +891 @@
             write(*,*) "> Updating ice table (dynamic method)"
             ice_porefilling_old = ice_porefilling
+            icetable_depth_old = icetable_depth
             allocate(porefill(nsoilmx_PEM))
             do ig = 1,ngrid
@@ -971 +976 @@
 !------------------------
     call recomp_tend_co2(ngrid,nslope,timelen,d_co2ice,d_co2ice_ini,co2_ice,emis,vmr_co2_PCM,vmr_co2_PEM_phys,ps_timeseries,ps_avg_global_old,ps_avg_global_new)
-    deallocate(vmr_co2_PEM_phys)
+    do ig = 1,ngrid
+        do islope = 1,nslope
+            call recomp_tend_h2o(icetable_depth_old(ig,islope),icetable_depth(ig,islope),tsurf_avg(ig,islope),tsoil_PEM_timeseries_old(ig,:,islope,:),tsoil_PEM_timeseries(ig,:,islope,:),d_h2oice(ig,islope))
+        enddo
+    enddo
+    deallocate(vmr_co2_PEM_phys,icetable_depth_old,tsoil_PEM_timeseries_old)
 
 !------------------------