Comment (by Roland Haas):

Summarizing my end of the conversation on mailing list (http://lists.einsteintoolkit.org/pipermail/users/2020-August/007565.html):

To make it fail I need more than 1 refinement level active. It fails on both 1 and 4 MPI ranks (that I have tested on).

Adding a bit of debug output:

std::cerr << "active1: " << active0 << " outputslab: " << outputslab << std::endl;

I get:

active1: bboxset<CCTK_INT4,3>(set<bbox>:{([12,12,12]:[136,136,136]:[4,4,4]/[3,3,3]:[34,
--

Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2451/reg-hdf5-output-of-quasilocalmeasures