Abstract. Routing streamflow through a river network is a fundamental requirement to verify lateral water fluxes simulated by hydrologic and land surface models. River routing is performed at diverse resolutions ranging from few kilometres to 1∘. The presented multiscale routing model mRM calculates streamflow at diverse spatial and temporal resolutions. mRM solves the kinematic wave equation using a finite difference scheme. An adaptive time stepping scheme fulfilling a numerical stability criterion is introduced…