Index: /issm/trunk/src/c/parallel/ControlInitialization.cpp
===================================================================
--- /issm/trunk/src/c/parallel/ControlInitialization.cpp	(revision 2018)
+++ /issm/trunk/src/c/parallel/ControlInitialization.cpp	(revision 2019)
@@ -102,8 +102,4 @@
 	VecDuplicatePatch(&ug_horiz,ug); FieldExtrudex( ug_horiz,fem_dh->elements,fem_dh->nodes, fem_dh->loads,fem_dh-> materials,"velocity",1);
 
-	//pressure
-	if(debug)_printf_("%s\n"," computing pressure according to MacAyeal...");
-	ComputePressurex( &pg,fem_dh->elements, fem_dh->nodes, fem_dh->loads,  fem_dh->materials, numberofnodes);
-
 	//vertical velocity
 	if(debug)_printf_("%s\n"," computing vertical velocities...");
@@ -141,9 +137,7 @@
 	VecFree(&ug);
 	diagnostic_core_nonlinear(&ug,NULL,NULL,NULL,fem_ds,inputs,DiagnosticAnalysisEnum(),StokesAnalysisEnum());
-	VecFree(&pg);
 
 	//Add in inputs
 	inputs->Add("velocity",ug,numberofdofspernode_ds,numberofnodes);
-	VecFree(&ug);
 
 	/*Assign output*/
@@ -152,4 +146,11 @@
 	/*Free ressources:*/
 	xfree((void**)&dofset);
+	VecFree(&ug);
+	VecFree(&ug_horiz);
+	VecFree(&ug_vert);
+	VecFree(&ug_stokes);
+	VecFree(&pg);
+	VecFree(&slopex);
+	VecFree(&slopey);
 
 }