A higher order Finite Volume resolution method for a system related to the inviscid primitive equations in a complex domain

We construct the cell-centered Finite Volume discretization of the two-dimensional inviscid primitive equations in a domain with topography. To compute the numerical fluxes, the so-called Upwind Scheme (US) and the Central-Upwind Scheme (CUS) are introduced. For the time discretization, we use the classical fourth order Runge–Kutta method. We verify, with our numerical simulations, that the US (or CUS) is a robust first (or second) order scheme, regardless of the shape or size of the topography and without any mesh refinement near the topography.