“…We numerically solve the system (1.1)-(1.2), (1.4) and its one-dimensional (1-D) version by the Riemann-problem-solver-free central-upwind (CU) schemes, which were originally introduced in [25,27,28] for general multidimensional hyperbolic systems of conservative laws, and then extended to nonconservative hyperbolic systems in [6], where path-conservative CU (PCCU) schemes were introduced. The PCCU schemes were extended to the flux globalization framework allowing to treat a wider variety of nonconservative systems in [4,5,7,26].…”