Index: /issm/trunk-jpl/src/m/plot/applyoptions.m
===================================================================
--- /issm/trunk-jpl/src/m/plot/applyoptions.m	(revision 24027)
+++ /issm/trunk-jpl/src/m/plot/applyoptions.m	(revision 24028)
@@ -51,5 +51,5 @@
 set(gca,'FontSize',getfieldvalue(options,'axisfontsize',fontsize));;
 if exist(options,'axis')
-	eval(['axis ' getfieldvalue(options,'axis')]);
+	axis(getfieldvalue(options,'axis'));
 else
 	if strcmp(domaintype(md.mesh),'3D'),
@@ -73,5 +73,5 @@
 %box
 if exist(options,'box')
-	eval(['box ' getfieldvalue(options,'box')]);
+	box(getfieldvalue(options,'box'));
 end
 
@@ -243,5 +243,4 @@
 else
 	%do nothing
-
 end
 
Index: /issm/trunk-jpl/src/m/plot/plot_elementnumbering.m
===================================================================
--- /issm/trunk-jpl/src/m/plot/plot_elementnumbering.m	(revision 24027)
+++ /issm/trunk-jpl/src/m/plot/plot_elementnumbering.m	(revision 24028)
@@ -25,7 +25,14 @@
 
 	%numbering
-	text(sum(x(elements(:,1:3)),2)/3,sum(y(elements(:,1:3)),2)/3,sum(z(elements(:,1:3)),2)/3,...
-		num2str(elementnumbers),...
-		'HorizontalAlignment','center','VerticalAlignment','middle');
+	X   = sum(x(elements(:,1:3)),2)/3;
+	Y   = sum(y(elements(:,1:3)),2)/3;
+	Z   = sum(z(elements(:,1:3)),2)/3;
+	NUM = elementnumbers;
+	if exist(options,'axis')
+		AXIS = getfieldvalue(options,'axis');
+		pos = find(X>AXIS(1) & X<AXIS(2) & Y>AXIS(3) & Y<AXIS(4));
+		X = X(pos); Y=Y(pos); Z=Z(pos); NUM=NUM(pos);
+	end
+	text(X,Y,Z,num2str(NUM),'HorizontalAlignment','center','VerticalAlignment','middle');
 else
 	%plot mesh 
@@ -47,7 +54,14 @@
 
 	%numbering
-	text(sum(x(elements(:,1:6)),2)/6,sum(y(elements(:,1:6)),2)/6,sum(z(elements(:,1:6)),2)/6,...
-		num2str(elementnumbers),...
-		'HorizontalAlignment','center','VerticalAlignment','middle');
+	X   = sum(x(elements(:,1:6)),2)/6;
+	Y   = sum(y(elements(:,1:6)),2)/6;
+	Z   = sum(z(elements(:,1:6)),2)/6;
+	NUM = elementnumbers;
+	if exist(options,'axis')
+		AXIS = getfieldvalue(options,'axis');
+		pos = find(X>AXIS(1) & X<AXIS(2) & Y>AXIS(3) & Y<AXIS(4));
+		X = X(pos); Y=Y(pos); Z=Z(pos); NUM=NUM(pos);
+	end
+	text(X,Y,Z,num2str(NUM),'HorizontalAlignment','center','VerticalAlignment','middle');
 end