Describe the bug
Within utils/grindintesect.py function "def _intersect_linestring_shapely2" the function shapely.multilinestring() is called in line 969.
However this function doesn't exist and should be changed to shapely.multilinestrings()

To Reproduce
Shapely with a version above 2 has to be installed.
In specific cases if multiple parts of a linestring are within one cell this function will be called
This is quite rare, that is the reason I think this bug is not noticed before, since it has been present for quite some releases of flopy.

Expected behavior
I expected a multilinestring as output, however I got the error as given below

Screenshots