Index: /issm/trunk-jpl/src/c/analyses/SealevelriseAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/SealevelriseAnalysis.cpp	(revision 25218)
+++ /issm/trunk-jpl/src/c/analyses/SealevelriseAnalysis.cpp	(revision 25219)
@@ -164,10 +164,8 @@
 				_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,9 +224,7 @@
 				_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);