Changeset 16291 for issm/trunk-jpl/src/c/modules/ModelProcessorx/Stressbalance/CreateConstraintsStressbalance.cpp
- Timestamp:
- 10/03/13 14:53:55 (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/modules/ModelProcessorx/Stressbalance/CreateConstraintsStressbalance.cpp
r16218 r16291 161 161 iomodel->FetchData(&spcvz,&Mz,&Nz,StressbalanceSpcvzEnum); 162 162 iomodel->FetchData(&nodeonSSA,NULL,NULL,FlowequationBorderSSAEnum); 163 if(iomodel-> dim==3)iomodel->FetchData(&nodeonHO,NULL,NULL,FlowequationBorderHOEnum);164 if(iomodel-> dim==3)iomodel->FetchData(&nodeonFS,NULL,NULL,FlowequationBorderFSEnum);165 if(iomodel-> dim==3)iomodel->FetchData(&nodeonbed,NULL,NULL,MeshVertexonbedEnum);166 if(iomodel-> dim==3)iomodel->FetchData(&groundedice_ls,NULL,NULL,MaskGroundediceLevelsetEnum);163 if(iomodel->meshtype==Mesh3DEnum)iomodel->FetchData(&nodeonHO,NULL,NULL,FlowequationBorderHOEnum); 164 if(iomodel->meshtype==Mesh3DEnum)iomodel->FetchData(&nodeonFS,NULL,NULL,FlowequationBorderFSEnum); 165 if(iomodel->meshtype==Mesh3DEnum)iomodel->FetchData(&nodeonbed,NULL,NULL,MeshVertexonbedEnum); 166 if(iomodel->meshtype==Mesh3DEnum)iomodel->FetchData(&groundedice_ls,NULL,NULL,MaskGroundediceLevelsetEnum); 167 167 iomodel->FetchData(&vertices_type,NULL,NULL,FlowequationVertexEquationEnum); 168 168 iomodel->FetchData(&surface,NULL,NULL,SurfaceEnum); … … 385 385 386 386 /*Constraint at the bedrock interface (v.n = vz = 0) (Coordinates will be updated according to the bed slope)*/ 387 if (iomodel-> dim==3) if(nodeonbed[i]>0. && groundedice_ls[i]>=0. && nodeonFS[i]>0.){387 if (iomodel->meshtype==Mesh3DEnum) if(nodeonbed[i]>0. && groundedice_ls[i]>=0. && nodeonFS[i]>0.){ 388 388 switch(reCast<int,IssmDouble>(vertices_type[i])){ 389 389 case SSAFSApproximationEnum: … … 410 410 iomodel->DeleteData(spcvz,StressbalanceSpcvzEnum); 411 411 iomodel->DeleteData(nodeonSSA,FlowequationBorderSSAEnum); 412 if(iomodel-> dim==3)iomodel->DeleteData(nodeonHO,FlowequationBorderHOEnum);413 if(iomodel-> dim==3)iomodel->DeleteData(nodeonFS,FlowequationBorderFSEnum);414 if(iomodel-> dim==3)iomodel->DeleteData(nodeonbed,MeshVertexonbedEnum);415 if(iomodel-> dim==3)iomodel->DeleteData(groundedice_ls,MaskGroundediceLevelsetEnum);412 if(iomodel->meshtype==Mesh3DEnum)iomodel->DeleteData(nodeonHO,FlowequationBorderHOEnum); 413 if(iomodel->meshtype==Mesh3DEnum)iomodel->DeleteData(nodeonFS,FlowequationBorderFSEnum); 414 if(iomodel->meshtype==Mesh3DEnum)iomodel->DeleteData(nodeonbed,MeshVertexonbedEnum); 415 if(iomodel->meshtype==Mesh3DEnum)iomodel->DeleteData(groundedice_ls,MaskGroundediceLevelsetEnum); 416 416 iomodel->DeleteData(vertices_type,FlowequationVertexEquationEnum); 417 417 iomodel->DeleteData(surface,SurfaceEnum);
Note:
See TracChangeset
for help on using the changeset viewer.