There's a new [advection component](https://github.com/landlab/landlab/pull/1582) in Landlab that can be used in the solution to the volcano problem.