Index: /issm/trunk/src/m/classes/public/loadresultsfromdisk.m
===================================================================
--- /issm/trunk/src/m/classes/public/loadresultsfromdisk.m	(revision 459)
+++ /issm/trunk/src/m/classes/public/loadresultsfromdisk.m	(revision 460)
@@ -30,14 +30,30 @@
 	u_g=results{4};
 
-	%Used to recover velocities
-	indx=1:2:gsize;
-	indy=2:2:gsize;
-	indx=indx(part);
-	indy=indy(part);
+	if strcmpi(md.type,'2d'),
+		%Used to recover velocities
+		indx=1:2:gsize;
+		indy=2:2:gsize;
+		indx=indx(part);
+		indy=indy(part);
 
-	%Recover velocity
-	md.vx=u_g(indx)*md.yts;
-	md.vy=u_g(indy)*md.yts;
-	md.vel=sqrt(md.vx.^2+md.vy.^2);
+		%Recover velocity
+		md.vx=u_g(indx)*md.yts;
+		md.vy=u_g(indy)*md.yts;
+		md.vel=sqrt(md.vx.^2+md.vy.^2);
+	else
+		%Used to recover velocities
+		indx=1:3:gsize;
+		indy=2:3:gsize;
+		indz=3:3:gsize;
+		indx=indx(part);
+		indy=indy(part);
+		indz=indz(part);
+
+		%Recover velocity
+		md.vx=u_g(indx)*md.yts;
+		md.vy=u_g(indy)*md.yts;
+		md.vz=u_g(indz)*md.yts;
+		md.vel=sqrt(md.vx.^2+md.vy.^2+md.vz.^2);
+	end
 
 elseif strcmpi(analysis_type,'control'),