Changeset 3126 for trunk/LMDZ.MARS/libf/phymars/soilwater.F90

Timestamp:

Nov 15, 2023, 3:13:25 PM (13 months ago)

Author:

llange

Message:

Mars PCM

• Update in soilwater: adding the possibility to run without adsorption, but with the possibility to run with seasonal frost forming in the subsurface
• THe choice of the isotherm for adsorption can be now done by setting the integer choice_ads in the callphys.def choice_ads = 1 adsorption rate is computed with the H2O thermal speed; choice_ads = 2 adsorption rate is computed based on exeperimental resutls, choice_ads =3 no adsorption

LL

File:

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

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

@@ -5 +5 @@
 
 
-      use comsoil_h, only: igcm_h2o_vap_soil, igcm_h2o_ice_soil, igcm_h2o_vap_ads, layer, mlayer
+      use comsoil_h, only: igcm_h2o_vap_soil, igcm_h2o_ice_soil, igcm_h2o_vap_ads, layer, mlayer, choice_ads
       use comcstfi_h
       use tracer_mod
@@ -204 +204 @@
 real*8, parameter:: Sm = 10.6D-20         ! Surface of the water molecule (m2) (only needed in the theoretical formula which is not used right now)
 
-integer, parameter :: choice_ads = 1      ! Choice of adsorption - desorption constants 3: no adsorption
 
 ! Reference values for choice_ads = 2
@@ -350 +349 @@
                               saturation_water_ice(ig, ik) = min(ice(ig, ik) / (rho_H2O_ice * porosity_ice_free(ig, ik)), 0.999D-0)
                               porosity(ig, ik) = porosity_ice_free(ig, ik) * (1.D0 - saturation_water_ice(ig, ik))
-                              
+
                               if (choice_ads.eq.1) then
                                     vth(ig, ik) = dsqrt(8.D0 * 8.314D0 * tsoil(ig, nsoil - 4) &
@@ -373 +372 @@
                                           
                                     Ka(ig, ik) = kinetic_factor * k_ads_eq(ig, ik) / (1.D0 + k_ads_eq(ig, ik) / porosity(ig, ik))
+        
                               endif
 
-                              adswater(ig, ik) = min(Ka(ig, ik) / Kd(ig, ik) * znsoil(ig, ik), adswater_sat)
+                              if(choice_ads .ne. 3) adswater(ig, ik) = min(Ka(ig, ik) / Kd(ig, ik) * znsoil(ig, ik), adswater_sat)
                               
                         else  ! in 3D simulations initialisation happens with newstart.F
@@ -381 +381 @@
                               ice(ig, ik) = pqsoil(ig, ik, igcm_h2o_ice_soil)
                               adswater(ig, ik) = pqsoil(ig, ik, igcm_h2o_vap_ads)
+                        endif
+
+                        if (choice_ads.eq.3) then ! no adsorption
+
+                               Ka(:, :) = 0.
+                               Kd(:, :) = 0.
+                               adswater(:,:) = 0.
+
                         endif
 
@@ -584 +592 @@
                         ! calculate the absorption coefficient
                         Ka(ig, ik) = kinetic_factor * k_ads_eq(ig, ik) / (1.D0 + k_ads_eq(ig, ik) / porosity(ig, ik))
+                  else  ! no ads
+
+                        Kd(ig, ik) = 0.
+
+                        Ka(ig, ik) = 0.
                   endif
                   
                   ! calculate the amount of water vapor at adorption saturation
-                  nsat(ig, ik) = adswater_sat * Kd(ig, ik) / Ka(ig, ik)
+                 
+                  if (choice_ads.ne.3) nsat(ig, ik) = adswater_sat * Kd(ig, ik) / Ka(ig, ik)
 
                   ! calculate C, E, and F coefficients for later calculations
@@ -1285 +1299 @@
                   
                   ! calculate how close the water vapor content is to saturizing the adsorbed water
-                  preduite(ig, ik) = znsoil(ig, ik) / nsat(ig, ik)
+                  if (choice_ads.ne.3) preduite(ig, ik) = znsoil(ig, ik) / nsat(ig, ik)
                   
                   ! write the results to the return variable