Index: /issm/trunk-jpl/src/m/exp/flowlines.py
===================================================================
--- /issm/trunk-jpl/src/m/exp/flowlines.py	(revision 24031)
+++ /issm/trunk-jpl/src/m/exp/flowlines.py	(revision 24032)
@@ -65,5 +65,4 @@
     flowpath['name'] = 'flowline{}'.format(0)
     for i in np.arange(1, len(flowpath['x'])):
-        print(i)
         flowpath['name'] = np.append(flowpath['name'], 'flowline{}'.format(i))
 
@@ -96,9 +95,10 @@
             while not all(done):
                 #find current triangle
-                queue = np.squeeze(np.where(done == 0))
+                queue = np.where(done == 0)[0]
                 tria = trifinder(X[queue], Y[queue])
                 #check that the point are actually inside a triangle of the mesh
-                outsiders = np.squeeze(np.where(tria < 0))
-                if len(outsiders) > 0:
+                outsiders = np.where(tria < 0)[0]
+                print('outsider is {} with queue {}'.format(outsiders, queue))
+                if np.size(outsiders) > 0:
                     for outsider in outsiders:
                         flowpath['x'][queue[outsider]] = np.delete(flowpath['x'][queue[outsider]][:], flowindex)