In this study, the numerical model developed by Hong et al.(2008) was improved to be applied to rapidly varying flows such as the inundation of dry land or flow transitions due to large gradients of the bathymetry. A numerical approximation was applied that was consistent with the conservation of momentum in flow expansions and with the Bernoulli equation in flow contractions. The approximation was second order, but the accuracy reduced to first order near extreme values by the use of a minmod limiter. The modified model was verified by acomparison with the theoretical critical depth of weir, and for sufficiently smooth conditions and a fine grid size, both approximations converged to the same solution. In terms of the grid size, it was more effective at obtaining solutions than the previous model and reproduced the inundation of dry land.