Two numerical models for 2D flood flows are presented. One model is first-order accurate and another is second-order accurate. Roe's numerical flux is used to develop the first-order accurate model, while second-order accuracy, in space and time, is obtained by using the Lax-Wendroff numerical flux. A simple operator splitting is found to yield the same results as that obtained by using more complicated, and thus, time consuming, operator splitting. Roe's approximate Jacobian is used for conservative properties and Harten and Hyman's procedure is followed for the entropy inequality condition. Flux limiter is used in the second-order accurate model that removes oscillations while maintaining the order of accuracy. The models are verified against available experimental data of a 2D flood wave due to partial dam-break. Numerical experiments are conducted to verify the models' ability to correctly predict behavior of the free surface, in addition to prediction of depth and velocity.