check if k needs to be changed based on current Q vector
update k
recalculate muskingum coefficients
recalculate the LHS matrix
calculate the iteration's rhs or b vector (existing)
solve (existing)
Need to determine the data structure for determining k or x values based on q thresholds
equations? require recomputing every time. Or only if q is >5% different than the previous Q or Q the last time it was calculate?
tables? what are the thresholds for changing k? can they be known a priori or must you iteratively solve to reach an equilibrium state when the thresholds stop changing (e.g. return periods)?
Allow routing with variable k (and/or x) values
On each iteration
Need to determine the data structure for determining k or x values based on q thresholds