Changeset 3272 for trunk/LMDZ.MARS

Timestamp:

Mar 19, 2024, 10:52:50 AM (8 months ago)

Author:

llange

Message:

Mars PCM
Adapt soil water to run with massive ice and not pore filling ice + constant diffusion coefficient
LL

Location:

trunk/LMDZ.MARS

Files:

• changelog.txt (modified) (2 diffs)
• libf/phymars/comsoil_h.F90 (modified) (2 diffs)
• libf/phymars/conf_phys.F (modified) (2 diffs)
• libf/phymars/dyn1d/init_testphys1d_mod.F90 (modified) (3 diffs)
• libf/phymars/soilwater.F90 (modified) (5 diffs)

trunk/LMDZ.MARS/changelog.txt

@@ -4549 +4549 @@
 Reversing -r 3250 (the ""bug"" I fixed was not a bug but my mystake)
 
-== 12/03/2024 ==Jliu
+== 12/03/2024 == Jliu
 Synchronize all default non-orographic gravity wave parameters in callphys.def with the ones in scheme. The parameters are as follows:
 #Eliassen-Palm FLux(only if calllott_nonoro=.true.)
@@ -4572 +4572 @@
 nonoro_gwimixing_vdl=1.5
 
-
+== 19/03/2024 == LL
+Adapt soil water to run with massive ice and not pore filling ice + constant diffusion coefficient

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

@@ -32 +32 @@
 logical, save                                  :: adsorption_soil      ! boolean to call adosrption (or not)
 real,    save                                  :: choice_ads           ! Choice for adsorption isotherm (3 means no adsorption, see soilwater.F90)
+logical, save                                  :: ads_const_D          ! boolean to have constant diffusion coefficient
+logical, save                                  :: ads_massive_ice      ! boolean to have massive subsurface ice
 real,    save, allocatable, dimension(:,:,:,:) :: qsoil                ! subsurface tracers (kg/m^3 of regol)
 integer, parameter :: nqsoil = 3                                       ! number of subsurface tracers, only three when working with water
@@ -39 +41 @@
 REAL, parameter :: porosity_reg  = 0.45
 !$OMP THREADPRIVATE(adsorption_soil,qsoil,choice_ads)
+!$OMP THREADPRIVATE(ads_const_D,ads_massive_ice)
 
 !=======================================================================

trunk/LMDZ.MARS/libf/phymars/conf_phys.F

@@ -40 +40 @@
      &                           lag_layer 
       use microphys_h, only: mteta
-      use comsoil_h, only: adsorption_soil, choice_ads
+      use comsoil_h, only: adsorption_soil, choice_ads,ads_const_D,
+     &                     ads_massive_ice
       use nonoro_gwd_mix_mod, only: calljliu_gwimix
 
@@ -1041 +1042 @@
          if (adsorption_soil) choice_ads = 1
          call getin_p("choice_ads",choice_ads) 
-         
+
+         ads_const_D = .false.
+         call getin_p("ads_const_D",ads_const_D)
+
+         ads_massive_ice = .false.
+         call getin_p("ads_massive_ice",ads_massive_ice)
+
          poreice_tifeedback = .false.
          call getin_p("poreice_tifeedback",poreice_tifeedback)

trunk/LMDZ.MARS/libf/phymars/dyn1d/init_testphys1d_mod.F90

@@ -22 +22 @@
 use comsoil_h,                only: volcapa, nsoilmx, inertiesoil, inertiedat, layer, mlayer, flux_geo, tsoil, &
                                     adsorption_soil, qsoil, igcm_h2o_ice_soil, porosity_reg, &
-                                    ini_comsoil_h_slope_var, end_comsoil_h_slope_var
+                                    ini_comsoil_h_slope_var, end_comsoil_h_slope_var, ads_massive_ice
 use comvert_mod,              only: ap, bp, aps, bps, pa, preff, presnivs, pseudoalt, scaleheight
 use dimradmars_mod,           only: tauvis, totcloudfrac, albedo, ini_dimradmars_mod_slope_var, end_dimradmars_mod_slope_var
@@ -688 +688 @@
             if(adsorption_soil) then
             ! add subsurface ice in qsoil if one runs with adsorption
-                qsoil(1,1,igcm_h2o_ice_soil,:) = (ice_depth-layer(2))/(layer(1) - layer(2))*porosity_reg*rho_H2O_ice
-                qsoil(1,2:,igcm_h2o_ice_soil,:) = porosity_reg*rho_H2O_ice
+               if(ads_massive_ice) then
+                   qsoil(1,1,igcm_h2o_ice_soil,:) = (ice_depth-layer(2))/(layer(1) - layer(2))*rho_H2O_ice ! assume no porosity in the ice
+                   qsoil(1,2:,igcm_h2o_ice_soil,:) = rho_H2O_ice
+               else
+                   qsoil(1,1,igcm_h2o_ice_soil,:) = (ice_depth-layer(2))/(layer(1) - layer(2))*porosity_reg*rho_H2O_ice
+                   qsoil(1,2:,igcm_h2o_ice_soil,:) = porosity_reg*rho_H2O_ice
+               endif            
             endif
         else ! searching for the ice/regolith boundary:
@@ -713 +718 @@
                 else
                    qsoil(1,iref+1:,igcm_h2o_ice_soil,:) = porosity_reg*rho_H2O_ice
-                   qsoil(1,iref,igcm_h2o_ice_soil,:) = (ice_depth-layer(iref))/(layer(iref-1) - layer(iref))*porosity_reg*rho_H2O_ice                
+                   qsoil(1,iref,igcm_h2o_ice_soil,:) = (ice_depth-layer(iref))/(layer(iref-1) - layer(iref))*porosity_reg*rho_H2O_ice
                 endif
             endif

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

@@ -4 +4 @@
 
 
-      use comsoil_h, only: igcm_h2o_vap_soil, igcm_h2o_ice_soil, igcm_h2o_vap_ads, layer, mlayer, choice_ads, porosity_reg
+      use comsoil_h, only: igcm_h2o_vap_soil, igcm_h2o_ice_soil, igcm_h2o_vap_ads, layer, mlayer, choice_ads, porosity_reg, &
+                           ads_const_D, ads_massive_ice 
       use comcstfi_h
       use tracer_mod
@@ -271 +272 @@
 real :: var_flux_soil(ngrid, nsoil)
 
+real :: default_diffcoeff = 4e-4
+
 ! Loop variables and counters
 integer :: ig, ik, i, j                   ! loop variables
@@ -346 +349 @@
                         
                         ! These properties are defined here in order to enable custom profiles
-                        porosity_ice_free(ig, ik) = porosity_reg
+                        if(ads_massive_ice) then
+                           porosity_ice_free(ig, ik) = 0.9999
+                        else
+                           porosity_ice_free(ig, ik) = porosity_reg
+                        endif
                         tortuosity(ig, ik) = 1.5D0
                         rho_soil(ig, ik) = 1.3D3 ! in kg/m3 of regolith (incl. porosity)
@@ -569 +576 @@
                   D_mid(ig, ik) = D0(ig, ik) * (1.D0 - saturation_water_ice(ig, ik)) ** 2.D0 * 1.D0 / (1.D0 / Dm(ig, ik) + 1.D0 / Dk(ig, ik))
 
+                  if(ads_const_D) D_mid(ig, ik) = default_diffcoeff
+
                   ! Midlayer diffusion coeff. (correlation by Rogers and Nielson, 1991)
                   !                             D_mid(ig, ik) = D0(ig, ik) * (1. - saturation_water_ice(ig, ik)) * exp(-6. * saturation_water_ice(ig, ik)  &
@@ -595 +604 @@
                   
                   D_inter(ig, ik) = D0(ig, ik) * (1.D0 - saturation_water_ice_inter(ig, ik)) ** 2.D0 * 1.D0 / (1.D0 / Dm_inter(ig, ik) + 1.D0 / Dk_inter(ig, ik))
-                  
+                 
+                  if(ads_const_D) D_inter(ig, ik) = default_diffcoeff
+ 
 !                  D_inter(ig, ik) = ( (layer(ik) - mlayer(ik - 1)) * D_mid(ig, ik + 1) &  ! old implementation with direct interpolation
 !                        + (mlayer(ik) - layer(ik)) * D_mid(ig, ik) ) / (mlayer(ik) - mlayer(ik - 1))