Changeset 3126 for trunk/LMDZ.MARS/libf

Timestamp:

Nov 15, 2023, 3:13:25 PM (16 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

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• comsoil_h.F90 (modified) (2 diffs)
• conf_phys.F (modified) (2 diffs)
• soil_settings.F (modified) (1 diff)
• soilwater.F90 (modified) (7 diffs)

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

@@ -32 +32 @@
 ! Subsurface tracers:
   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)
   integer, parameter :: nqsoil = 3            ! number of subsurface tracers, only three when working with water
   real,save,allocatable :: qsoil(:,:,:,:)     ! subsurface tracers (kg/m^3 of regol)
@@ -37 +38 @@
   integer, parameter :: igcm_h2o_ice_soil = 2
   integer, parameter :: igcm_h2o_vap_ads  = 3
-!$OMP THREADPRIVATE(adsorption_soil,qsoil)
+!$OMP THREADPRIVATE(adsorption_soil,qsoil,choice_ads)
 
 

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

@@ -40 +40 @@
      &                           lag_layer 
       use microphys_h, only: mteta
-      use comsoil_h, only: adsorption_soil
+      use comsoil_h, only: adsorption_soil, choice_ads
 
       IMPLICIT NONE
@@ -1017 +1017 @@
      &        "Adsorption must be used with water = true",1)
          endif  
+
+         if(adsorption_soil) choice_ads = 1
+         call getin_p("choice_ads",choice_ads) 
+           
 c        ----------------------------------------------------------
 

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

@@ -450 +450 @@
              call get_field(txt,qsoil(:,:,igcm_h2o_vap_soil,:),found,
      &                                                 indextime)
-             write(*,*) 'found',found
-             write(*,*) 'igcm_',igcm_h2o_vap_soil
-             write(*,*) 'q=',qsoil(:,:,:,:)
              if (.not.found) then
                write(*,*) "phyetat0: Failed loading <",trim(txt),">"

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