Changeset 25219

Timestamp:

07/06/20 19:32:50 (5 years ago)

Author:

Eric.Larour

Message:

CHG: fixed bug in the mme processing.

File:

• issm/trunk-jpl/src/c/analyses/SealevelriseAnalysis.cpp (modified) (2 diffs)

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

@@ -164 +164 @@
                                 _assert_(iomodel->elements);
                                 for(int k=0;k<numvertices;k++){
-                                        vertexsids[k] =reCast<int>(iomodel->elements[numvertices*element->Sid()+k]); //ids for vertices are in the elements array from Matlab
-                                        vertexlids[k]=iomodel->my_vertices_lids[vertexsids[k]-1];
+                                        vertexsids[k] =reCast<int>(iomodel->elements[numvertices*element->Sid()+k]-1); //ids for vertices are in the elements array from Matlab
+                                        vertexlids[k]=iomodel->my_vertices_lids[vertexsids[k]];
                                 }
 
-                                //element->GetVerticesLidList(vertexlids);
-                                //element->GetVerticesSidList(vertexsids);
                                 IssmDouble* values=xNew<IssmDouble>(numvertices);
 
@@ -226 +224 @@
                                 _assert_(iomodel->elements);
                                 for(int k=0;k<numvertices;k++){
-                                        vertexsids[k] =reCast<int>(iomodel->elements[numvertices*element->Sid()+k]); //ids for vertices are in the elements array from Matlab
-                                        vertexlids[k]=iomodel->my_vertices_lids[vertexsids[k]-1];
+                                        vertexsids[k] =reCast<int>(iomodel->elements[numvertices*element->Sid()+k]-1); //ids for vertices are in the elements array from Matlab
+                                        vertexlids[k]=iomodel->my_vertices_lids[vertexsids[k]];
                                 }
-                                //element->GetVerticesLidList(vertexlids);
-                                //element->GetVerticesSidList(vertexsids);
                                 
                                 IssmDouble* values=xNew<IssmDouble>(numvertices);