Changeset 26079

Timestamp:

03/11/21 17:57:54 (4 years ago)

Author:

Eric.Larour

Message:

CHG: avoid looking for vertex on base if you are a Mesh 3D Surface domain type

File:

• issm/trunk-jpl/src/c/analyses/MasstransportAnalysis.cpp (modified) (3 diffs)

issm/trunk-jpl/src/c/analyses/MasstransportAnalysis.cpp

@@ -62 +62 @@
         IssmDouble *nodeonbase=NULL;
         iomodel->FetchData(&vertex_pairing,&numvertex_pairing,NULL,"md.masstransport.vertex_pairing");
-        if(iomodel->domaintype!=Domain2DhorizontalEnum) iomodel->FetchData(&nodeonbase,NULL,NULL,"md.mesh.vertexonbase");
+        if(iomodel->domaintype!=Domain2DhorizontalEnum && iomodel->domaintype!=Domain3DsurfaceEnum) iomodel->FetchData(&nodeonbase,NULL,NULL,"md.mesh.vertexonbase");
 
         for(int i=0;i<numvertex_pairing;i++){
@@ -100 +100 @@
 
         /*Create Nodes either DG or CG depending on stabilization*/
-        if(iomodel->domaintype!=Domain2DhorizontalEnum) iomodel->FetchData(2,"md.mesh.vertexonbase","md.mesh.vertexonsurface");
+        if(iomodel->domaintype!=Domain2DhorizontalEnum && iomodel->domaintype!=Domain3DsurfaceEnum) iomodel->FetchData(2,"md.mesh.vertexonbase","md.mesh.vertexonsurface");
         if(stabilization!=3){
                 ::CreateNodes(nodes,iomodel,MasstransportAnalysisEnum,FINITEELEMENT,isamr);
@@ -221 +221 @@
         }
 
-        if(iomodel->domaintype!=Domain2DhorizontalEnum){
+        if(iomodel->domaintype!=Domain2DhorizontalEnum && iomodel->domaintype!=Domain3DsurfaceEnum){
                 iomodel->FetchDataToInput(inputs,elements,"md.mesh.vertexonbase",MeshVertexonbaseEnum);
                 iomodel->FetchDataToInput(inputs,elements,"md.mesh.vertexonsurface",MeshVertexonsurfaceEnum);