Index: /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp	(revision 17243)
+++ /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp	(revision 17244)
@@ -2056,5 +2056,5 @@
 	/*Fetch number of nodes and dof for this finite element*/
 	int numnodes = element->GetNumberOfNodes();
-	int numdof   = numnodes*2;
+	int numdof   = numnodes*(dim-1);
 
 	/*Initialize Element matrix and vectors*/
@@ -2099,5 +2099,5 @@
 
 	/*Transform Coordinate System*/
-	element->TransformStiffnessMatrixCoord(Ke,XYEnum);
+	if(dim==3) element->TransformStiffnessMatrixCoord(Ke,XYEnum);
 
 	/*Clean up and return*/
@@ -2185,5 +2185,5 @@
 
 	/*Transform Coordinate System*/
-	element->TransformStiffnessMatrixCoord(Ke,XYEnum);
+	if(dim==3) element->TransformStiffnessMatrixCoord(Ke,XYEnum);
 
 	/*Clean up and return*/
@@ -2250,5 +2250,5 @@
 
 	/*Transform coordinate system*/
-	element->TransformLoadVectorCoord(pe,XYEnum);
+	if(dim==3) element->TransformLoadVectorCoord(pe,XYEnum);
 
 	/*Clean up and return*/
@@ -2308,5 +2308,6 @@
 		gauss->GaussPoint(ig);
 		surface_input->GetInputValue(&surface,gauss);
-		z=element->GetZcoord(gauss);
+		if(dim==3) z=element->GetZcoord(gauss);
+		else       z=element->GetYcoord(gauss);
 		element->NodalFunctions(basis,gauss);
 		element->JacobianDeterminantSurface(&Jdet,xyz_list_front,gauss);
@@ -2323,5 +2324,5 @@
 
 	/*Transform coordinate system*/
-	element->TransformLoadVectorCoord(pe,XYEnum);
+	if(dim==3)element->TransformLoadVectorCoord(pe,XYEnum);
 
 	/*Clean up and return*/
@@ -2506,5 +2507,5 @@
 
 	/*Transform solution in Cartesian Space*/
-	element->TransformSolutionCoord(&values[0],XYEnum);
+	if(dim==3) element->TransformSolutionCoord(&values[0],XYEnum);
 
 	/*Ok, we have vx and vy in values, fill in vx and vy arrays: */
Index: /issm/trunk-jpl/src/c/classes/Materials/Matice.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Materials/Matice.cpp	(revision 17243)
+++ /issm/trunk-jpl/src/c/classes/Materials/Matice.cpp	(revision 17244)
@@ -406,5 +406,5 @@
 	}
 	else{
-		if((epsilon[0]==0) && (epsilon[1]==0) && (epsilon[2]==0)){
+		if((epsilon[0]==0) && (epsilon[1]==0)){
 			viscosity=0.5*pow(10.,14);
 		}