Index: /issm/trunk-jpl/src/m/exp/flowlines.m
===================================================================
--- /issm/trunk-jpl/src/m/exp/flowlines.m	(revision 14439)
+++ /issm/trunk-jpl/src/m/exp/flowlines.m	(revision 14440)
@@ -38,6 +38,6 @@
 
 %check seed points
-%tria=TriaSearch(index,x,y,x0,y0);
-tria=tsearch(x,y,index,x0,y0);
+tria=TriaSearch(index,x,y,x0,y0);
+%tria=tsearch(x,y,index,x0,y0);
 pos=find(isnan(tria));
 x0(pos)=[];
@@ -70,6 +70,6 @@
 	%find current triangle
 	queue=find(~done);
-	%tria=TriaSearch(index,x,y,X(queue),Y(queue));
-	tria=tsearch(x,y,index,X(queue),Y(queue));
+	tria=TriaSearch(index,x,y,X(queue),Y(queue));
+	%tria=tsearch(x,y,index,X(queue),Y(queue));
 
 	%check that the point is actually inside a triangle of the mesh
@@ -112,4 +112,5 @@
 
 %same process but reverse (vel=-vel) to have a vcomplete flow line
+queue=[];
 counter=1;
 X=x0; Y=y0;
@@ -120,6 +121,6 @@
 	%find current triangle
 	queue=find(~done);
-	%tria=TriaSearch(index,x,y,X(queue),Y(queue));
-	tria=tsearch(x,y,index,X(queue),Y(queue));
+	tria=TriaSearch(index,x,y,X(queue),Y(queue));
+	%tria=tsearch(x,y,index,X(queue),Y(queue));
 
 	%check that the point is actually inside a triangle of the mesh