Index: /issm/trunk/src/c/modules/TriaSearchx/TriaSearchx.cpp
===================================================================
--- /issm/trunk/src/c/modules/TriaSearchx/TriaSearchx.cpp	(revision 5390)
+++ /issm/trunk/src/c/modules/TriaSearchx/TriaSearchx.cpp	(revision 5391)
@@ -33,5 +33,4 @@
 
 	// read background mesh 
-	if (verbose) printf("Reading mesh\n");
 	Mesh Th(index,x,y,nods,nel); 
 	Th.CreateSingleVertexToTriangleConnectivity();
@@ -48,20 +47,10 @@
 
 		// internal point 
-		if (tb.det>0){ 
-			//Area coordinate
-			areacoord[0]= (double) dete[0]/ tb.det;
-			areacoord[1]= (double) dete[1] / tb.det;
-			areacoord[2]= (double) dete[2] / tb.det;
-			//3 vertices of the triangle
-			i0=Th.GetId(tb[0]);
-			i1=Th.GetId(tb[1]);
-			i2=Th.GetId(tb[2]);
-			//triangle number
-			tria[i]=(double)Th.GetId(tb);
-		}
+		if (Th.GetId(tb)<nel)tria[i]=(double)Th.GetId(tb);
 		//external point
 		else tria[i]=NAN;
 	}
 
+
 	/*Assign output pointers:*/
 	*ptria=tria;