Index: /issm/trunk-jpl/src/c/classes/FemModel.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/FemModel.cpp	(revision 21390)
+++ /issm/trunk-jpl/src/c/classes/FemModel.cpp	(revision 21391)
@@ -2497,5 +2497,5 @@
 	IssmDouble	load_love_k2 = -0.30922675; //degree 2 load Love number 
 	IssmDouble	m1, m2, m3; 
-	IssmDouble	lati, longi, value; 
+	IssmDouble	lati, longi, radi, value; 
 
 	/*Serialize vectors from previous iteration:*/
@@ -2547,14 +2547,16 @@
 	 */
 	for(int i=0;i<vertices->Size();i++){
+		int sid;
 		//Vertex* vertex=(Vertex*)vertices->GetObjectByOffset(i);
 		Vertex* vertex=xDynamicCast<Vertex*>(vertices->GetObjectByOffset(i));
-		
-		lati=latitude[i]/180*PI;	longi=longitude[i]/180*PI;
+		sid=vertex->Sid();
+
+		lati=latitude[sid]/180*PI;	longi=longitude[sid]/180*PI; radi=radius[sid];
 
 		/*only first order terms are considered now: */ 
-		value=((1.0+tide_love_k-tide_love_h)/9.81)*pow(omega*radius[i],2.0)*
+		value=((1.0+tide_love_k-tide_love_h)/9.81)*pow(omega*radi,2.0)*
 						(-m3/6.0 + 0.5*m3*cos(2.0*lati) - 0.5*sin(2.*lati)*(m1*cos(longi)+m2*sin(longi))); 
 	
-		pSgo_rot->SetValue(vertex->Sid(),value,INS_VAL); //INS_VAL ensures that you don't add several times
+		pSgo_rot->SetValue(sid,value,INS_VAL); //INS_VAL ensures that you don't add several times
 	}