Changeset 24032

Timestamp:

06/20/19 05:31:04 (6 years ago)

Author:

bdef

Message:

BUG: minor fix

File:

• issm/trunk-jpl/src/m/exp/flowlines.py (modified) (2 diffs)

issm/trunk-jpl/src/m/exp/flowlines.py

@@ -65 +65 @@
     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 +95 @@
             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)