Index: /issm/trunk-jpl/src/m/geometry/FlagElements.m
===================================================================
--- /issm/trunk-jpl/src/m/geometry/FlagElements.m	(revision 15989)
+++ /issm/trunk-jpl/src/m/geometry/FlagElements.m	(revision 15990)
@@ -2,5 +2,5 @@
 %FLAGELEMENTS - flag the elements in an region
 %
-%   The region can be given with an exp file, a list of elements.
+%   The region can be given with an exp file, a list of elements or vertices
 %
 %   Usage: 
@@ -46,9 +46,12 @@
 		end
 	elseif isfloat(region) | islogical(region),
-		if size(region,1)~=md.mesh.numberofelements,
+		if size(region,1)==md.mesh.numberofelements,
+			flag=region;
+		elseif size(region,1)==md.mesh.numberofvertices,
+			flag=logical(sum(region(md.mesh.elements)>0,2)==size(md.mesh.elements,2));
+		else
 			help FlagElements
 			error('Flaglist for region must be of same size as number of elements in model');
 		end
-		flag=region;
 	else
 		error('Invalid region option');
Index: /issm/trunk-jpl/src/m/geometry/FlagElements.py
===================================================================
--- /issm/trunk-jpl/src/m/geometry/FlagElements.py	(revision 15989)
+++ /issm/trunk-jpl/src/m/geometry/FlagElements.py	(revision 15990)
@@ -10,5 +10,5 @@
 	FLAGELEMENTS - flag the elements in an region
 
-	   The region can be given with an exp file, a list of elements.
+	   The region can be given with an exp file, a list of elements or vertices
 
 	   Usage: 
@@ -54,7 +54,10 @@
 
 	elif isinstance(region,numpy.ndarray) or isinstance(region,bool):
-		if not numpy.size(region,0)==md.mesh.numberofelements:
+		if numpy.size(region,0)==md.mesh.numberofelements:
+			flag=region
+		else if numpy.size(region,0)==md.mesh.numberofvertices:
+			raise TypeError("look at FlagElements.m and implement the same...")
+		else:
 			raise TypeError("Flaglist for region must be of same size as number of elements in model.")
-		flag=region
 
 	else: