Index: /issm/trunk-jpl/src/c/bamg/Mesh.cpp
===================================================================
--- /issm/trunk-jpl/src/c/bamg/Mesh.cpp	(revision 21820)
+++ /issm/trunk-jpl/src/c/bamg/Mesh.cpp	(revision 21821)
@@ -3076,8 +3076,8 @@
 				IssmDouble area_3 =((bv.r.x -(*tcvj)(1)->r.x)*((*tcvj)(0)->r.y-(*tcvj)(1)->r.y)
 						- (bv.r.y -(*tcvj)(1)->r.y)*((*tcvj)(0)->r.x-(*tcvj)(1)->r.x));
-			//	if(area_1<0 || area_2<0 || area_3<0){
-			//		pointsoutside = true;
-			//		continue;
-			//	}
+				if(area_1<0 || area_2<0 || area_3<0){
+					pointsoutside = true;
+					continue;
+				}
 				if(!bv.GeomEdgeHook){
 					vertices[nbv].r              = bv.r;