Index: /issm/trunk-jpl/src/m/plot/plot_BC.m
===================================================================
--- /issm/trunk-jpl/src/m/plot/plot_BC.m	(revision 14285)
+++ /issm/trunk-jpl/src/m/plot/plot_BC.m	(revision 14286)
@@ -7,29 +7,34 @@
 
 %plot dirichlets
-h1=plot3(...
-	md.mesh.x(find(~isnan(md.diagnostic.spcvx(1:md.mesh.numberofvertices,1)))),...
-	md.mesh.y(find(~isnan(md.diagnostic.spcvx(1:md.mesh.numberofvertices,1)))),...
-	md.mesh.z(find(~isnan(md.diagnostic.spcvx(1:md.mesh.numberofvertices,1)))),...
-	'ro','MarkerSize',14,'MarkerFaceColor','r');
-h2=plot3(...
-	md.mesh.x(find(~isnan(md.diagnostic.spcvy(1:md.mesh.numberofvertices,1)))),...
-	md.mesh.y(find(~isnan(md.diagnostic.spcvy(1:md.mesh.numberofvertices,1)))),...
-	md.mesh.z(find(~isnan(md.diagnostic.spcvy(1:md.mesh.numberofvertices,1)))),...
-	'bo','MarkerSize',10,'MarkerFaceColor','b');
-h3=plot3(...
-	md.mesh.x(find(~isnan(md.diagnostic.spcvz(1:md.mesh.numberofvertices,1)))),...
-	md.mesh.y(find(~isnan(md.diagnostic.spcvz(1:md.mesh.numberofvertices,1)))),...
-	md.mesh.z(find(~isnan(md.diagnostic.spcvz(1:md.mesh.numberofvertices,1)))),...
-	'yo','MarkerSize',6 ,'MarkerFaceColor','y');
+dirichleton=getfieldvalue(options,'dirichlet','on');
+if strcmpi(dirichleton,'on'),
+	h1=plot3(...
+		md.mesh.x(find(~isnan(md.diagnostic.spcvx(1:md.mesh.numberofvertices,1)))),...
+		md.mesh.y(find(~isnan(md.diagnostic.spcvx(1:md.mesh.numberofvertices,1)))),...
+		md.mesh.z(find(~isnan(md.diagnostic.spcvx(1:md.mesh.numberofvertices,1)))),...
+		'ro','MarkerSize',14,'MarkerFaceColor','r');
+	h2=plot3(...
+		md.mesh.x(find(~isnan(md.diagnostic.spcvy(1:md.mesh.numberofvertices,1)))),...
+		md.mesh.y(find(~isnan(md.diagnostic.spcvy(1:md.mesh.numberofvertices,1)))),...
+		md.mesh.z(find(~isnan(md.diagnostic.spcvy(1:md.mesh.numberofvertices,1)))),...
+		'bo','MarkerSize',10,'MarkerFaceColor','b');
+	h3=plot3(...
+		md.mesh.x(find(~isnan(md.diagnostic.spcvz(1:md.mesh.numberofvertices,1)))),...
+		md.mesh.y(find(~isnan(md.diagnostic.spcvz(1:md.mesh.numberofvertices,1)))),...
+		md.mesh.z(find(~isnan(md.diagnostic.spcvz(1:md.mesh.numberofvertices,1)))),...
+		'yo','MarkerSize',6 ,'MarkerFaceColor','y');
+end
 
 %update legend
 [legend_h,object_h,plot_h,text_strings]=legend();
 legend('off');
-text_strings{end+1}='vx Dirichlet';
-text_strings{end+1}='vy Dirichlet';
-if h3, text_strings{end+1}='vz Dirichlet'; end
-plot_h(end+1)=h1;
-plot_h(end+1)=h2;
-if h3, plot_h(end+1)=h3; end
+if strcmpi(dirichleton,'on'),
+	text_strings{end+1}='vx Dirichlet';
+	text_strings{end+1}='vy Dirichlet';
+	if h3, text_strings{end+1}='vz Dirichlet'; end
+	plot_h(end+1)=h1;
+	plot_h(end+1)=h2;
+	if h3, plot_h(end+1)=h3; end
+end
 legend(plot_h,text_strings,'location','NorthEast')