I am really interested in using the `DifferentialEquations.jl`

framework to solve the shallow water equations over a 2D domain efficiently instead of coding my own solution.

How would you formulate the problem in `DifferentialEquations.jl`

ignoring all frictional, viscous and Coriolis terms, and only considering gravity? The corresponding system of equations for this set of assumptions in conservative form is in the Wikipedia link above.

In terms of the implementation: given a 2D map (or Julia array) of the landscape `H`

and a map indicating the depth of water `h_t`

on top of the landscape at time `t`

, I would like to advance in time and obtain `h_{t+1}`

iteratively until reach a time span [0,T].

I have a finite volume solution that is working for now, but I am not happy reinventing the wheel. It would be great to compare it with the solution obtained with `DifferentialEquations.jl`

and use that instead.

@ChrisRackauckas can you give a hand?