Fix(_mg_resync): Added checks to reset the modelgrid resync

autotest/t007_test.py

@@ -650,6 +650,11 @@ def test_sr():
     if mm.modelgrid.proj4 != "test test test":
         raise AssertionError()
 
+    mm.dis.top = 5000
+
+    if not np.allclose(mm.dis.top.array, mm.modelgrid.top):
+        raise AssertionError("modelgrid failed dynamic update test")
+
 
 def test_dis_sr():
 
@@ -683,6 +688,42 @@ def test_dis_sr():
     np.testing.assert_almost_equal(y, yul)
 
 
+def test_mf6_modelgrid_update():
+    base_dir = os.path.join("..", "examples", "data", "mf6")
+    # dis
+    model_ws = os.path.join(base_dir, "test001a_Tharmonic")
+    sim = flopy.mf6.MFSimulation.load(sim_ws=model_ws)
+    gwf = sim.get_model("flow15")
+
+    mg = gwf.modelgrid
+    gwf.dis.top = 12
+
+    if not np.allclose(gwf.dis.top.array, gwf.modelgrid.top):
+        raise AssertionError("StructuredGrid failed dynamic update test")
+
+    # disv
+    model_ws = os.path.join(base_dir, "test003_gwfs_disv")
+    sim = flopy.mf6.MFSimulation.load(sim_ws=model_ws)
+    gwf = sim.get_model("gwf_1")
+
+    mg = gwf.modelgrid
+    gwf.disv.top = 6.12
+
+    if not np.allclose(gwf.disv.top.array, gwf.modelgrid.top):
+        raise AssertionError("VertexGrid failed dynamic update test")
+
+    # disu
+    model_ws = os.path.join(base_dir, "test006_gwf3")
+    sim = flopy.mf6.MFSimulation.load(sim_ws=model_ws)
+    gwf = sim.get_model("gwf_1")
+
+    mg = gwf.modelgrid
+    gwf.disu.top = 101
+
+    if not np.allclose(gwf.disu.top.array, gwf.modelgrid.top):
+        raise AssertionError("UnstructuredGrid failed dynamic update test")
+
+
 def test_twri_mg():
     name = "twri.nam"
     ml = flopy.modflow.Modflow.load(name, model_ws=pth, check=False)
@@ -1666,6 +1707,7 @@ def main():
     # test_tricontour_NaN()
     # test_export_contourf()
     test_sr()
+    # test_mf6_modelgrid_update()
     # test_shapefile_polygon_closed()
     # test_mapview_plot_bc()
     # test_crosssection_plot_bc()

flopy/mf6/mfpackage.py

@@ -1614,6 +1614,14 @@ def __setattr__(self, name, value):
                         package=self._get_pname(),
                     )
                 return
+
+        if all(
+            hasattr(self, attr) for attr in ["model_or_sim", "_package_type"]
+        ):
+            if hasattr(self.model_or_sim, "_mg_resync"):
+                if not self.model_or_sim._mg_resync:
+                    self.model_or_sim._mg_resync = self._mg_resync
+
         super().__setattr__(name, value)
 
     def __repr__(self):

flopy/pakbase.py

@@ -78,6 +78,12 @@ def plottable(self):
     def has_stress_period_data(self):
         return self.__dict__.get("stress_period_data", None) is not None
 
+    @property
+    def _mg_resync(self):
+        if self.package_type.lower()[:4] in ("dis", "bas"):
+            return True
+        return False
+
     @staticmethod
     def _check_thresholds(chk, array, active, thresholds, name):
         """Checks array against min and max threshold values."""
@@ -605,6 +611,10 @@ def __setattr__(self, key, value):
                             )
                         value = new_list
 
+        if all(hasattr(self, attr) for attr in ["parent", "_name"]):
+            if not self.parent._mg_resync:
+                self.parent._mg_resync = self._mg_resync
+
         super().__setattr__(key, value)
 
     @property