Index: /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 17669)
+++ /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 17670)
@@ -165,4 +165,5 @@
 	IssmDouble	sigma_yz[NUMVERTICES]={0,0,0};
 	GaussTria*  gauss=NULL;
+	int dim=2;
 
 	/* Get node coordinates and dof list: */
@@ -181,5 +182,5 @@
 		/*Compute strain rate viscosity and pressure: */
 		this->StrainRateSSA(&epsilon[0],&xyz_list[0][0],gauss,vx_input,vy_input);
-		this->ViscositySSA(&viscosity,2,&xyz_list[0][0],gauss,vx_input,vy_input);
+		this->ViscositySSA(&viscosity,dim,&xyz_list[0][0],gauss,vx_input,vy_input);
 		pressure_input->GetInputValue(&pressure,gauss);