“…Many shock capturing schemes with explicit time discretizations have been developed to solve the SWEs with source term (1.1), including high order finite difference [20,36,58,59,62], finite volume [1,10,34,38,43,44,64], residual distribution methods [46,47] and discontinuous Galerkin schemes [60,61,65,67], and many references therein. When solving the SWEs with source term numerically, it is important to preserve the exact conservation property (C-property) [35], namely, the nonzero flux gradient should be exactly balanced by the source term in the case of a stationary water.…”