Index: ../trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp
===================================================================
--- ../trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp	(revision 18783)
+++ ../trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp	(revision 18784)
@@ -271,7 +271,7 @@
 
 		/* 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){
 			case 0:
@@ -297,10 +297,10 @@
 			case 2:
 				/* 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,
 							D,dim,dim,0,
