Minor bikeshedding

[ranocha]

I found it helpful to return NamedTuples instead of separating values and names. For example, you could return a NamedTuple from integrals(analysis_callback), e.g.,

integrals(analysis_callback) = (waterheight_total = first_vector, velocity = second_vector, entropy=  third_vector)

Then, you can iterate over the pairs(integrals(analysis_callback)) for plotting

It's just an idea - your current setup works fine of course

[JoshuaLampert]

Sounds reasonable. Will do that later.