An accurate numerical solution to the Saint-Venant-Hirano model for mixed-sediment morphodynamics in rivers