Index: /issm/trunk/src/m/solutions/cielo/diagnostic_core.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/diagnostic_core.m	(revision 517)
+++ /issm/trunk/src/m/solutions/cielo/diagnostic_core.m	(revision 518)
@@ -91,4 +91,5 @@
 
 		displaystring(debug,'\n%s',['update boundary conditions for stokes using velocities previously computed...']);
+		m_ds.y_g=zeros(m_ds.nodesets.gsize,1);
 		m_ds.y_g(dofsetgen([1,2,3],4,m_ds.nodesets.gsize))=u_g;
 		[m_ds.ys m_ds.ys0]=Reducevectorgtos(m_ds.y_g,m_ds.nodesets);