Index: /issm/trunk-jpl/src/m/mesh/triangle.m
===================================================================
--- /issm/trunk-jpl/src/m/mesh/triangle.m	(revision 19948)
+++ /issm/trunk-jpl/src/m/mesh/triangle.m	(revision 19949)
@@ -41,19 +41,22 @@
 
 %check that all the created nodes belong to at least one element
-orphan=find(~ismember([1:length(x)],sort(unique(elements(:)))));
-for i=1:length(orphan),
-	disp('WARNING: removing orphans');
-	%get rid of the orphan node i
-	%update x and y
-	x=[x(1:orphan(i)-(i-1)-1); x(orphan(i)-(i-1)+1:end)];
-	y=[y(1:orphan(i)-(i-1)-1); y(orphan(i)-(i-1)+1:end)];
-	%update elements
-	pos=find(elements>orphan(i)-(i-1));
-	elements(pos)=elements(pos)-1;
-	%update segments
-	pos1=find(segments(:,1)>orphan(i)-(i-1));
-	pos2=find(segments(:,2)>orphan(i)-(i-1));
-	segments(pos1,1)=segments(pos1,1)-1;
-	segments(pos2,2)=segments(pos2,2)-1;
+removeorphans=1;
+if removeorphans,
+	orphan=find(~ismember([1:length(x)],sort(unique(elements(:)))));
+	for i=1:length(orphan),
+		disp('WARNING: removing orphans');
+		%get rid of the orphan node i
+		%update x and y
+		x=[x(1:orphan(i)-(i-1)-1); x(orphan(i)-(i-1)+1:end)];
+		y=[y(1:orphan(i)-(i-1)-1); y(orphan(i)-(i-1)+1:end)];
+		%update elements
+		pos=find(elements>orphan(i)-(i-1));
+		elements(pos)=elements(pos)-1;
+		%update segments
+		pos1=find(segments(:,1)>orphan(i)-(i-1));
+		pos2=find(segments(:,2)>orphan(i)-(i-1));
+		segments(pos1,1)=segments(pos1,1)-1;
+		segments(pos2,2)=segments(pos2,2)-1;
+	end
 end