Index: /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp	(revision 16928)
+++ /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp	(revision 16929)
@@ -3308,4 +3308,5 @@
 	/*Intermediaries*/
 	int         i,j,approximation;
+	int         dim=3;
 	IssmDouble  Jdet,Jdet2d,FSreconditioning;
 	IssmDouble	bed_normal[3];
@@ -3364,5 +3365,5 @@
 
 		element->NormalBase(&bed_normal[0],xyz_list_tria);
-		element->ViscosityFS(&viscosity,3,xyz_list,gauss,vx_input,vy_input,vz_input);
+		element->ViscosityFS(&viscosity,dim,xyz_list,gauss,vx_input,vy_input,vz_input);
 		friction->GetAlpha2(&alpha2_gauss, gauss,vx_input,vy_input,vz_input);
 
@@ -3392,4 +3393,5 @@
 	/*Intermediaries */
 	int         i,approximation;
+	int         dim=3;
 	IssmDouble  viscosity,Jdet,FSreconditioning;
 	IssmDouble  dw[3];
@@ -3435,5 +3437,5 @@
 		element->NodalFunctionsP1Derivatives(&dbasis[0][0],xyz_list,gauss);
 		
-		element->ViscosityFS(&viscosity,3,xyz_list,gauss,vx_input,vy_input,vz_input);
+		element->ViscosityFS(&viscosity,dim,xyz_list,gauss,vx_input,vy_input,vz_input);
 		vzHO_input->GetInputDerivativeValue(&dw[0],xyz_list,gauss);