Index: /issm/trunk/src/m/solutions/cielo/diagnostic.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/diagnostic.m	(revision 133)
+++ /issm/trunk/src/m/solutions/cielo/diagnostic.m	(revision 134)
@@ -22,5 +22,17 @@
 	disp(sprintf('\n%s',['computing horizontal velocities...']));
 
-	u_g=diagnostic_core_nonlinear(m_dh,{});
+	%plug existing velocity in inputs
+	if ~isnan(md.vx) & ~isnan(md.vy),
+		indx=m_dh.part(1:2:end);
+		indy=m_dh.part(2:2:end);
+		velocity=zeros(m_dh.nodesets.gsize,1);
+		velocity(indx)=md.vx/md.yts;
+		velocity(indy)=md.vy/md.yts;
+		inputs=struct('velocity',velocity);
+	else
+		inputs={};
+	end
+
+	u_g=diagnostic_core_nonlinear(m_dh,inputs);
 
 	if strcmpi(md.type,'3d'),