Changeset 3568 for trunk/LMDZ.MARS/libf/phymars

Timestamp:

Jan 6, 2025, 5:13:01 PM (3 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

File:

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

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