Index: /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp	(revision 16927)
+++ /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp	(revision 16928)
@@ -3364,5 +3364,5 @@
 
 		element->NormalBase(&bed_normal[0],xyz_list_tria);
-		element->ViscosityFS(&viscosity,xyz_list,gauss,vx_input,vy_input,vz_input);
+		element->ViscosityFS(&viscosity,3,xyz_list,gauss,vx_input,vy_input,vz_input);
 		friction->GetAlpha2(&alpha2_gauss, gauss,vx_input,vy_input,vz_input);
 
@@ -3435,5 +3435,5 @@
 		element->NodalFunctionsP1Derivatives(&dbasis[0][0],xyz_list,gauss);
 		
-		element->ViscosityFS(&viscosity,xyz_list,gauss,vx_input,vy_input,vz_input);
+		element->ViscosityFS(&viscosity,3,xyz_list,gauss,vx_input,vy_input,vz_input);
 		vzHO_input->GetInputDerivativeValue(&dw[0],xyz_list,gauss);
 
@@ -3456,5 +3456,4 @@
 	return pe;
 }/*}}}*/
-
 void StressbalanceAnalysis::GetBSSAHO(IssmDouble* B,Element* element,IssmDouble* xyz_list,Gauss* gauss){/*{{{*/
 	/*Compute B  matrix. B=[B1 B2 B3 B4 B5 B6] where Bi is of size 3*NDOF2.