Simulating flood wave propagation in open channels requires solving the non-linear Saint-Venant equations numerically. The Preissmann Implicit Finite-Difference Scheme is widely used for 1D channel routing due to its unconditional numerical stability. In the Preissmann scheme, dependent variables (f) and their spatial/temporal derivatives are discretized on a four-point computational grid cell (i, i+1) between time steps (n, n+1): $f(x,t) \approx \theta \cdot \frac{f_{i}^{n+1} + f_{i+1}^{n+1}}{2} + (1-\theta) \cdot \frac{f_{i}^n + f_{i+1}^n}{2}$ $\frac{\partial f}{\partial x} \approx \theta \cdot \frac{f_{i+1}^{n+1} - f_{i}^{n+1}}{\Delta x} + (1-\theta) \cdot \frac{f_{i+1}^n - f_{i}^n}{\Delta x}$ $\frac{\partial f}{\partial t} \approx \frac{f_{i}^{n+1} + f_{i+1}^{n+1} - f_{i}^n - f_{i+1}^n}{2 \cdot \Delta t}$ Where $\theta$ is a weighting factor $(0.5 \le \theta \le 1.0).$ Setting $\theta \ge 0.55$ ensures numerical damping of high-frequency oscillations. Urban flooding in majo...