For datasets in slab geometry or with no geometry applied, no regions are created. This might break some methods which in principle could still work if the whole grid was just treated as one region? Maybe apply_geometry should create a single, whole-grid region if no regions are created otherwise?
For datasets in slab geometry or with no geometry applied, no regions are created. This might break some methods which in principle could still work if the whole grid was just treated as one region? Maybe
apply_geometryshould create a single, whole-grid region if no regions are created otherwise?