Modifying the bounding box metric so that it uses the bounding box calculated in populate_state__compute_domain_statistics consistent with the periodic case.  This is then recalculated after every adapt rather than once per simulation.  Previously in adaptive parallel simulations the domain could resize (in my case quite significantly) after the load balancing.  This lead to strange constraints on my metric and the resolution suddenly changing over large areas of the domain.  As the load balancing occurs three times per adapt this fix isn't necessarily sufficient but hopefully, as the already bounded metric is interpolated and transferred, cutting this corner won't make much difference.