Index: /issm/trunk-jpl/src/c/modules/SurfaceMassBalancex/Gembx.cpp
===================================================================
--- /issm/trunk-jpl/src/c/modules/SurfaceMassBalancex/Gembx.cpp	(revision 21227)
+++ /issm/trunk-jpl/src/c/modules/SurfaceMassBalancex/Gembx.cpp	(revision 21228)
@@ -1239,5 +1239,4 @@
     
 	int        D_size;
-	int         i;
 
 	/*outputs:*/
@@ -1566,5 +1565,5 @@
         IssmDouble m_new = m[X2] + m[X1];
         T[X1] = (T[X2]*m[X2] + T[X1]*m[X1]) / m_new;
-        a[X1] = (a[X2]*m[i] + a[X1]*m[X1]) / m_new;
+        a[X1] = (a[X2]*m[X2] + a[X1]*m[X1]) / m_new;
         re[X1] = (re[X2]*m[X2] + re[X1]*m[X1]) / m_new;
         gdn[X1] = (gdn[X2]*m[X2] + gdn[X1]*m[X1]) / m_new;
@@ -1608,5 +1607,5 @@
 	}
 	
-	i=0;
+	int i=0;
 	while(i<=X){
 		if (dz [i] > dzMin *2){
@@ -1727,5 +1726,7 @@
 	if(D)xDelete<int>(D);
 	if(M)xDelete<IssmDouble>(M);
-	
+    xDelete<IssmDouble>(Zcum);
+    xDelete<IssmDouble>(dzMin2);
+    
 	/*Assign output pointers:*/
 	*pM=sumM;