Merge pull request #673 from rjleveque/make_shoreline_noplot

suppress plot made by make_shoreline_xy

Author: rjleveque

src/python/geoclaw/topotools.py

@@ -1716,13 +1716,15 @@ def make_shoreline_xy(self, sea_level=0):
         x = self.x
         y = self.y
         Z = self.Z
+        fig = plt.figure()
         c = plt.contour(x,y,Z,[sea_level])
         # c is the level 0 contour as list of arrays, one for each segement
         # catenate these together separated by array([nan,nan]):
         shoreline_xy = c.allsegs[0][0]  # first segment
         for k in range(1,len(c.allsegs[0])):
             shoreline_xy = numpy.vstack((shoreline_xy, \
                            numpy.array([numpy.nan,numpy.nan]), c.allsegs[0][k]))
+        plt.close(fig)
         return shoreline_xy