Changeset 3248 for trunk

Timestamp:

Mar 1, 2024, 3:58:11 PM (9 months ago)

Author:

llange

Message:

Mars PCM
Following -r 3223, re-introduce the possibility to run the subsurface water ice evolution without adsorption
LL

Location:

trunk/LMDZ.MARS

Files:

• changelog.txt (modified) (1 diff)
• libf/phymars/soilwater.F90 (modified) (2 diffs)
• libf/phymars/vdifc_mod.F (modified) (1 diff)

trunk/LMDZ.MARS/changelog.txt

@@ -4531 +4531 @@
 Fixing bugs in soilwater (array size and indexes in loop)
 
+== 01/03/2024 == LL
+Following -r 3223, re-introduce the possibility to run the subsurface water ice
+evolution without  adsorption
+

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

@@ -671 +671 @@
                         ! calculate the absorption coefficient
                         Ka2(ig, ik) = kinetic_factor * k_ads_eq2(ig, ik) / (1.D0 + k_ads_eq2(ig, ik) / porosity(ig, ik)) 
+                  
+                  elseif (choice_ads.eq.0)  then!No adsorption/desorption
+                  
+                        Ka(ig, ik) = 0.
+                        Kd(ig, ik) = 0.
+                        Ka2(ig, ik) = 0.
+                        Kd2(ig, ik) = 0.
                   endif
                   
                   ! calculate the amount of water vapor at monolayer saturation ! modified 2020
-                  nsat(ig, ik) = adswater_sat_mono * Kd(ig, ik) / Ka(ig, ik) ! 
+                  if(choice_ads.ne.0) nsat(ig, ik) = adswater_sat_mono * Kd(ig, ik) / Ka(ig, ik)  
 
                   ! calculate the intercept of the second line in the bilinear function ! added 2020; corrected 2024 
@@ -1455 +1462 @@
                   
                   ! 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.0) preduite(ig, ik) = znsoil(ig, ik) / nsat(ig, ik)
                   
                   ! write the results to the return variable

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

@@ -1100 +1100 @@
                      endif  ! of "if exchange = true"
                  endif  ! of "if not.watercaptag"
-              endif ! adsorption
+              endif ! adsorption          
 
              if(.not.watercaptag(ig).and.(.not.adsorption_soil)) then