Index: /issm/trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp	(revision 18783)
+++ /issm/trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp	(revision 18784)
@@ -272,5 +272,5 @@
 		/* Stabilization */
 		vel=0.;
-		for(i=0;i<dim;i++) vel+=v[i]*v[i];
+		for(i=0;i<dim;i++) vel+=w[i]*w[i];
 		vel=sqrt(vel)+1.e-14;
 		switch(stabilization){
@@ -298,8 +298,8 @@
 				/* Streamline Upwinding */
 				basalelement->ElementSizes(&hx,&hy,&hz);
-				h=sqrt( pow(hx*v[0]/vel,2) + pow(hy*v[1]/vel,2) );
+				h=sqrt( pow(hx*w[0]/vel,2) + pow(hy*w[1]/vel,2) );
 				for(row=0;row<dim;row++) 
 					for(col=0;col<dim;col++) 
-						D[row*dim+col] = D_scalar*h/(2.*vel)*v[row]*v[col];
+						D[row*dim+col] = D_scalar*h/(2.*vel)*w[row]*w[col];
 
 				TripleMultiply(Bprime,dim,numnodes,1,