Index: /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp	(revision 17620)
+++ /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp	(revision 17621)
@@ -3555,4 +3555,6 @@
 	xDelete<IssmDouble>(D);
 	xDelete<IssmDouble>(tau);
+	xDelete<IssmDouble>(vdbasis);
+	xDelete<IssmDouble>(tbasis);
 	return pe;
 }/*}}}*/
@@ -4130,11 +4132,12 @@
 
 		/*Clean up*/
+		delete gauss;
 		xDelete<IssmDouble>(xyz_list);
-		xDelete<IssmDouble>(epsxx);
-		xDelete<IssmDouble>(epsyy);
-		xDelete<IssmDouble>(epszz);
-		xDelete<IssmDouble>(epsxy);
-		xDelete<IssmDouble>(epsxz);
-		xDelete<IssmDouble>(epsyz);
+		xDelete<IssmDouble>(epsxx); xDelete<IssmDouble>(sigmapxx);
+		xDelete<IssmDouble>(epsyy); xDelete<IssmDouble>(sigmapyy);
+		xDelete<IssmDouble>(epszz); xDelete<IssmDouble>(sigmapzz);
+		xDelete<IssmDouble>(epsxy); xDelete<IssmDouble>(sigmapxy);
+		xDelete<IssmDouble>(epsxz); xDelete<IssmDouble>(sigmapxz);
+		xDelete<IssmDouble>(epsyz); xDelete<IssmDouble>(sigmapyz);
 	}
 
@@ -4418,5 +4421,6 @@
 			epsyz_input=element->GetInput(StrainRateyzEnum); _assert_(epsyz_input);
 		}
-		Gauss* gauss = element->NewGauss();
+		delete gauss;
+		gauss = element->NewGauss();
 		for(int ig=0;ig<tnumnodes;ig++){
 			gauss->GaussNode(P1DGEnum,ig);
@@ -4469,4 +4473,5 @@
 
 		/*Clean up and */
+		delete gauss;
 		xDelete<IssmDouble>(xyz_list);
 		xDelete<IssmDouble>(tbasis);
@@ -4479,7 +4484,4 @@
 		xDelete<IssmDouble>(pe_yz); xDelete<IssmDouble>(d_yz);   xDelete<IssmDouble>(tau_yz);
 	}
-
-	delete gauss;
-
 }/*}}}*/