Index: /issm/trunk-jpl/src/m/plot/plot_BC.m
===================================================================
--- /issm/trunk-jpl/src/m/plot/plot_BC.m	(revision 21704)
+++ /issm/trunk-jpl/src/m/plot/plot_BC.m	(revision 21705)
@@ -7,4 +7,8 @@
 
 [x y z elements is2d isplanet]=processmesh(md,[],options);
+spcvx=processdata(md,md.stressbalance.spcvx,options);
+spcvy=processdata(md,md.stressbalance.spcvy,options);
+spcvz=processdata(md,md.stressbalance.spcvz,options);
+nbv = numel(x);
 
 %plot dirichlets
@@ -12,17 +16,17 @@
 if strcmpi(dirichleton,'on'),
 	h1=plot3(...
-		x(find(~isnan(md.stressbalance.spcvx(1:md.mesh.numberofvertices,1)))),...
-		y(find(~isnan(md.stressbalance.spcvx(1:md.mesh.numberofvertices,1)))),...
-		z(find(~isnan(md.stressbalance.spcvx(1:md.mesh.numberofvertices,1)))),...
+		x(find(~isnan(spcvx(1:nbv,1)))),...
+		y(find(~isnan(spcvx(1:nbv,1)))),...
+		z(find(~isnan(spcvx(1:nbv,1)))),...
 		'ro','MarkerSize',14,'MarkerFaceColor','r');
 	h2=plot3(...
-		x(find(~isnan(md.stressbalance.spcvy(1:md.mesh.numberofvertices,1)))),...
-		y(find(~isnan(md.stressbalance.spcvy(1:md.mesh.numberofvertices,1)))),...
-		z(find(~isnan(md.stressbalance.spcvy(1:md.mesh.numberofvertices,1)))),...
+		x(find(~isnan(spcvy(1:nbv,1)))),...
+		y(find(~isnan(spcvy(1:nbv,1)))),...
+		z(find(~isnan(spcvy(1:nbv,1)))),...
 		'bo','MarkerSize',10,'MarkerFaceColor','b');
 	h3=plot3(...
-		x(find(~isnan(md.stressbalance.spcvz(1:md.mesh.numberofvertices,1)))),...
-		y(find(~isnan(md.stressbalance.spcvz(1:md.mesh.numberofvertices,1)))),...
-		z(find(~isnan(md.stressbalance.spcvz(1:md.mesh.numberofvertices,1)))),...
+		x(find(~isnan(spcvz(1:nbv,1)))),...
+		y(find(~isnan(spcvz(1:nbv,1)))),...
+		z(find(~isnan(spcvz(1:nbv,1)))),...
 		'yo','MarkerSize',6 ,'MarkerFaceColor','y');
 end
Index: /issm/trunk-jpl/src/m/plot/plot_icefront.m
===================================================================
--- /issm/trunk-jpl/src/m/plot/plot_icefront.m	(revision 21704)
+++ /issm/trunk-jpl/src/m/plot/plot_icefront.m	(revision 21705)
@@ -19,5 +19,5 @@
 elementzeroice=sum(zeroice(md.mesh.elements),2);
 
-if dimension(md.mesh)==2,
+if is2d,
 	icefront=(elementice & elementnoice) & ~(elementice==2 & elementzeroice);
 
Index: /issm/trunk-jpl/src/m/plot/plot_manager.m
===================================================================
--- /issm/trunk-jpl/src/m/plot/plot_manager.m	(revision 21704)
+++ /issm/trunk-jpl/src/m/plot/plot_manager.m	(revision 21705)
@@ -21,5 +21,5 @@
 			plot_boundaries(md,options,subplotwidth,i);
 			return;
-		case 'BC',
+		case {'BC','bc'},
 			plot_BC(md,options,subplotwidth,i,data);
 			return;