Index: /issm/trunk/src/c/ModelProcessorx/Balancedthickness2/CreateElementsNodesAndMaterialsBalancedthickness2.cpp
===================================================================
--- /issm/trunk/src/c/ModelProcessorx/Balancedthickness2/CreateElementsNodesAndMaterialsBalancedthickness2.cpp	(revision 3597)
+++ /issm/trunk/src/c/ModelProcessorx/Balancedthickness2/CreateElementsNodesAndMaterialsBalancedthickness2.cpp	(revision 3598)
@@ -135,4 +135,5 @@
 	xfree((void**)&iomodel->gridonicesheet);
 	xfree((void**)&iomodel->gridoniceshelf);
+	xfree((void**)&iomodel->elements);
 
 	/*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/io/ParameterInputsInit.cpp
===================================================================
--- /issm/trunk/src/c/io/ParameterInputsInit.cpp	(revision 3597)
+++ /issm/trunk/src/c/io/ParameterInputsInit.cpp	(revision 3598)
@@ -92,3 +92,3 @@
 void ParameterInputs::Init( void* vinput_handle){
 } 
-#endif //#ifdef _SERIAL-
+#endif //#ifdef _SERIAL_
Index: /issm/trunk/src/c/parallel/balancedthickness.cpp
===================================================================
--- /issm/trunk/src/c/parallel/balancedthickness.cpp	(revision 3597)
+++ /issm/trunk/src/c/parallel/balancedthickness.cpp	(revision 3598)
@@ -25,5 +25,4 @@
 	Model* model=NULL;
 
-	Vec     u_g=NULL;
 	double* u_g_serial=NULL;
 	double* melting_g=NULL;
@@ -67,5 +66,4 @@
 	/*Initialize model structure: */
 	MPI_Barrier(MPI_COMM_WORLD); start_init=MPI_Wtime();
-	model=new Model();
 
 	/*Open handle to data on disk: */
@@ -138,4 +136,7 @@
 
 	/*Free ressources:*/
+	xfree((void**)&u_g_serial);
+	xfree((void**)&melting_g);
+	xfree((void**)&accumulation_g);
 	delete processedresults;
 	delete results;
Index: /issm/trunk/src/c/parallel/balancedthickness2.cpp
===================================================================
--- /issm/trunk/src/c/parallel/balancedthickness2.cpp	(revision 3597)
+++ /issm/trunk/src/c/parallel/balancedthickness2.cpp	(revision 3598)
@@ -69,5 +69,4 @@
 	/*Initialize model structure: */
 	MPI_Barrier(MPI_COMM_WORLD); start_init=MPI_Wtime();
-	model=new Model();
 
 	/*Open handle to data on disk: */
Index: /issm/trunk/src/c/parallel/balancedvelocities.cpp
===================================================================
--- /issm/trunk/src/c/parallel/balancedvelocities.cpp	(revision 3597)
+++ /issm/trunk/src/c/parallel/balancedvelocities.cpp	(revision 3598)
@@ -64,5 +64,4 @@
 	/*Initialize model structure: */
 	MPI_Barrier(MPI_COMM_WORLD); start_init=MPI_Wtime();
-	model=new Model();
 
 	/*Open handle to data on disk: */
@@ -135,4 +134,7 @@
 
 	/*Free ressources:*/
+	xfree((void**)&u_g_serial);
+	xfree((void**)&melting_g);
+	xfree((void**)&accumulation_g);
 	delete processedresults;
 	delete results;
Index: /issm/trunk/src/c/parallel/diagnostic_core.cpp
===================================================================
--- /issm/trunk/src/c/parallel/diagnostic_core.cpp	(revision 3597)
+++ /issm/trunk/src/c/parallel/diagnostic_core.cpp	(revision 3598)
@@ -81,7 +81,5 @@
 	fem_ds->FindParam(&numberofdofspernode_ds,"numberofdofspernode");
 
-	//for qmu analysis, be sure the velocity input we are starting from  is the one in the parameters, otherwise, 
-	//we will create sensitivities just because we are starting from different velocities each time (and if tolerance is 
-	//not set tight enough): */
+	//for qmu analysis, be sure the velocity input we are starting from  is the one in the parameters: */
 	if(qmu_analysis){
 		model->FindParam(&u_g_initial,NULL,NULL,"u_g",DiagnosticAnalysisEnum,HorizAnalysisEnum);