Catch negative heights produced by solver in WavePropagation
Move boundary type (outflow or reflecting) decision to scenarios
Add ReflectedShock scenario