Changeset 3673 for trunk/LMDZ.COMMON/libf/evolution/pemetat0.F90

Timestamp:

Mar 7, 2025, 10:01:32 AM (4 months ago)

Author:

jbclement

Message:

PEM:

• Correction for the update of average pressure at the end.
• Addition of pore-filled ice information for 'stratum' in the layering algorithm.
• Deactivating the soil temperature shifting because of bad values given by interpolation.

JBC

File:

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

trunk/LMDZ.COMMON/libf/evolution/pemetat0.F90 ¶

@@ -65 +65 @@
 real, dimension(ngrid,nsoil_PEM,nslope), intent(inout) :: watersoil_avg       ! surface water ice density, yearly averaged (kg/m^3)
 ! local
-real, dimension(ngrid,nsoil_PEM,nslope) :: tsoil_startpem           ! soil temperature saved in the start [K]
-real, dimension(ngrid,nsoil_PEM,nslope) :: TI_startPEM              ! soil thermal inertia saved in the start [SI]
-logical                                 :: found, found2            ! check if variables are found in the start
-integer                                 :: iloop, ig, islope, isoil ! index for loops
-real                                    :: delta                    ! Depth of the interface regolith-breccia, breccia -bedrock [m]
-character(2)                            :: num                      ! intermediate string to read PEM start sloped variables
-logical                                 :: startpem_file            ! boolean to check if we read the startfile or not
-real, dimension(:,:,:,:), allocatable   :: stratif_array            ! Array for stratification (layerings)
+real, dimension(ngrid,nsoil_PEM,nslope) :: tsoil_startpem    ! soil temperature saved in the start [K]
+real, dimension(ngrid,nsoil_PEM,nslope) :: TI_startPEM       ! soil thermal inertia saved in the start [SI]
+logical                                 :: found, found2     ! check if variables are found in the start
+integer                                 :: iloop, ig, islope ! index for loops
+real                                    :: delta             ! Depth of the interface regolith-breccia, breccia -bedrock [m]
+character(2)                            :: num               ! intermediate string to read PEM start sloped variables
+logical                                 :: startpem_file     ! boolean to check if we read the startfile or not
+real, dimension(:,:,:,:), allocatable   :: stratif_array     ! Array for stratification (layerings)
 
 #ifdef CPP_STD
@@ -143 +143 @@
             nb_str_max = int(inquire_dimension_length('nb_str_max'))
         endif
-        allocate(stratif_array(ngrid,nslope,nb_str_max,6))
+        allocate(stratif_array(ngrid,nslope,nb_str_max,7))
         stratif_array = 0.
         do islope = 1,nslope
@@ -156 +156 @@
             call get_field('stratif_slope'//num//'_dust_volfrac',stratif_array(:,islope,:,5),found2)
             found = found .or. found2
-            call get_field('stratif_slope'//num//'_air_volfrac',stratif_array(:,islope,:,6),found2)
+            call get_field('stratif_slope'//num//'_pore_volfrac',stratif_array(:,islope,:,6),found2)
+            found = found .or. found2
+            call get_field('stratif_slope'//num//'_icepore_volfrac',stratif_array(:,islope,:,7),found2)
             found = found .or. found2
             if (.not. found) then
@@ -267 +269 @@
                 call compute_tsoil_pem(ngrid,nsoil_PEM,.false.,TI_PEM(:,:,islope),timestep,tsurf_avg_yr2(:,islope),tsoil_startpem(:,:,islope))
 
-                do iloop = nsoil_PCM + 1,nsoil_PEM
-                    tsoil_PEM(:,iloop,islope) = tsoil_startpem(:,iloop,islope)
-                enddo
+                tsoil_PEM(:,nsoil_PCM + 1:nsoil_PEM,islope) = tsoil_startpem(:,nsoil_PCM + 1:nsoil_PEM,islope)
             endif !found
 
-            do isoil = nsoil_PCM + 1,nsoil_PEM
-                watersoil_avg(:,isoil,islope) = exp(beta_clap_h2o/tsoil_PEM(:,isoil,islope) + alpha_clap_h2o)/tsoil_PEM(:,isoil,islope)*mmol(igcm_h2o_vap)/(mugaz*r)
-            enddo
+            watersoil_avg(:,nsoil_PCM + 1:nsoil_PEM,islope) = exp(beta_clap_h2o/tsoil_PEM(:,nsoil_PCM + 1:nsoil_PEM,islope) + alpha_clap_h2o)/tsoil_PEM(:,nsoil_PCM + 1:nsoil_PEM,islope)*mmol(igcm_h2o_vap)/(mugaz*r)
         enddo ! islope
         write(*,*) 'PEMETAT0: TSOIL done'
@@ -431 +429 @@
 
 ! First raw initialization
-            do isoil = nsoil_PCM + 1,nsoil_PEM
-                do ig = 1,ngrid
-                    watersoil_avg(ig,isoil,islope) = exp(beta_clap_h2o/tsoil_PEM(ig,isoil,islope) + alpha_clap_h2o)/tsoil_PEM(ig,isoil,islope)*mmol(igcm_h2o_vap)/(mugaz*r)
-                enddo
-            enddo
+            watersoil_avg(:,nsoil_PCM + 1:nsoil_PEM,islope) = exp(beta_clap_h2o/tsoil_PEM(:,nsoil_PCM + 1:nsoil_PEM,islope) + alpha_clap_h2o)/tsoil_PEM(:,nsoil_PCM + 1:nsoil_PEM,islope)*mmol(igcm_h2o_vap)/(mugaz*r)
         enddo !islope
         write(*,*) 'PEMETAT0: TSOIL done'