Index: /issm/trunk/src/c/ModelProcessorx/DiagnosticStokes/CreateLoadsDiagnosticStokes.cpp
===================================================================
--- /issm/trunk/src/c/ModelProcessorx/DiagnosticStokes/CreateLoadsDiagnosticStokes.cpp	(revision 3623)
+++ /issm/trunk/src/c/ModelProcessorx/DiagnosticStokes/CreateLoadsDiagnosticStokes.cpp	(revision 3624)
@@ -68,13 +68,13 @@
 	IoModelFetchData(&iomodel->gridonicesheet,NULL,NULL,iomodel_handle,"gridonicesheet");
 	IoModelFetchData(&iomodel->gridonstokes,NULL,NULL,iomodel_handle,"gridonstokes");
+	iomodel->CreateSingleNodeToElementConnectivity;
 	
 	for (i=0;i<iomodel->numberofvertices;i++){
-		/*keep only this partition's nodes:*/
+		
 		if((iomodel->my_vertices[i]==1)){
 
 			if ((iomodel->gridonbed[i]) && (iomodel->gridonicesheet[i]) && (iomodel->gridonstokes[i])){
-
-				loads->AddObject(new Pengrid(count+1, i+1,iomodel->numberofelements+1,1, 0, iomodel->penalty_offset,1,iomodel->stabilize_constraints));
-				count++;
+				
+				loads->AddObject(new Pengrid(count,i,iomodel));
 			}
 		} //if((iomodel->my_vertices[i]==1))
@@ -84,4 +84,5 @@
 	xfree((void**)&iomodel->gridonstokes);
 	xfree((void**)&iomodel->gridonicesheet);
+	xfree((void**)&iomodel->singlenodetoelementconnectivity);
 
 	/*All our datasets are already order by ids. Set presort flag so that later on, when sorting is requested on these 
Index: /issm/trunk/src/c/ModelProcessorx/Melting/CreateLoadsMelting.cpp
===================================================================
--- /issm/trunk/src/c/ModelProcessorx/Melting/CreateLoadsMelting.cpp	(revision 3623)
+++ /issm/trunk/src/c/ModelProcessorx/Melting/CreateLoadsMelting.cpp	(revision 3624)
@@ -13,5 +13,4 @@
 
 	int i;
-	int count=0;
 
 	DataSet*    loads    = NULL;
@@ -25,18 +24,18 @@
 	//create penalties for grids: no grid can have a temperature over the melting point
 	IoModelFetchData(&iomodel->gridonbed,NULL,NULL,iomodel_handle,"gridonbed");
+	iomodel->CreateSingleNodeToElementConnectivity;
 
-	count=1; //matlab indexing
 	for (i=0;i<iomodel->numberofvertices;i++){
-		/*keep only this partition's nodes:*/
+		
 		if((iomodel->my_vertices[i]==1)){
 
 			if (iomodel->gridonbed[i]){ 
-
-				loads->AddObject( new Pengrid(count+1, i+1,iomodel->numberofelements+1, 1, 0, iomodel->penalty_offset, 1, iomodel->stabilize_constraints));
-				count++;
+				
+				loads->AddObject(new Pengrid(i+1,i,iomodel));
 			}
-		} //if((iomodel->my_vertices[i]==1))
+		}
 	}
 	xfree((void**)&iomodel->gridonbed);
+	xfree((void**)&iomodel->singlenodetoelementconnectivity);
 
 	/*All our datasets are already order by ids. Set presort flag so that later on, when sorting is requested on these 
Index: /issm/trunk/src/c/ModelProcessorx/Thermal/CreateLoadsThermal.cpp
===================================================================
--- /issm/trunk/src/c/ModelProcessorx/Thermal/CreateLoadsThermal.cpp	(revision 3623)
+++ /issm/trunk/src/c/ModelProcessorx/Thermal/CreateLoadsThermal.cpp	(revision 3624)
@@ -13,5 +13,4 @@
 
 	int i;
-	int count=0;
 
 	DataSet*    loads    = NULL;
@@ -26,7 +25,8 @@
 	//create penalties for grids: no grid can have a temperature over the melting point
 	IoModelFetchData(&iomodel->spctemperature,NULL,NULL,iomodel_handle,"spctemperature");
+	iomodel->CreateSingleNodeToElementConnectivity;
 
-	count=1; //matlab indexing
 	for (i=0;i<iomodel->numberofvertices;i++){
+	
 		/*keep only this partition's nodes:*/
 		if((iomodel->my_vertices[i]==1)){
@@ -34,10 +34,11 @@
 			if (!iomodel->spctemperature[2*i]){ //No penalty applied on spc grids!
 
-				loads->AddObject( new Pengrid(count+1, i+1,iomodel->numberofelements+1,1, 0, iomodel->penalty_offset,1,iomodel->stabilize_constraints));
-				count++;
+				loads->AddObject(new Pengrid(i+1,i,iomodel));
+						
 			}
-		} //if((iomodel->my_vertices[i]==1))
+		}
 	}
 	xfree((void**)&iomodel->spctemperature);
+	xfree((void**)&iomodel->singlenodetoelementconnectivity);
 
 	/*All our datasets are already order by ids. Set presort flag so that later on, when sorting is requested on these 
@@ -51,2 +52,4 @@
 
 }
+				
+