Hello Ian, Erik,
Could the routine GetAllActive be replaced by just accessing level_dboxes::active?
It could yes. GetAllActive computes level_dboxes::active (it dpulicates code from dh.cc, expanding the code in one of the bboxset constructors taking a vector and a pointer-to-member as arguments). I am not sure though if this information is actually stored anywhere outside of dh::regrid (in level_boxes and the various references to parts of it). Erik: is the "active" information stored permanently anywhere?
Yours, Roland