“…Bermudez and Vazquez [3] first introduced a concept of the "exact C-property". Since then, a number of well-balanced numerical methods have been developed for the SWEs, e.g., finite volume methods [1,3,18,38], finite difference/volume WENO methods [20,28,29,30], and discontinuous Galerkin (DG) methods [9,10,19,27,29,30,31,32]. DG methods have the advantages of highorder accuracy, high parallel efficiency, and flexibility for hp-adaptivity and arbitrary geometry and meshes.…”