Calc.Greene.Qe: Calculate the Greene Qe value using numerical methods
In seanpor/phoslam: Helper Functions for Phosphate Load Apportionment

Calculate the Greene Qe value using numerical methods

1	Calc.Greene.Qe(Greene, adf, sDEBUG = NULL)

`Greene`	a Greene object XXX
`adf`	a dataframe containing XXX
`sDEBUG`	a parameter not normally used except in debugging, set to be an integer >5 for maximum debug information printed.

Here Qe has to be calculated numerically - so the flow value at which Point is bigger than (Diffuse + Hysteresis)

So to figure out Qe, need to find the first zero point of d1$GreeneDH.L-d1$GreenePoint.L in the Q range that is actually there, i.e. Qe > 0 and Qe<max(Q) but we can get this from the function itself... by substituting the coefficients we have into the function itself... i.e. we want to find the zero for the function (B*Q + C*Q^2) - A/Q in the range 0, to max(d1$Q) nb. should limit Q to be in the range 0, Q value of min(d1$GreenePoint.L-d1$GreeneDH.L)

seanpor/phoslam documentation built on May 16, 2019, 9:11 p.m.