We present a new coupled numerical solution to the Saint-Venant-Hirano model.
We develop a concept of balancedness for the vertical sediment mass fluxes.
We propose a fully-mathematical benchmark based on linearised wave dynamics.
We model streamwise and vertical sorting by comparison with laboratory data.
Our accurate solution allows us to identify strengths and weaknesses in the model.