Skip to content
Snippets Groups Projects
Commit 4279e0ac authored by Victor Stepanenko's avatar Victor Stepanenko
Browse files

Bug fixed in methane oxidation

parent 02c313d2
No related branches found
No related tags found
No related merge requests found
Show whitespace changes
Inline Side-by-side
source/model/methane_mod.f90
+ 2
− 2
View file @ 4279e0ac
...@@ -1154,8 +1154,8 @@ endif ifdeep ...@@ -1154,8 +1154,8 @@ endif ifdeep
! (implicit scheme) ! (implicit scheme)
x = (k_ch40 + qwater(i))*(k_ch40 + qwater(i)) + & x = (k_ch40 + qwater(i))*(k_ch40 + qwater(i)) + &
& 2.*(k_ch40 - qwater(i))*Vmaxw*dt + Vmaxw*dt*Vmaxw*dt & 2.*(k_ch40 - qwater(i))*Vmaxw*dt + Vmaxw*dt*Vmaxw*dt
qwater(i) = 0.5*( (qwater(i) - Vmax*dt - k_ch40) + sqrt(x) ) qwater(i) = 0.5*( (qwater(i) - Vmaxw*dt - k_ch40) + sqrt(x) )
x = dt*Vmax*qwater(i)/(k_ch40 + qwater(i)) x = dt*Vmaxw*qwater(i)/(k_ch40 + qwater(i))
oxyg(i) = oxyg(i) - 2.*x oxyg(i) = oxyg(i) - 2.*x
DIC(i) = DIC(i) + x DIC(i) = DIC(i) + x
case(3) case(3)
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment