From 4279e0ac36f622232b2bc93215308f321c43cca4 Mon Sep 17 00:00:00 2001 From: Victor Stepanenko <stepanen@srcc.msu.ru> Date: Mon, 22 Feb 2021 19:58:01 +0300 Subject: [PATCH] Bug fixed in methane oxidation --- source/model/methane_mod.f90 | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/source/model/methane_mod.f90 b/source/model/methane_mod.f90 index 660f5b2..5d20cb9 100644 --- a/source/model/methane_mod.f90 +++ b/source/model/methane_mod.f90 @@ -1154,8 +1154,8 @@ endif ifdeep ! (implicit scheme) x = (k_ch40 + qwater(i))*(k_ch40 + qwater(i)) + & & 2.*(k_ch40 - qwater(i))*Vmaxw*dt + Vmaxw*dt*Vmaxw*dt - qwater(i) = 0.5*( (qwater(i) - Vmax*dt - k_ch40) + sqrt(x) ) - x = dt*Vmax*qwater(i)/(k_ch40 + qwater(i)) + qwater(i) = 0.5*( (qwater(i) - Vmaxw*dt - k_ch40) + sqrt(x) ) + x = dt*Vmaxw*qwater(i)/(k_ch40 + qwater(i)) oxyg(i) = oxyg(i) - 2.*x DIC(i) = DIC(i) + x case(3) -- GitLab