# Pump and system operating point
# Scope: Quadratic pump-head/system-head intersection; coefficients supplied.
# Status: illustrative worked example; engineering review pending.
# Inputs are example values. Checks are editable example requirements, not standards.

shutoff = 20 m
staticHead = 5 m
pumpCoefficient = 100000000 s^2/m^5
systemCoefficient = 50000000 s^2/m^5

flow := sqrt((shutoff-staticHead)/(pumpCoefficient+systemCoefficient)) [m^3/s]
head := staticHead+systemCoefficient*flow^2 [m]

check positive available head: shutoff > staticHead

scenario alternative: shutoff = 16 m