- Timestamp:
- 04/08/14 20:39:20 (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/analyses/SmoothedSurfaceSlopeYAnalysis.cpp
r17674 r17686 6 6 7 7 /*Model processing*/ 8 int SmoothedSurfaceSlopeYAnalysis::DofsPerNode(int** doflist,int meshxdim,int approximation){/*{{{*/8 int SmoothedSurfaceSlopeYAnalysis::DofsPerNode(int** doflist,int domaintype,int approximation){/*{{{*/ 9 9 return 1; 10 10 }/*}}}*/ … … 26 26 iomodel->FetchDataToInput(elements,BaseEnum); 27 27 iomodel->FetchDataToInput(elements,MaskIceLevelsetEnum); 28 if(iomodel-> meshxdim==Mesh3DEnum){28 if(iomodel->domaintype==Mesh3DEnum){ 29 29 iomodel->FetchDataToInput(elements,MeshElementonbaseEnum); 30 30 iomodel->FetchDataToInput(elements,MeshElementonsurfaceEnum); 31 31 } 32 if(iomodel-> meshxdim==Mesh2DverticalEnum){32 if(iomodel->domaintype==Mesh2DverticalEnum){ 33 33 iomodel->FetchDataToInput(elements,MeshVertexonbaseEnum); 34 34 } 35 35 }/*}}}*/ 36 36 void SmoothedSurfaceSlopeYAnalysis::CreateNodes(Nodes* nodes,IoModel* iomodel){/*{{{*/ 37 if(iomodel-> meshxdim==Mesh3DEnum) iomodel->FetchData(1,MeshVertexonbaseEnum);37 if(iomodel->domaintype==Mesh3DEnum) iomodel->FetchData(1,MeshVertexonbaseEnum); 38 38 ::CreateNodes(nodes,iomodel,SmoothedSurfaceSlopeYAnalysisEnum,P1Enum); 39 if(iomodel-> meshxdim==Mesh3DEnum) iomodel->DeleteData(1,MeshVertexonbaseEnum);39 if(iomodel->domaintype==Mesh3DEnum) iomodel->DeleteData(1,MeshVertexonbaseEnum); 40 40 }/*}}}*/ 41 41 void SmoothedSurfaceSlopeYAnalysis::CreateConstraints(Constraints* constraints,IoModel* iomodel){/*{{{*/ … … 105 105 106 106 /*Intermediaries*/ 107 int meshxdim;107 int domaintype; 108 108 Element* basalelement; 109 109 110 110 /*Get basal element*/ 111 element->FindParam(& meshxdim,MeshXDimEnum);112 switch( meshxdim){111 element->FindParam(&domaintype,DomainTypeEnum); 112 switch(domaintype){ 113 113 case Mesh2DhorizontalEnum: 114 114 basalelement = element; … … 118 118 basalelement = element->SpawnBasalElement(); 119 119 break; 120 default: _error_("mesh "<<EnumToStringx( meshxdim)<<" not supported yet");120 default: _error_("mesh "<<EnumToStringx(domaintype)<<" not supported yet"); 121 121 } 122 122 … … 169 169 xDelete<IssmDouble>(basis); 170 170 delete gauss; 171 if( meshxdim!=Mesh2DhorizontalEnum){basalelement->DeleteMaterials(); delete basalelement;};171 if(domaintype!=Mesh2DhorizontalEnum){basalelement->DeleteMaterials(); delete basalelement;}; 172 172 return pe; 173 173 }/*}}}*/
Note:
See TracChangeset
for help on using the changeset viewer.