Index: /issm/trunk-jpl/src/c/classes/FemModel.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/FemModel.cpp	(revision 24034)
+++ /issm/trunk-jpl/src/c/classes/FemModel.cpp	(revision 24035)
@@ -1314,9 +1314,4 @@
 void FemModel::GetLocalVectorWithClonesVertices(IssmDouble** plocal_vector,Vector<IssmDouble> *vector){/*{{{*/
 
-	/*recover my_rank:*/
-	ISSM_MPI_Status status;
-	int my_rank   = IssmComm::GetRank();
-	int num_procs = IssmComm::GetSize();
-
 	/*retrieve vertex info*/
 	int localsize         = this->vertices->NumberOfVerticesLocalAll();
@@ -1334,4 +1329,17 @@
 	xMemCpy<IssmDouble>(local_vector,local_vector_masters,localsize_masters);
 	xDelete<IssmDouble>(local_vector_masters);
+
+	/*Now send and receive vector for vertices on partition edge*/
+	SyncLocalVectorWithClonesVertices(local_vector);
+
+	/*Assign output pointer*/
+	*plocal_vector = local_vector;
+}/*}}}*/
+void FemModel::SyncLocalVectorWithClonesVertices(IssmDouble* local_vector){/*{{{*/
+
+	/*recover my_rank:*/
+	ISSM_MPI_Status status;
+	int my_rank   = IssmComm::GetRank();
+	int num_procs = IssmComm::GetSize();
 
 	/*Now send and receive vector for vertices on partition edge*/
@@ -1366,7 +1374,4 @@
 	}
 	xDelete<IssmDouble>(buffer);
-
-	/*Assign output pointer*/
-	*plocal_vector = local_vector;
 }/*}}}*/
 void FemModel::GetLocalVectorWithClonesNodes(IssmDouble** plocal_vector,Vector<IssmDouble> *vector){/*{{{*/
Index: /issm/trunk-jpl/src/c/classes/FemModel.h
===================================================================
--- /issm/trunk-jpl/src/c/classes/FemModel.h	(revision 24034)
+++ /issm/trunk-jpl/src/c/classes/FemModel.h	(revision 24035)
@@ -98,4 +98,5 @@
 		void GetLocalVectorWithClonesGset(IssmDouble** plocal_ug,Vector<IssmDouble> *ug);
 		void GetLocalVectorWithClonesVertices(IssmDouble** plocal_vector,Vector<IssmDouble> *vector);
+		void SyncLocalVectorWithClonesVertices(IssmDouble* local_vector);
 		void GetLocalVectorWithClonesNodes(IssmDouble** plocal_vector,Vector<IssmDouble> *vector);
 		void GroundedAreax(IssmDouble* pV, bool scaled);