Unsteady non-uniform open channel flows, such as surge waves created by sudden sluice gate movements, are modeled by transforming the 1D hyperbolic Saint-Venant partial differential equations into ordinary differential equations using the Method of Characteristics (MOC). The dynamic flow field yields two characteristic velocity paths $(C^+ and C^-):$ $\frac{dx}{dt} = v \pm c = v \pm \sqrt{g \cdot y}$ Where $v$ is mean flow velocity, $y$ is flow depth, and $c$ is shallow-water wave celerity. Along these characteristic lines, the Riemann invariants $(I_+ and I_-)$ remain constant in frictionless rectangular channels: $I_+ = v + 2 \cdot \sqrt{g \cdot y} = \text{constant along } C^+$ $I_- = v - 2 \cdot \sqrt{g \cdot y} = \text{constant along } C^-$ Evaluating positive surge waves (dam breaks or sudden canal gate closures) along major headworks across India requires accurate tracking of wave front arrival times to prevent embankment overtopping. Modern canal automation systems comb...