Randall J. LeVeque

This is a big enough job that I think we should leave it for the next release, and put a note in that plotting with Python3 may have some problems...

Does this also work for making a colorbar when all patches are completely masked? This happens in what we're now testing with using the adjoint to flag for refinement and...

Also need to add documentation on using these tools.

@mjberger: Sorry, I was using mass computation as a simple example but the actual problem is to compute e.g. potential energy of a tsunami, which requires summing `eta**2 * dx*dy`...

This could be incorporated in the `gridtools.py` module I am working on to deal with AMR data.

This is timely since I've been working on making fgmax output a bit more flexible. In particular, I have a version that only writes one or two values per fgmax...

Regarding using `GeoClawData.read` on `manning_break`, the best way to handle reading this as a list might be to define a `read` function for this class that calls `clawutil.data.ClawData.read` and then...

@piyueh: Yes, I think you're right that there are some inconsistencies in the documentation, and probably unintended behavior in the code in some cases. Thanks for pointing this out. This...

Looking at the earlier lines in `qinit_module.f90`: https://github.com/clawpack/geoclaw/blob/3431f6596985ecf39041ec798e67a4f86611288a/src/2d/shallow/qinit_module.f90#L164-L170 I see that for `qinit_type < 4` (i.e. when specifying a perturbation to `h, hu,` or `hv`) the perturbation is only apply...

Oops -- thinking more about this, I think that we actually wanted `qinit_type == 4` to mean the *value* of eta is specified, not a *perturbation* as in the other...