Changeset 3568

Timestamp:

Jan 6, 2025, 5:13:01 PM (5 weeks ago)

Author:

llange

Message:

Mars PCM
Fixing bug in soilwater: the specific surface area was missing in the
expression of the adsorption rate when choice_ads = 2.
LL

Location:

trunk/LMDZ.MARS

Files:

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

trunk/LMDZ.MARS/changelog.txt

@@ -4745 +4745 @@
 == 12/12/2024 == JBC
 In the 1D model, merging under the 'paleomars' flag of the modifications for orbital parameters taken from "callphys.def".
+
+== 06/01/2025 == LL
+Fixing bug in soilwater: the specific surface area was missing in the
+expression of the adsorption rate when choice_ads = 2.

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

@@ -681 +681 @@
                         ! first segment of the bilinear function (before monolayer saturation)
                         ! calculate the equilibrium adsorption coefficient
-                        k_ads_eq(ig, ik) = 8.314D0 * ptsoil(ig,ik)/18.D-3 *Kref * dexp(DeltaQ_ads / 8.314D0 *  &
+                        k_ads_eq(ig, ik) = 8.314D0 * rho_soil(ig, ik) * S * ptsoil(ig,ik)/18.D-3 *Kref * dexp(DeltaQ_ads / 8.314D0 *  &
                               (1.D0 / ptsoil(ig, ik) - 1.D0 / Tref)) !  Changed enthalpy_ads to DeltaQ_ads to ensure correct dependance/behaviour of k_ads_eq with temperature
 
@@ -692 +692 @@
                         ! second segment of the bilinear function (after monolayer saturation) ! added 2020
                         ! calculate the equilibrium adsorption coefficient
-                        k_ads_eq2(ig, ik) = 8.314D0 * ptsoil(ig,ik)/18.D-3 * Kref2 * dexp(DeltaQ_ads2 / 8.314D0 *  &
+                        k_ads_eq2(ig, ik) = 8.314D0 * rho_soil(ig, ik) * S * ptsoil(ig,ik)/18.D-3 * Kref2 * dexp(DeltaQ_ads2 / 8.314D0 *  &
                               (1.D0 / ptsoil(ig, ik) - 1.D0 / Tref)) ! Changed enthalpy_ads2 to DeltaQ_ads2 to ensure correct dependance/behaviour of k_ads_eq with temperature