“…The success of Godunov-type flood models can be attributed to approximate Riemann solvers [37,75] which are also embedded in discontinuous Galerkin finite element schemes [4,26,45] and Boussinesq models that account for non-hydrostatic flow effects [47]. Godunov-type models have generally assumed either a structured mesh of quadrilateral cells [3,13,29,35,42,86] or an unstructured mesh of triangular cells [7,15,41,70,84]. The latter mandates greater overhead to track the neighborhood of data around each cell, and makes it more challenging to compute gradients in the solution because data points do not fall on a regular grid [7,41], but the unstructured mesh is very appealing for the ease with which meshes can be generated and tailored to the unique geometry of application sites and the ability to locally refine the mesh around areas of interest [11,21,54,57,66,67,76].…”