Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/ElementsAndVerticesPartitioning.cpp
===================================================================
--- /issm/trunk-jpl/src/c/modules/ModelProcessorx/ElementsAndVerticesPartitioning.cpp	(revision 23503)
+++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/ElementsAndVerticesPartitioning.cpp	(revision 23504)
@@ -17,16 +17,8 @@
 void  ElementsAndVerticesPartitioning(bool** pmy_elements, int** pmy_vertices, IoModel* iomodel){
 
-	int i,j;
-
-	const int RIFTINFOSIZE = 12;
 	int numberofelements2d;
 	int numberofvertices2d;
 	int numlayers;
-	int numrifts;
 	int numvertex_pairing;
-
-	/*output: */
-	bool *my_elements = NULL;
-	int  *my_vertices = NULL;
 
 	/*intermediary: */
@@ -36,6 +28,4 @@
 	int         el1,el2;
 	int        *elements2d     = NULL;
-	int        *vertex_pairing = NULL;
-	IssmDouble *riftinfo       = NULL;
 
 	/*Get my_rank:*/
@@ -43,9 +33,6 @@
 	int num_procs = IssmComm::GetSize();
 
-	/*Fetch parameters: */
-	iomodel->FindConstant(&numrifts,"md.rifts.numrifts");
-
 	/*First, check that partitioning has not yet been carryed out. Just check whether my_elements pointers is not already assigned a value: */
-	if(*pmy_elements)return;
+	if(*pmy_elements) return;
 
 	/*Number of vertices per elements, needed to correctly retrieve data: */
@@ -75,13 +62,17 @@
 	}
 
+	/*Partition and free resouces*/
 	MeshPartitionx(&epart,&npart,iomodel->numberofelements,iomodel->numberofvertices,iomodel->elements,numberofelements2d,numberofvertices2d,elements2d,numlayers,elements_width,iomodel->meshelementtype,num_procs);
-
-	/*Free elements2d: */
 	xDelete<int>(elements2d);
+	xDelete<int>(npart);
 
 	/*Deal with rifts, they have to be included into one partition only, not several: */
+	int numrifts;
+	iomodel->FindConstant(&numrifts,"md.rifts.numrifts");
 	if(numrifts){
+		IssmDouble *riftinfo = NULL;
 		iomodel->FetchData(&riftinfo,&numrifts,NULL,"md.rifts.riftstruct");
-		for(i=0;i<numrifts;i++){
+		for(int i=0;i<numrifts;i++){
+			const int RIFTINFOSIZE = 12;
 			el1=reCast<int>(*(riftinfo+RIFTINFOSIZE*i+2))-1; //matlab indexing to c indexing
 			el2=reCast<int>(*(riftinfo+RIFTINFOSIZE*i+3))-1; //matlab indexing to c indexing
@@ -91,10 +82,10 @@
 	}
 
-	/*Used later on: */
-	my_vertices=xNewZeroInit<int>(iomodel->numberofvertices);
-	my_elements=xNewZeroInit<bool>(iomodel->numberofelements);
+	/*Create my_vertices and my_elements, used by each partition */
+	bool *my_elements = xNewZeroInit<bool>(iomodel->numberofelements);
+	int  *my_vertices = xNewZeroInit<int>(iomodel->numberofvertices);
 
 	/*Start figuring out, out of the partition, which elements belong to this cpu: */
-	for(i=0;i<iomodel->numberofelements;i++){
+	for(int i=0;i<iomodel->numberofelements;i++){
 
 		/*!All elements have been partitioned above, only deal with elements for this cpu: */
@@ -105,5 +96,5 @@
 			 into the vertices coordinates. If we start plugging 1 into my_vertices for each index[n][i] (i=0:2), then my_vertices 
 			 will hold which vertices belong to this partition*/
-			for(j=0;j<elements_width;j++){
+			for(int j=0;j<elements_width;j++){
 				my_vertices[iomodel->elements[elements_width*i+j]-1]=1;
 			}
@@ -114,6 +105,7 @@
 	 * penpair has 2 nodes that are poointing toward 2 vertices.
 	 * The 2 vertices must be in the same cpu as the penpair*/
+	int *vertex_pairing = NULL;
 	iomodel->FetchData(&vertex_pairing,&numvertex_pairing,NULL,"md.stressbalance.vertex_pairing");
-	for(i=0;i<numvertex_pairing;i++){
+	for(int i=0;i<numvertex_pairing;i++){
 		if(my_vertices[vertex_pairing[2*i+0]-1] && !my_vertices[vertex_pairing[2*i+1]-1]){
 			my_vertices[vertex_pairing[2*i+1]-1]=2; //to know that these elements are not on the partition
@@ -122,5 +114,5 @@
 	xDelete<int>(vertex_pairing);
 	iomodel->FetchData(&vertex_pairing,&numvertex_pairing,NULL,"md.masstransport.vertex_pairing");
-	for(i=0;i<numvertex_pairing;i++){
+	for(int i=0;i<numvertex_pairing;i++){
 		if(my_vertices[vertex_pairing[2*i+0]-1] && !my_vertices[vertex_pairing[2*i+1]-1]){
 			my_vertices[vertex_pairing[2*i+1]-1]=2; //to know that these elements are not on the partition
@@ -129,9 +121,6 @@
 	xDelete<int>(vertex_pairing);
 
-	/*Free ressources:*/
-	xDelete<int>(npart);
+	/*cleanup and assign output pointer*/
 	xDelete<int>(epart);
-
-	/*Assign output pointers:*/
 	*pmy_elements=my_elements;
 	*pmy_vertices=my_vertices;