Index: /issm/trunk-jpl/src/c/classes/FemModel.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/FemModel.cpp	(revision 16034)
+++ /issm/trunk-jpl/src/c/classes/FemModel.cpp	(revision 16035)
@@ -455,4 +455,5 @@
 	/*Get vector size and number of nodes*/
 	int numnodes            = nodes->NumberOfNodes(configuration_type);
+	int localnumnodes       = nodes->Size();
 	int numberofdofspernode = nodes->MaxNumDofs(configuration_type,GsetEnum);
 	int M                   = nodes->NumberOfDofs(configuration_type,set1enum);
@@ -516,5 +517,5 @@
 
 	/*OK now count number of dofs and flag each nodes for each node i*/
-	bool *flags                  = xNew<bool>(numnodes);
+	bool *flags                  = xNew<bool>(localnumnodes);
 	int  *d_connectivity         = xNewZeroInit<int>(numnodes);
 	int  *o_connectivity         = xNewZeroInit<int>(numnodes);
@@ -523,4 +524,5 @@
 
 	/*Create connectivity vector*/
+	MPI_Barrier(MPI_COMM_WORLD); double start_core2=MPI_Wtime( );
 	for(i=0;i<nodes->Size();i++){
 		Node* node=dynamic_cast<Node*>(nodes->GetObjectByOffset(i));
@@ -528,5 +530,5 @@
 
 			/*Reinitialize flags to 0*/
-			for(j=0;j<numnodes;j++) flags[j]=false;
+			for(j=0;j<localnumnodes;j++) flags[j]=false;
 
 			/*Loop over elements that hold node number i*/
@@ -567,4 +569,6 @@
 	xDelete<int>(head_l);
 	xDelete<int>(next_l);
+	MPI_Barrier(MPI_COMM_WORLD); double finish_core2=MPI_Wtime( );
+	_printf0_("\nTotal ellapsed time "<<finish_core2-start_core2<<" seconds  \n");
 
 	/*sum over all cpus*/
Index: /issm/trunk-jpl/src/c/classes/Loads/Numericalflux.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Loads/Numericalflux.cpp	(revision 16034)
+++ /issm/trunk-jpl/src/c/classes/Loads/Numericalflux.cpp	(revision 16035)
@@ -379,8 +379,8 @@
 	for(int i=0;i<this->GetNumberOfNodes();i++){
 
-		if(!flags[this->nodes[i]->Sid()]){
+		if(!flags[this->nodes[i]->Lid()]){
 
 			/*flag current node so that no other element processes it*/
-			flags[this->nodes[i]->Sid()]=true;
+			flags[this->nodes[i]->Lid()]=true;
 
 			/*if node is clone, we have an off-diagonal non-zero, else it is a diagonal non-zero*/
Index: /issm/trunk-jpl/src/c/classes/Loads/Pengrid.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Loads/Pengrid.cpp	(revision 16034)
+++ /issm/trunk-jpl/src/c/classes/Loads/Pengrid.cpp	(revision 16035)
@@ -298,8 +298,8 @@
 	int o_nz = 0;
 
-	if(!flags[this->node->Sid()]){
+	if(!flags[this->node->Lid()]){
 
 		/*flag current node so that no other element processes it*/
-		flags[this->node->Sid()]=true;
+		flags[this->node->Lid()]=true;
 
 		/*if node is clone, we have an off-diagonal non-zero, else it is a diagonal non-zero*/
Index: /issm/trunk-jpl/src/c/classes/Loads/Penpair.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Loads/Penpair.cpp	(revision 16034)
+++ /issm/trunk-jpl/src/c/classes/Loads/Penpair.cpp	(revision 16035)
@@ -217,8 +217,8 @@
 	for(int i=0;i<NUMVERTICES;i++){
 
-		if(!flags[this->nodes[i]->Sid()]){
+		if(!flags[this->nodes[i]->Lid()]){
 
 			/*flag current node so that no other element processes it*/
-			flags[this->nodes[i]->Sid()]=true;
+			flags[this->nodes[i]->Lid()]=true;
 
 			/*if node is clone, we have an off-diagonal non-zero, else it is a diagonal non-zero*/
Index: /issm/trunk-jpl/src/c/classes/Loads/Riftfront.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Loads/Riftfront.cpp	(revision 16034)
+++ /issm/trunk-jpl/src/c/classes/Loads/Riftfront.cpp	(revision 16035)
@@ -394,8 +394,8 @@
 	for(int i=0;i<NUMVERTICES;i++){
 
-		if(!flags[this->nodes[i]->Sid()]){
+		if(!flags[this->nodes[i]->Lid()]){
 
 			/*flag current node so that no other element processes it*/
-			flags[this->nodes[i]->Sid()]=true;
+			flags[this->nodes[i]->Lid()]=true;
 
 			/*if node is clone, we have an off-diagonal non-zero, else it is a diagonal non-zero*/