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

Bug fixed in oxygen-saturated methane oxidation

parent c28a8cf5
No related branches found
No related tags found
No related merge requests found
......@@ -1152,9 +1152,9 @@ endif ifdeep
case(2)
! Michaelis-Menten with methane concentrations only (oxygen-independent or oxygen-saturated)
! (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
qwater(i) = 0.5*( (qwater(i) + Vmax*dt - k_ch40) + sqrt(x) )
qwater(i) = 0.5*( (qwater(i) - Vmax*dt - k_ch40) + sqrt(x) )
x = dt*Vmax*qwater(i)/(k_ch40 + qwater(i))
oxyg(i) = oxyg(i) - 2.*x
DIC(i) = DIC(i) + x
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment