Fix(features_to_shapefile): fix missing bracket around linestring for…

… shapefile (#1410) * Fix(features_to_shapefile): fix missing bracket around linestring for shapefile * remove debugging code * update line representation for consistency with __geo_interface__ and flopy geometery LineString input geometry documentation * update(gridgen examples): update linestring type examples to be consistent with flopy geometry objects and standard spatial __geo_reference__ formats * Close #1405
modflowpy · May 6, 2022 · cadd216 · cadd216
1 parent 45a9e57
commit cadd216
Show file tree

Hide file tree

Showing 5 changed files with 78 additions and 98 deletions.
diff --git a/autotest/t061_test_gridgen.py b/autotest/t061_test_gridgen.py
@@ -196,7 +196,7 @@ def test_gridgen():
         assert n == 308, msg
 
         # test the gridgen line intersection
-        line = [[[(Lx, Ly), (Lx, 0.0)]]]
+        line = [[(Lx, Ly), (Lx, 0.0)]]
         cells = g.intersect(line, "line", 0)
         nlist = [n for n in cells["nodenumber"]]
         nlist2 = [

diff --git a/examples/Notebooks/flopy3_gridgen.ipynb b/examples/Notebooks/flopy3_gridgen.ipynb
diff --git a/examples/Notebooks/flopy3_modpath7_unstructured_example.ipynb b/examples/Notebooks/flopy3_modpath7_unstructured_example.ipynb
diff --git a/examples/common/setup_pmv_demo.py b/examples/common/setup_pmv_demo.py
@@ -190,7 +190,7 @@ def run():
         gwf, recharge=0.005, auxiliary=[("iface",)], aux={0: [6]}
     )
     # riv
-    riverline = [[[(Lx - 1.0, Ly), (Lx - 1.0, 0.0)]]]
+    riverline = [[(Lx - 1.0, Ly), (Lx - 1.0, 0.0)]]
     rivcells = g.intersect(riverline, "line", 0)
     rivspd = [
         [(0, icpl), 320.0, 100000.0, 318] for icpl in rivcells["nodenumber"]

diff --git a/flopy/utils/gridgen.py b/flopy/utils/gridgen.py
@@ -78,7 +78,7 @@ def features_to_shapefile(features, featuretype, filename):
         wr = shapefile.Writer(filename, shapeType=shapefile.POLYLINE)
         wr.field("SHAPEID", "N", 20, 0)
         for i, line in enumerate(features):
-            wr.line(line.__geo_interface__["coordinates"])
+            wr.line([line.__geo_interface__["coordinates"]])
             wr.record(i)
 
     elif featuretype.lower() == "point":