Index: /issm/trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp	(revision 17370)
+++ /issm/trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp	(revision 17371)
@@ -98,7 +98,6 @@
 	IssmDouble Jdet, dt, D_scalar;
 	IssmDouble h,hx,hy,hz;
-	IssmDouble vel;
-	IssmDouble calvingrate;
-	IssmDouble dlsf[dim], norm_dlsf, normal[dim];
+	IssmDouble vel, calvingrate;
+	IssmDouble norm_dlsf;
 	IssmDouble* xyz_list = NULL;
 
@@ -115,4 +114,5 @@
 	IssmDouble*    w        = xNew<IssmDouble>(dim);
 	IssmDouble*    c        = xNew<IssmDouble>(dim);
+	IssmDouble*    dlsf     = xNew<IssmDouble>(dim);
 
 	/*Retrieve all inputs and parameters*/
@@ -156,8 +156,6 @@
 		norm_dlsf=sqrt(norm_dlsf);
 
-		if(norm_dlsf>1.e-10){
-			for(i=0;i<dim;i++) normal[i]=dlsf[i]/norm_dlsf;
-			for(i=0;i<dim;i++) c[i]=calvingrate*normal[i];
-		}
+		if(norm_dlsf>1.e-10)
+			for(i=0;i<dim;i++) c[i]=calvingrate*dlsf[i]/norm_dlsf;
 		else
 			for(i=0;i<dim;i++) c[i]=0.;
@@ -230,4 +228,5 @@
 	xDelete<IssmDouble>(w);
 	xDelete<IssmDouble>(c);
+	xDelete<IssmDouble>(dlsf);
 	delete gauss;
 	return Ke;