Changeset 3274 for trunk/LMDZ.MARS/libf

Timestamp:

Mar 20, 2024, 9:50:21 AM (10 months ago)

Author:

llange

Message:

Mars PCM
Small update of soil properties (porosity) when running soilwater to simulate the evolution of massive water ice
LL

File:

• trunk/LMDZ.MARS/libf/phymars/soilwater.F90 (modified) (3 diffs)

trunk/LMDZ.MARS/libf/phymars/soilwater.F90

@@ -267 +267 @@
 ! variables for the 1D case
 real*8, parameter :: mmr_h2o = 0.6D-4     ! Water vapor mass mixing ratio (for initialization only)
-
 real*8 :: diff(ngrid, nsoil)              ! difference between total water content and ice + vapor content 
                                           ! (only used for output, inconstistent: should just be adswater)
-real :: var_flux_soil(ngrid, nsoil)
-
-real :: default_diffcoeff = 4e-4
-
+real :: var_flux_soil(ngrid, nsoil)       ! Output of the flux between soil layers (kg/m^3/s)
+real :: default_diffcoeff = 4e-4          ! Diffusion coefficient by default (no variations with Temperature and pressure (m/s^2)
+real :: tol_massiveice = 50.              ! Tolerance to account for massive ice (kg/m^3)
 ! Loop variables and counters
 integer :: ig, ik, i, j                   ! loop variables
@@ -350 +348 @@
                         ! These properties are defined here in order to enable custom profiles
                         if(ads_massive_ice) then
-                           porosity_ice_free(ig, ik) = 0.9999
+                              if(pqsoil(ig, ik, igcm_h2o_ice_soil).gt.tol_massiveice) then
+                                    porosity_ice_free(ig, ik) = 0.999999
+                              else
+                                    porosity_ice_free(ig, ik) = porosity_reg
+                              endif
                         else
                            porosity_ice_free(ig, ik) = porosity_reg
@@ -521 +523 @@
 endif  ! of "if firstcall_soil"
 
+
+
+
+! Update properties in case of the sublimation of massive ice
+
+if(ads_massive_ice) then
+      do ig = 1, ngrid      
+            do ik = 1, nsoil  
+                  if(pqsoil(ig, ik, igcm_h2o_ice_soil).gt.tol_massiveice) then
+                        porosity_ice_free(ig, ik) = 0.999999
+                  else
+                        porosity_ice_free(ig, ik) = porosity_reg
+                  endif
+            enddo
+      enddo
+endif     
 
 ! Computation of new (new step) transport coefficients :