should be ncpl not nnodes

modflowpy · Apr 15, 2024 · 95248e6 · 95248e6
1 parent 5f60c09
commit 95248e6
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/flopy/utils/gridintersect.py b/flopy/utils/gridintersect.py
@@ -443,14 +443,14 @@ def _vtx_grid_to_geoms_cellids(self):
                         )
                     )
                 )
-                for node in range(self.mfgrid.nnodes)
+                for node in range(self.mfgrid.ncpl)
             ]
         else:
             geoms = [
                 shapely.polygons(self.mfgrid.get_cell_vertices(node))
-                for node in range(self.mfgrid.nnodes)
+                for node in range(self.mfgrid.ncpl)
             ]
-        return np.array(geoms), np.arange(self.mfgrid.nnodes)
+        return np.array(geoms), np.arange(self.mfgrid.ncpl)
 
     def _rect_grid_to_shape_list(self):
         """internal method, list of shapely polygons for structured grid cells.