fix(get_active) include the cbd layers when checking layer thickness …

…before BAS is loaded (#1051)
modflowpy · Feb 8, 2021 · 3ccdd57 · 3ccdd57
1 parent 5bfdf71
commit 3ccdd57
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/flopy/utils/check.py b/flopy/utils/check.py
@@ -436,7 +436,8 @@ def get_active(self, include_cbd=False):
         """
         mg = self.model.modelgrid
         if mg.grid_type == "structured":
-            inds = (mg.nlay, mg.nrow, mg.ncol)
+            nlaycbd = mg._StructuredGrid__laycbd.sum() if include_cbd else 0
+            inds = (mg.nlay + nlaycbd, mg.nrow, mg.ncol)
         elif mg.grid_type == "vertex":
             inds = (mg.nlay, mg.ncpl)
         else: