Changeset 23907

Timestamp:

05/03/19 14:08:30 (6 years ago)

Author:

Mathieu Morlighem

Message:

BUG: pressure nodes start now at iomodel->numberofvertices+iomodel->numberofedges+iomodel->numberofverticalfaces

File:

• issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp (modified) (1 diff)

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

@@ -159 +159 @@
                                                         if(iomodel->my_vertices[i]){
                                                                 if(IoCodeToEnumVertexEquation(reCast<int>(vertices_type[i]))==NoneApproximationEnum){
-                                                                        constraints->AddObject(new SpcStatic(count+1,iomodel->numberofvertices+iomodel->numberofedges+iomodel->numberoffaces+i+1,0,g*rho_ice*(surface[i]-z[i])/FSreconditioning,StressbalanceAnalysisEnum));
+                                                                        constraints->AddObject(new SpcStatic(count+1,iomodel->numberofvertices+iomodel->numberofedges+iomodel->numberofverticalfaces+i+1,0,g*rho_ice*(surface[i]-z[i])/FSreconditioning,StressbalanceAnalysisEnum));
                                                                         count++;
                                                                 }