Bathymetry support! The current implementation is a first order well balanced hydrostatic reconstruction scheme with Piecewise Constant reconstruction for the bottom topography. It preserves the lake at rest for Flat Bottom, No-Bottom (links to the old numerical solver without any sources) and Gaussian Hill bottom.