Add sill example

This adds an example and test of a pure Python version of a shallow water fwave Riemann solver (see PR https://github.com/clawpack/riemann/pull/82).