Changeset 677

Timestamp:

06/01/09 16:09:27 (16 years ago)

Author:

seroussi

Message:

velocity average on horizontal velocity only

File:

• issm/trunk/src/m/solutions/cielo/prognostic.m (modified) (2 diffs)

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

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