Changeset 518

Timestamp:

05/20/09 09:50:53 (16 years ago)

Author:

seroussi

Message:

initialise y_g with zeros

File:

• issm/trunk/src/m/solutions/cielo/diagnostic_core.m (modified) (1 diff)

issm/trunk/src/m/solutions/cielo/diagnostic_core.m

@@ -91 +91 @@
 
                 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);