Changeset 3274 for trunk/LMDZ.MARS

Timestamp:

Mar 20, 2024, 9:50:21 AM (8 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

Location:

trunk/LMDZ.MARS

Files:

• changelog.txt (modified) (1 diff)
• libf/phymars/soilwater.F90 (modified) (3 diffs)

trunk/LMDZ.MARS/changelog.txt

@@ -4576 +4576 @@
 Small fixes in vdifc for the interaction between the subsurface and frost (the flux from the surface frost to the subsurface was oo large if one removes all the frost at the surface)
 
+== 20/03/2024 == LL
+Small update of soil properties (porosity) when running soilwater to simulate the evolution of massive water ice

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 :