Index: /issm/trunk/src/m/classes/public/modelextract.m
===================================================================
--- /issm/trunk/src/m/classes/public/modelextract.m	(revision 2986)
+++ /issm/trunk/src/m/classes/public/modelextract.m	(revision 2987)
@@ -34,35 +34,5 @@
 
 %first
-if ischar(area),
-	if isempty(area),
-		flag_elem=zeros(md1.numberofelements,1);
-		invert=0;
-	elseif strcmpi(area,'all')
-		flag_elem=ones(md1.numberofelements,1);
-		invert=0;
-	else
-		%make sure that we actually don't want the elements outside the domain outline!
-		if strcmpi(area(1),'~'),
-			area=area(2:length(area));
-			invert=1;
-		else
-			invert=0;
-		end
-
-		%ok, flag_elem elements
-		A=expread(area,1);
-		flag_elem=ContourToMesh(md1.elements(:,1:3),md1.x,md1.y,A,'element',1);
-	end
-	if invert, flag_elem=~flag_elem; end
-
-elseif isfloat(area),
-	if size(area,1)~=md1.numberofelements,
-		setelementstypeusage();
-		error('Flags must be of same size as number of elements in model');
-	end
-	flag_elem=area;
-else
-	error('Invalide option');
-end
+flag_elem=FlagElements(md,area);
 
 %kick out all elements with 3 dirichlets