Changeset 1703
- Timestamp:
- 08/14/09 11:49:39 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/m/solutions/cielo/ControlInitialization.m
r1654 r1703 38 38 %horizontal velocities 39 39 displaystring(debug,'\n%s',['computing horizontal velocities...']); 40 u_g=diagnostic_core_nonlinear(m_dh,inputs, 'diagnostic','horiz');40 u_g=diagnostic_core_nonlinear(m_dh,inputs,DiagnosticAnalysisEnum(),HorizAnalysisEnum()); 41 41 displaystring(debug,'\n%s',['extruding horizontal velocities...']); 42 42 u_g_horiz=FieldExtrude(m_dh.elements,m_dh.nodes,m_dh.loads,m_dh.materials,u_g,'velocity',1); … … 45 45 displaystring(debug,'\n%s',['computing vertical velocities...']); 46 46 inputs=add(inputs,'velocity',u_g_horiz,'doublevec',m_dh.parameters.numberofdofspernode,m_dh.parameters.numberofnodes); 47 u_g_vert=diagnostic_core_linear(m_dv,inputs, 'diagnostic','vert');47 u_g_vert=diagnostic_core_linear(m_dv,inputs,DiagnosticAnalysisEnum(),VertAnalysisEnum()); 48 48 49 49 %create 3d u_g … … 70 70 %Compute Stokes velocities once to have a reasonably good ug in input 71 71 displaystring(debug,'\n%s',['computing stokes velocities and pressure ...']); 72 u_g=diagnostic_core_nonlinear(m_ds,inputs, 'diagnostic','stokes');72 u_g=diagnostic_core_nonlinear(m_ds,inputs,DiagnosticAnalysisEnum(),StokesAnalysisEnum()); 73 73 inputs=add(inputs,'velocity',u_g,'doublevec',4,m_ds.parameters.numberofnodes); 74 74
Note:
See TracChangeset
for help on using the changeset viewer.