Index: /issm/trunk/src/m/solutions/cielo/prognostic.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/prognostic.m	(revision 676)
+++ /issm/trunk/src/m/solutions/cielo/prognostic.m	(revision 677)
@@ -12,6 +12,8 @@
 	md.analysis_type='prognostic'; m_p=CreateFemModel(md);
 
+	%Take only the first two dofs of m_p.parameters.u_g
+	u_g=m_p.parameters.u_g(dofsetgen([1,2],3,m_p.parameters.numberofnodes*3));
 	displaystring(md.debug,'\n%s',['depth averaging velocity...']);
-	u_g=VelocityDepthAverage(m_p.elements,m_p.nodes,m_p.loads,m_p.materials,m_p.parameters.u_g);
+	u_g=VelocityDepthAverage(m_p.elements,m_p.nodes,m_p.loads,m_p.materials,u_g);
 
 	% figure out number of dof: just for information purposes.
@@ -21,5 +23,5 @@
 	displaystring(md.debug,'\n%s',['setup inputs...']);
 	inputs=inputlist;
-	inputs=add(inputs,'velocity_average',u_g,'doublevec',3,m_p.parameters.numberofnodes);
+	inputs=add(inputs,'velocity_average',u_g,'doublevec',2,m_p.parameters.numberofnodes);
 	inputs=add(inputs,'thickness',m_p.parameters.h_g,'doublevec',1,m_p.parameters.numberofnodes);
 	inputs=add(inputs,'melting',m_p.parameters.m_g,'doublevec',1,m_p.parameters.numberofnodes);