Changeset 17681
- Timestamp:
- 04/08/14 17:34:23 (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/classes/Elements/Tria.cpp
r17680 r17681 165 165 IssmDouble sigma_yz[NUMVERTICES]={0,0,0}; 166 166 GaussTria* gauss=NULL; 167 int mesh type,dim=2;167 int meshxdim,dim=2; 168 168 169 169 /* Get node coordinates and dof list: */ … … 171 171 172 172 /*Retrieve all inputs we will be needing: */ 173 this->FindParam(&mesh type,MeshTypeEnum);174 if(mesh type==Mesh2DhorizontalEnum) _error_("stress tensor calculation not supported for mesh of type " <<EnumToStringx(meshtype)<<", extrude mesh or call ComputeDeviatoricStressTensor");173 this->FindParam(&meshxdim,MeshTypeEnum); 174 if(meshxdim==Mesh2DhorizontalEnum) _error_("stress tensor calculation not supported for mesh of type " <<EnumToStringx(meshxdim)<<", extrude mesh or call ComputeDeviatoricStressTensor"); 175 175 Input* pressure_input=inputs->GetInput(PressureEnum); _assert_(pressure_input); 176 176 Input* vx_input=inputs->GetInput(VxEnum); _assert_(vx_input); … … 218 218 IssmDouble tau_yz[NUMVERTICES]={0,0,0}; 219 219 GaussTria* gauss=NULL; 220 int mesh type,dim=2;220 int meshxdim,dim=2; 221 221 222 222 /* Get node coordinates and dof list: */ … … 224 224 225 225 /*Retrieve all inputs we will be needing: */ 226 this->FindParam(&mesh type,MeshTypeEnum);227 if(mesh type!=Mesh2DhorizontalEnum) _error_("deviatoric stress tensor calculation not implemented for mesh of type " <<EnumToStringx(meshtype));226 this->FindParam(&meshxdim,MeshTypeEnum); 227 if(meshxdim!=Mesh2DhorizontalEnum) _error_("deviatoric stress tensor calculation not implemented for mesh of type " <<EnumToStringx(meshxdim)); 228 228 Input* vx_input=inputs->GetInput(VxEnum); _assert_(vx_input); 229 229 Input* vy_input=inputs->GetInput(VyEnum); _assert_(vy_input);
Note:
See TracChangeset
for help on using the changeset viewer.