Use ellipses q[0,...] for depth, which is q[0,:] in 1d or q[0,:,:] in 2d.
Use eta = q[-1,...] to work in 1d or 2d, and also with 2d bouss where q will have additional components for hu and hv correction terms. Assumes that eta is always the final component of the q array read in from geoclaw output.
@mandli: Not sure what sort of example you mean since we don't have any examples in visclaw repo now. For old examples in geoclaw it should work the same as before, I hope.
Use ellipses q[0,...] for depth, which is q[0,:] in 1d or q[0,:,:] in 2d.
Use eta = q[-1,...] to work in 1d or 2d, and also with 2d bouss where q will have additional components for hu and hv correction terms. Assumes that eta is always the final component of the q array read in from geoclaw output.