Index: /issm/trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp	(revision 26433)
+++ /issm/trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp	(revision 26434)
@@ -321,12 +321,11 @@
 			case 5:{
 				/*SUPG*/
-				mf_vx_input->GetInputAverage(&w[0]);
-				mf_vy_input->GetInputAverage(&w[1]);
-				vel=sqrt(w[0]*w[0]+w[1]*w[1])+1.e-8;
+				IssmDouble vx,vy;
+				mf_vx_input->GetInputAverage(&vx);
+				mf_vy_input->GetInputAverage(&vy);
+				vel=sqrt(vx*vx+vy*vy)+1.e-8;
 				IssmPDouble xi=1.;
 				IssmDouble  tau=xi*h/(2*vel);
 
-				IssmDouble vx,vy;
-				vx = w[0]; vy = w[1];
 				/*Mass matrix - part 2*/
 				for(int i=0;i<numnodes;i++){
@@ -371,5 +370,5 @@
 	/*Fetch number of nodes and dof for this finite element*/
 	int numnodes = basalelement->GetNumberOfNodes();
-	basalelement->FindParam(&stabilization,MasstransportStabilizationEnum);
+	basalelement->FindParam(&stabilization,LevelsetStabilizationEnum);
 
 	/*Initialize Element vector*/
@@ -406,5 +405,5 @@
 			mf_vy_input->GetInputAverage(&vy);
 			vel=sqrt(vx*vx+vy*vy)+1.e-8;
-			IssmPDouble xi=1;
+			IssmPDouble xi=1.;
 			IssmDouble  tau=xi*h/(2*vel);