We propose a new algorithm to solve the problem of flow routing across local depressions in the topography, one of the main computational bottlenecks in landscape evolution models. Our solution is more efficient than the state-of-the-art algorithms, with an optimal linear asymptotic complexity. The algorithm has been designed specifically to be used within landscape evolution models, and also suits more generally the efficient treatment of large digital elevation models.