Kyle Mandli

First idea sounds good to me. Maybe a broader consistency check but we kind of do that already I suppose with `num_eqn`, etc. The `aux` arrays may not be able...

We did have something that would effectively call a Riemann solver in a mode that would be an actual Riemann problem. Something like that would be useful although I am...

@rgommers pointed out some issues in #252 and more from a discussion we had recently. I think `reconstruct.c` was specifically identified as being an issue.

It definitely looks like something wrong in the patch boundaries. I doubt it's related to #653 as that only happens if the aux arrays are time dependent. I have seen...

To look at the bathymetry on the entire patch, `b4step` is a good place to do that. I have a sneaking suspicion that the representation of the bathymetry is probably...

Which version of the code are you working with?

Possibly, there were some changes we made due to weirdness around how ghost cells and aux were filled. I had made a comment above about that only being relevant to...

There are a lot of parameters that can be changed. Can you provide more detail as to how/where things are being over estimated?

Thanks @sourabhmet. Can you provide the location of one of the points that you are seeing the 12 m surge? Also, can you provide a plot of the surface height...

Well that does make it difficult. From the Harvey example there's no where that I can see that produces more than 2 meters of surge (and it's more like 1...