Index: /issm/trunk/src/m/solutions/cielo/diagnostic.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/diagnostic.m	(revision 128)
+++ /issm/trunk/src/m/solutions/cielo/diagnostic.m	(revision 129)
@@ -1,3 +1,7 @@
 function md=diagnostic(md)
+%DIAGNOSTIC - diagnostic solution sequence
+%
+%   Usage:
+%      md=diagnostic(md)
 
 	%timing
@@ -38,9 +42,8 @@
 		md.pressure=md.rho_ice*md.g*(md.surface-md.z);
 	else
-		
 		%load results onto model
 		md.vx=u_g(1:2:end)*md.yts;
 		md.vy=u_g(2:2:end)*md.yts;
-		md.vz=zeros(m_dv.nodesets.gsize,1);
+		md.vz=zeros(md.numberofgrids,1);
 		md.vel=sqrt(md.vx.^2+md.vy.^2+md.vz.^2);