Index: /issm/trunk/test/Miscellaneous/Bump/Bump.par
===================================================================
--- /issm/trunk/test/Miscellaneous/Bump/Bump.par	(revision 8823)
+++ /issm/trunk/test/Miscellaneous/Bump/Bump.par	(revision 8824)
@@ -28,6 +28,6 @@
 %Deal with boundary conditions:
 md=SetIceSheetBC(md);
-md.spcvelocity(:,4)=zeros(md.numberofnodes,1);
-md.spcvelocity(:,5)=200*ones(md.numberofnodes,1);
+md.spcvx=zeros(md.numberofnodes,1);
+md.spcvy=200*ones(md.numberofnodes,1);
 
 %Parallel options
Index: /issm/trunk/test/Miscellaneous/Regularization/Square.par
===================================================================
--- /issm/trunk/test/Miscellaneous/Regularization/Square.par	(revision 8823)
+++ /issm/trunk/test/Miscellaneous/Regularization/Square.par	(revision 8824)
@@ -25,10 +25,7 @@
 md=SetIceSheetBC(md);
 pos=find(md.gridonboundary);
-md.spcvelocity(pos,1)=1;
-md.spcvelocity(pos,2)=1;
-md.spcvelocity(pos,3)=1;
-md.spcvelocity(pos,4)=0;
-md.spcvelocity(pos,5)=200;
-md.spcvelocity(pos,6)=0;
+md.spcvx(pos)=0;
+md.spcvy(pos)=200;
+md.spcvz(pos)=0;
 
 md.eps_res=10^-5;
Index: /issm/trunk/test/NightlyRun/test1101.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1101.m	(revision 8823)
+++ /issm/trunk/test/NightlyRun/test1101.m	(revision 8824)
@@ -21,7 +21,11 @@
 
 	%Create dirichlet on the bed only
-	md.spcvelocity=zeros(md.numberofnodes,6);
+	md.spcvx=NaN*ones(md.numberofnodes,1);
+	md.spcvy=NaN*ones(md.numberofnodes,1);
+	md.spcvz=NaN*ones(md.numberofnodes,1);
+
 	pos=find(md.nodeonbed);
-	md.spcvelocity(pos,1:2)=1;
+	md.spcvx(pos)=0;
+	md.spcvy(pos)=0;
 
 	%Create MPCs to have periodic boundary conditions
Index: /issm/trunk/test/NightlyRun/test1102.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1102.m	(revision 8823)
+++ /issm/trunk/test/NightlyRun/test1102.m	(revision 8824)
@@ -28,7 +28,8 @@
 
 	%Create dirichlet on the bed only
-	%md.spcvelocity=zeros(md.numberofnodes,6);
 	pos=find(md.nodeonbed);
-	md.spcvelocity(pos,1:3)=1;
+	md.spcvx(pos)=0;
+	md.spcvy(pos)=0;
+	md.spcvz(pos)=0;
 
 	%Create MPCs to have periodic boundary conditions
Index: /issm/trunk/test/NightlyRun/test1103.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1103.m	(revision 8823)
+++ /issm/trunk/test/NightlyRun/test1103.m	(revision 8824)
@@ -19,7 +19,10 @@
 
 	%Create dirichlet on the bed only
-	md.spcvelocity=zeros(md.numberofnodes,6);
+	md.spcvx=NaN*ones(md.numberofnodes,1);
+	md.spcvy=NaN*ones(md.numberofnodes,1);
+	md.spcvz=NaN*ones(md.numberofnodes,1);
 	pos=find(md.nodeonbed);
-	md.spcvelocity(pos,1:2)=1;
+	md.spcvx(pos)=0;
+	md.spcvy(pos)=0;
 
 	%Create MPCs to have periodic boundary conditions
Index: /issm/trunk/test/NightlyRun/test1104.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1104.m	(revision 8823)
+++ /issm/trunk/test/NightlyRun/test1104.m	(revision 8824)
@@ -18,7 +18,11 @@
 
 	%Create dirichlet on the bed only
-	md.spcvelocity=zeros(md.numberofnodes,6);
+	md.spcvx=NaN*ones(md.numberofnodes,1);
+	md.spcvy=NaN*ones(md.numberofnodes,1);
+	md.spcvz=NaN*ones(md.numberofnodes,1);
+
 	pos=find(md.nodeonbed);
-	md.spcvelocity(pos,1:2)=1;
+	md.spcvx(pos)=0;
+	md.spcvy(pos)=0;
 
 	%Create MPCs to have periodic boundary conditions
Index: /issm/trunk/test/NightlyRun/test1105.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1105.m	(revision 8823)
+++ /issm/trunk/test/NightlyRun/test1105.m	(revision 8824)
@@ -21,5 +21,7 @@
 
 	%Create MPCs to have periodic boundary conditions
-	md.spcvelocity=zeros(md.numberofnodes,6);
+	md.spcvx=NaN*ones(md.numberofnodes,1);
+	md.spcvy=NaN*ones(md.numberofnodes,1);
+	md.spcvz=NaN*ones(md.numberofnodes,1);
 
 	posx=find(md.x==0 & md.y~=0 & md.y~=L);
@@ -33,28 +35,29 @@
 	%Add spc on the corners
 	pos=find((md.x==0 | md.x==L) & (md.y==0 | md.y==L) & md.nodeonbed);
-	md.spcvelocity(pos,1:2)=1;
+	md.spcvx(pos)=0;
+	md.spcvy(pos)=0;
 	if(L==5000),
-		md.spcvelocity(pos,4)=15.66;
-		md.spcvelocity(pos,5)=-0.1967;
+		md.spcvx(pos)=15.66;
+		md.spcvy(pos)=-0.1967;
 	elseif(L==10000),
-		md.spcvelocity(pos,4)=16.04;
-		md.spcvelocity(pos,5)=-0.1977;
+		md.spcvx(pos)=16.04;
+		md.spcvy(pos)=-0.1977;
 	elseif(L==20000),
-		md.spcvelocity(pos,4)=16.53;
-		md.spcvelocity(pos,5)=-1.27;
+		md.spcvx(pos)=16.53;
+		md.spcvy(pos)=-1.27;
 	elseif(L==40000),
-		md.spcvelocity(pos,4)=17.23;
-		md.spcvelocity(pos,5)=-3.17;
+		md.spcvx(pos)=17.23;
+		md.spcvy(pos)=-3.17;
 	elseif(L==80000),
-		md.spcvelocity(pos,4)=16.68;
-		md.spcvelocity(pos,5)=-2.69;
+		md.spcvx(pos)=16.68;
+		md.spcvy(pos)=-2.69;
 	elseif(L==160000),
-		md.spcvelocity(pos,4)=16.03;
-		md.spcvelocity(pos,5)=-1.27;
+		md.spcvx(pos)=16.03;
+		md.spcvy(pos)=-1.27;
 	end
 	
 	%Spc the bed at zero for vz
 	pos=find(md.nodeonbed);
-	md.spcvelocity(pos,3)=1;
+	md.spcvz(pos)=0;
 
 	%Compute the diagnostic
Index: /issm/trunk/test/NightlyRun/test1106.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1106.m	(revision 8823)
+++ /issm/trunk/test/NightlyRun/test1106.m	(revision 8824)
@@ -18,5 +18,4 @@
 
 	%Create MPCs to have periodic boundary conditions
-	%md.spcvelocity=zeros(md.numberofnodes,6);
 
 	%posx=find(md.x==0);
Index: /issm/trunk/test/NightlyRun/test1107.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1107.m	(revision 8823)
+++ /issm/trunk/test/NightlyRun/test1107.m	(revision 8824)
@@ -21,5 +21,7 @@
 
 	%We need one grd on dirichlet: the 4 corners are set to zero
-	md.spcvelocity=zeros(md.numberofnodes,6);
+	md.spcvx=NaN*ones(md.numberofnodes,1);
+	md.spcvy=NaN*ones(md.numberofnodes,1);
+	md.spcvz=NaN*ones(md.numberofnodes,1);
 
 	%Create MPCs to have periodic boundary conditions
@@ -34,23 +36,23 @@
 	%Add spc on the corners
 	pos=find((md.x==0 | md.x==L) & (md.y==0 | md.y==L) & md.nodeonbed);
-	md.spcvelocity(:,2)=1;
-	md.spcvelocity(pos,1:2)=1;
+	md.spcvy(:)=0;
+	md.spcvx(pos)=0;
 	if(L==5000),
-		md.spcvelocity(pos,4)=16.0912;
+		md.spcvx(pos)=16.0912;
 	elseif(L==10000),
-		md.spcvelocity(pos,4)=16.52;
+		md.spcvx(pos)=16.52;
 	elseif(L==20000),
-		md.spcvelocity(pos,4)=17.77;
+		md.spcvx(pos)=17.77;
 	elseif(L==40000),
-		md.spcvelocity(pos,4)=19.88;
+		md.spcvx(pos)=19.88;
 	elseif(L==80000),
-		md.spcvelocity(pos,4)=18.65;
+		md.spcvx(pos)=18.65;
 	elseif(L==160000),
-		md.spcvelocity(pos,4)=16.91;
+		md.spcvx(pos)=16.91;
 	end
 	
 	%Spc the bed at zero for vz
 	pos=find(md.nodeonbed);
-	md.spcvelocity(pos,3)=1;
+	md.spcvz(pos,3)=0;
 
 	%Compute the diagnostic
Index: /issm/trunk/test/NightlyRun/test1108.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1108.m	(revision 8823)
+++ /issm/trunk/test/NightlyRun/test1108.m	(revision 8824)
@@ -18,7 +18,12 @@
 
 	%We need one grd on dirichlet: the 4 corners are set to zero
-	%md.spcvelocity=zeros(md.numberofnodes,6);
+	%md.spcvx=NaN*ones(md.numberofnodes,1);
+	%md.spcvy=NaN*ones(md.numberofnodes,1);
+	%md.spcvz=NaN*ones(md.numberofnodes,1);
+	
 	%pos=find((md.x==0 | md.x==max(md.x)) & (md.y==0 | md.y==max(md.y)));
-	%md.spcvelocity(pos,1:3)=1;
+	%md.spcvx(pos)=0;
+	%md.spcvy(pos)=0;
+	%md.spcvz(pos)=0;
 
 	%%Create MPCs to have periodic boundary conditions
Index: /issm/trunk/test/NightlyRun/test1109.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1109.m	(revision 8823)
+++ /issm/trunk/test/NightlyRun/test1109.m	(revision 8824)
@@ -32,11 +32,17 @@
 	%Create spcs on the bed 
 	pos=find(md.nodeonbed);
-	md.spcvelocity=zeros(md.numberofnodes,6);
-	md.spcvelocity(pos,1:3)=1;
+	md.spcvx=NaN*ones(md.numberofnodes,1);
+	md.spcvy=NaN*ones(md.numberofnodes,1);
+	md.spcvz=NaN*ones(md.numberofnodes,1);
+	md.spcvx(pos)=0;
+	md.spcvy(pos)=0;
+	md.spcvz(pos)=0;
 
 	%Remove the spc where there is some sliding (case 3 and 4):
 	if i==3 | i==4,
 		pos=find(md.y/max(md.y)>=0.44 & md.y/max(md.y)<=0.5);
-		md.spcvelocity(pos,1:3)=0;
+		md.spcvx(pos)=0;
+		md.spcvy(pos)=0;
+		md.spcvz(pos)=0;
 	end
 
Index: /issm/trunk/test/NightlyRun/test1110.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1110.m	(revision 8823)
+++ /issm/trunk/test/NightlyRun/test1110.m	(revision 8824)
@@ -20,14 +20,18 @@
 	end
 
-	md.spcvelocity(:,:)=0;
+	md.spcvx=NaN*ones(md.numberofnodes,1);
+	md.spcvy=NaN*ones(md.numberofnodes,1);
+	md.spcvz=NaN*ones(md.numberofnodes,1);
 	if (i==1 | i==3),
 		%Create dirichlet on the bed if no slip
 		pos=find(md.nodeonbed);
-		md.spcvelocity(pos,1:3)=1;
+		md.spcvx(pos)=0;
+		md.spcvy(pos)=0;
+		md.spcvz(pos)=0;
 	else
 		pos=find(md.nodeonbed & (md.x==0 | md.x==max(md.x)) & (md.y==0 | md.y==max(md.y)));
-		md.spcvelocity(pos,1:3)=1;
-		md.spcvelocity(pos,4)=100; %because we need a dirichlet somewhere
-
+		md.spcvx(pos)=100; %because we need a dirichlet somewhere
+		md.spcvy(pos)=0;
+		md.spcvz(pos)=0;
 	end
 	pos=find(~md.nodeonbed);
Index: /issm/trunk/test/NightlyRun/test120.m
===================================================================
--- /issm/trunk/test/NightlyRun/test120.m	(revision 8823)
+++ /issm/trunk/test/NightlyRun/test120.m	(revision 8824)
@@ -4,4 +4,5 @@
 md=setelementstype(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
+
 md=solve(md,TransientSolutionEnum);
 
Index: /issm/trunk/test/NightlyRun/test1201.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1201.m	(revision 8823)
+++ /issm/trunk/test/NightlyRun/test1201.m	(revision 8824)
@@ -18,8 +18,10 @@
 
 	%analytical test
-	md.spcvelocity(:,1:3)=1;
-	md.spcvelocity(:,4)=0;
-	md.spcvelocity(:,5)=-400;
-	md.spcvelocity(1,1)=0;
+	md.spcvx=zeros(md.numberofnodes,1);
+	md.spcvy=zeros(md.numberofnodes,1);
+	md.spcvz=zeros(md.numberofnodes,1);
+	md.spcvx(:)=0;
+	md.spcvy(:)=-400;
+	md.spcvx(1)=0;
 	if type==1,
 		md.artificial_diffusivity=1;
Index: /issm/trunk/test/NightlyRun/test1203.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1203.m	(revision 8823)
+++ /issm/trunk/test/NightlyRun/test1203.m	(revision 8824)
@@ -11,5 +11,5 @@
 %Impose a non zero velocity on the upper boundary condition (y=max(y))
 pos=find(md.y==max(md.y));
-md.spcvelocity(pos,5)=400*(((md.x(pos)-100000)/25000).^2-ones(size(pos,1),1)).*heaviside((1+eps)*ones(size(pos,1),1)-((md.x(pos)-100000)/25000).^2);
+md.spcvy(pos)=400*(((md.x(pos)-100000)/25000).^2-ones(size(pos,1),1)).*heaviside((1+eps)*ones(size(pos,1),1)-((md.x(pos)-100000)/25000).^2);
 
 %Compute solution for MacAyeal's model 
Index: /issm/trunk/test/NightlyRun/test1204.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1204.m	(revision 8823)
+++ /issm/trunk/test/NightlyRun/test1204.m	(revision 8824)
@@ -11,5 +11,5 @@
 %Impose a non zero velocity on the upper boundary condition (y=max(y))
 pos=find(md.y==max(md.y));
-md.spcvelocity(pos,5)=400*(((md.x(pos)-100000)/25000).^2-ones(size(pos,1),1)).*heaviside((1+eps)*ones(size(pos,1),1)-((md.x(pos)-100000)/25000).^2);
+md.spcvy(pos)=400*(((md.x(pos)-100000)/25000).^2-ones(size(pos,1),1)).*heaviside((1+eps)*ones(size(pos,1),1)-((md.x(pos)-100000)/25000).^2);
 
 %Compute solution for MacAyeal's model 
Index: /issm/trunk/test/NightlyRun/test1205.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1205.m	(revision 8823)
+++ /issm/trunk/test/NightlyRun/test1205.m	(revision 8824)
@@ -23,5 +23,7 @@
 %Spc the nodes on the bed
 pos=find(md.nodeonbed);
-md.spcvelocity(pos,1:3)=1;
+md.spcvx(pos)=0;
+md.spcvy(pos)=0;
+md.spcvz(pos)=0;
 
 %Now we can solve the problem 
Index: /issm/trunk/test/NightlyRun/test1206.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1206.m	(revision 8823)
+++ /issm/trunk/test/NightlyRun/test1206.m	(revision 8824)
@@ -23,5 +23,7 @@
 %Spc the nodes on the bed
 pos=find(md.nodeonbed);
-md.spcvelocity(pos,1:3)=1;
+md.spcvx(pos)=0;
+md.spcvy(pos)=0;
+md.spcvz(pos)=0;
 
 %Now we can solve the problem 
Index: /issm/trunk/test/NightlyRun/test1207.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1207.m	(revision 8823)
+++ /issm/trunk/test/NightlyRun/test1207.m	(revision 8824)
@@ -23,5 +23,7 @@
 %Spc the nodes on the bed
 pos=find(md.nodeonbed);
-md.spcvelocity(pos,1:3)=1;
+md.spcvx(pos)=0;
+md.spcvy(pos)=0;
+md.spcvz(pos)=0;
 
 %Now we can solve the problem 
Index: /issm/trunk/test/NightlyRun/test1208.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1208.m	(revision 8823)
+++ /issm/trunk/test/NightlyRun/test1208.m	(revision 8824)
@@ -14,5 +14,7 @@
 %Spc the nodes on the bed
 pos=find(md.nodeonbed);
-md.spcvelocity(pos,1:3)=1;
+md.spcvx(pos)=0;
+md.spcvy(pos)=0;
+md.spcvz(pos)=0;
 
 %Adapt the time steps to the resolution
Index: /issm/trunk/test/NightlyRun/test233.m
===================================================================
--- /issm/trunk/test/NightlyRun/test233.m	(revision 8823)
+++ /issm/trunk/test/NightlyRun/test233.m	(revision 8824)
@@ -38,15 +38,17 @@
 
 %Boundary conditions:
-md.spcvelocity=zeros(md.numberofnodes,6);
+md.spcvx=NaN*ones(md.numberofnodes,1);
+md.spcvy=NaN*ones(md.numberofnodes,1);
+md.spcvz=NaN*ones(md.numberofnodes,1);
 
 %constrain flanks to 0 normal velocity
 pos=find(md.x==xmin | md.x==xmax);
-md.spcvelocity(pos,1)=1;
-md.spcvelocity(pos,3)=0;
+md.spcvx(pos)=0;
+md.spcvz(pos)=NaN;
 
 %constrain grounding line to 0 velocity
 pos=find(md.y==ymin);
-md.spcvelocity(pos,1:2)=1;
-md.spcvelocity(pos,3:4)=0;
+md.spcvx(pos)=0;
+md.spcvy(pos)=0;
 
 %icefront
Index: /issm/trunk/test/NightlyRun/test234.m
===================================================================
--- /issm/trunk/test/NightlyRun/test234.m	(revision 8823)
+++ /issm/trunk/test/NightlyRun/test234.m	(revision 8824)
@@ -38,15 +38,17 @@
 
 %Boundary conditions:
-md.spcvelocity=zeros(md.numberofnodes,6);
+md.spcvx=NaN*ones(md.numberofnodes,1);
+md.spcvy=NaN*ones(md.numberofnodes,1);
+md.spcvz=NaN*ones(md.numberofnodes,1);
 
 %constrain flanks to 0 normal velocity
 pos=find(md.x==xmin | md.x==xmax);
-md.spcvelocity(pos,1)=1;
-md.spcvelocity(pos,3)=0;
+md.spcvx(pos)=0;
+md.spcvz(pos)=NaN;
 
 %constrain grounding line to 0 velocity
 pos=find(md.y==ymin);
-md.spcvelocity(pos,1:2)=1;
-md.spcvelocity(pos,3:4)=0;
+md.spcvx(pos)=0;
+md.spcvy(pos)=0;
 
 %icefront
Index: /issm/trunk/test/NightlyRun/test235.m
===================================================================
--- /issm/trunk/test/NightlyRun/test235.m	(revision 8823)
+++ /issm/trunk/test/NightlyRun/test235.m	(revision 8824)
@@ -35,15 +35,17 @@
 
 %Boundary conditions:
-md.spcvelocity=zeros(md.numberofnodes,6);
+md.spcvx=NaN*ones(md.numberofnodes,1);
+md.spcvy=NaN*ones(md.numberofnodes,1);
+md.spcvz=NaN*ones(md.numberofnodes,1);
 
 %constrain flanks to 0 normal velocity
 pos=find(md.x==xmin | md.x==xmax);
-md.spcvelocity(pos,1)=1;
-md.spcvelocity(pos,3)=0;
+md.spcvx(pos)=0;
+md.spcvz(pos)=NaN;
 
 %constrain grounding line to 0 velocity
 pos=find(md.y==ymin);
-md.spcvelocity(pos,1:2)=1;
-md.spcvelocity(pos,3:4)=0;
+md.spcvx(pos)=0;
+md.spcvy(pos)=0;
 
 %icefront
Index: /issm/trunk/test/NightlyRun/test236.m
===================================================================
--- /issm/trunk/test/NightlyRun/test236.m	(revision 8823)
+++ /issm/trunk/test/NightlyRun/test236.m	(revision 8824)
@@ -35,15 +35,17 @@
 
 %Boundary conditions:
-md.spcvelocity=zeros(md.numberofnodes,6);
+md.spcvx=NaN*ones(md.numberofnodes,1);
+md.spcvy=NaN*ones(md.numberofnodes,1);
+md.spcvz=NaN*ones(md.numberofnodes,1);
 
 %constrain flanks to 0 normal velocity
 pos=find(md.x==xmin | md.x==xmax);
-md.spcvelocity(pos,1)=1;
-md.spcvelocity(pos,3)=0;
+md.spcvx(pos)=0;
+md.spcvz(pos)=NaN;
 
 %constrain grounding line to 0 velocity
 pos=find(md.y==ymin);
-md.spcvelocity(pos,1:2)=1;
-md.spcvelocity(pos,3:4)=0;
+md.spcvx(pos)=0;
+md.spcvy(pos)=0;
 
 %icefront
Index: /issm/trunk/test/NightlyRun/test427.m
===================================================================
--- /issm/trunk/test/NightlyRun/test427.m	(revision 8823)
+++ /issm/trunk/test/NightlyRun/test427.m	(revision 8824)
@@ -8,8 +8,7 @@
 
 %constrain all velocities to 1 m/yr, in the y-direction
-md.spcvelocity(:,1:3)=1;
-md.spcvelocity(:,4)=0;
-md.spcvelocity(:,5)=1;
-md.spcvelocity(:,6)=0;
+md.spcvx(:)=0;
+md.spcvy(:)=1;
+md.spcvz(:)=0;
 
 %Dakota options
Index: /issm/trunk/test/NightlyRun/test428.m
===================================================================
--- /issm/trunk/test/NightlyRun/test428.m	(revision 8823)
+++ /issm/trunk/test/NightlyRun/test428.m	(revision 8824)
@@ -8,8 +8,7 @@
 
 %constrain all velocities to 1 m/yr, in the y-direction
-md.spcvelocity(:,1:3)=1;
-md.spcvelocity(:,4)=0;
-md.spcvelocity(:,5)=1;
-md.spcvelocity(:,6)=0;
+md.spcvx=zeros(md.numberofnodes,1);
+md.spcvy=ones(md.numberofnodes,1);
+md.spcvz=zeros(md.numberofnodes,1);
 
 %Dakota options
Index: /issm/trunk/test/NightlyRun/test434.m
===================================================================
--- /issm/trunk/test/NightlyRun/test434.m	(revision 8823)
+++ /issm/trunk/test/NightlyRun/test434.m	(revision 8824)
@@ -8,8 +8,7 @@
 
 %constrain all velocities to 1 m/yr, in the y-direction
-md.spcvelocity(:,1:3)=1;
-md.spcvelocity(:,4)=0;
-md.spcvelocity(:,5)=1;
-md.spcvelocity(:,6)=0;
+md.spcvx(:)=0;
+md.spcvy(:)=1;
+md.spcvz(:)=0;
 
 %Dakota options
Index: /issm/trunk/test/Par/ISMIPF.par
===================================================================
--- /issm/trunk/test/Par/ISMIPF.par	(revision 8823)
+++ /issm/trunk/test/Par/ISMIPF.par	(revision 8824)
@@ -23,5 +23,5 @@
 %Create node on boundary fist (because we cannot use mesh)
 md=SetIceSheetBC(md);
-md.spcvelocity(:,4)=100*ones(md.numberofnodes,1);
+md.spcvx=100*ones(md.numberofnodes,1);
 md.vx=zeros(md.numberofnodes,1);
 md.vy=zeros(md.numberofnodes,1);
Index: /issm/trunk/test/Par/Pig.par
===================================================================
--- /issm/trunk/test/Par/Pig.par	(revision 8823)
+++ /issm/trunk/test/Par/Pig.par	(revision 8824)
@@ -2,5 +2,5 @@
 
 %Geometry and observation
-load('./../Data/Pig.data','-mat');
+load('../Data/Pig.data','-mat');
 md.vx_obs   =InterpFromMeshToMesh2d(index,x,y,vx_obs,md.x,md.y);
 md.vy_obs   =InterpFromMeshToMesh2d(index,x,y,vy_obs,md.x,md.y);
Index: /issm/trunk/test/Par/RoundSheetEISMINT.par
===================================================================
--- /issm/trunk/test/Par/RoundSheetEISMINT.par	(revision 8823)
+++ /issm/trunk/test/Par/RoundSheetEISMINT.par	(revision 8824)
@@ -48,6 +48,7 @@
 md.x(pos)=0; md.y(pos)=0; %the closest node to the center is changed to be exactly at the center
 
-md.spcvelocity(pos,1:3)=1;
-md.spcvelocity(pos,4:6)=0;
+md.spcvx(pos)=0;
+md.spcvy(pos)=0;
+md.spcvz(pos)=0;
 
 %parallel options
Index: /issm/trunk/test/Par/RoundSheetShelf.par
===================================================================
--- /issm/trunk/test/Par/RoundSheetShelf.par	(revision 8823)
+++ /issm/trunk/test/Par/RoundSheetShelf.par	(revision 8824)
@@ -76,7 +76,12 @@
 
 %Deal with boundary conditions:
-md.spcvelocity=zeros(md.numberofnodes,6);
+md.spcvx=NaN*ones(md.numberofnodes,1);
+md.spcvy=NaN*ones(md.numberofnodes,1);
+md.spcvz=NaN*ones(md.numberofnodes,1);
+
 pos=find(md.x==0 & md.y==0);
-md.spcvelocity(pos,1:2)=1;
+md.spcvx(pos)=0;
+md.spcvy(pos)=0;
+
 pressureload=md.segments;
 pressureload=[pressureload WaterEnum*md.elementoniceshelf(pressureload(:,end))+AirEnum*md.elementonicesheet(pressureload(:,end))];
Index: /issm/trunk/test/Par/RoundSheetStaticEISMINT.par
===================================================================
--- /issm/trunk/test/Par/RoundSheetStaticEISMINT.par	(revision 8823)
+++ /issm/trunk/test/Par/RoundSheetStaticEISMINT.par	(revision 8824)
@@ -54,4 +54,5 @@
 md.x(pos)=0; md.y(pos)=0; %the closest node to the center is changed to be exactly at the center
 
-md.spcvelocity(pos,1:3)=1;
-md.spcvelocity(pos,4:6)=0;
+md.spcvx(pos)=0;
+md.spcvy(pos)=0;
+md.spcvz(pos)=0;