GreeneQe: Solve the Greene equation for Qe

In seanpor/phoslam: Helper Functions for Phosphate Load Apportionment

Description

Qe is when the larger contribution switches from point to diffuse - there is always point contribution but at high flows it is outweighed by the diffuse contribution. GreeneQe takes a set of three Greene coefficient, A,B,C and solves the Greene equation for the $TRP = A/Q + B*Q + C*Q^2$, i.e. TRP = Point + Diffuse + Hysteresis flow rate Q at which Point = Diffuse + Hysteresis

Usage

GreeneQe(A, B, C, rightedge = 1, sDEBUG = NULL)

Arguments

A	Greene equation coefficient
B	Greene equation coefficient
C	Greene equation coefficient
rightedge	The flow rate parameter which is expected to be greater than Qe
sDEBUG	a parameter not normally used except in debugging, set to be an integer >5 for maximum debug information printed.

Details

Note that sometimes this has problems converging hence there is a further parater right edge which limits the range of Q at the upper end.

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