add inactive inhoms bounds check

some tiny notebook improvements
mbakker7 · Dec 11, 2024 · 00bfa45 · 00bfa45
1 parent 6280815
commit 00bfa45
Show file tree

Hide file tree

Showing 3 changed files with 15 additions and 213 deletions.
diff --git a/docs/03examples/1d_elements.ipynb b/docs/03examples/1d_elements.ipynb
diff --git a/docs/03examples/1d_inhom_elements.ipynb b/docs/03examples/1d_inhom_elements.ipynb
diff --git a/ttim/model.py b/ttim/model.py
@@ -922,6 +922,14 @@ def check_inhoms(self):
             raise ValueError(
                 f"Number of aquifers does not match {self.aq.naq}:\n{naqs}"
             )
+        # # shared boundary check
+        # # NOTE: does not deal with nested inhoms
+        # xcoords = np.concatenate(
+        #     [(inhom.x1, inhom.x2) for inhom in self.aq.inhomdict.values()]
+        # )
+        # xcoords.sort()
+        # if not np.all(np.diff(xcoords[1:-1])[::2] < 1e-10):
+        #     raise ValueError("Not all inhomogeneities have shared boundaries.")
 
     def initialize(self):
         self.check_inhoms()