Abstract
A novel finite volume method has been presented to solve the shallow water equations on triangular unstructured grid. Keeping both PV and VIA as the model variables and updating the boundary PVs through a local Riemann solver make the formulation much easier to be extended to unstructured grids with local reconstruction over single grid cell. We present a two dimensional formulation for shallow water equations on triangular unstructured grid of high order accuracy. Numerical oscillation can be effectively eliminated through the use of the reconstructionsv of oscillation free interpolation functions.