Index: /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp	(revision 16989)
+++ /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp	(revision 16990)
@@ -3199,8 +3199,20 @@
 
 	/*compute all stiffness matrices for this element*/
-	ElementMatrix* Ke1=CreateKMatrixSSA(element);
+	ElementMatrix* Ke1=CreateKMatrixSSA3d(element);
 	ElementMatrix* Ke2=CreateKMatrixHO(element);
 	ElementMatrix* Ke3=CreateKMatrixCouplingSSAHO(element);
 	ElementMatrix* Ke =new ElementMatrix(Ke1,Ke2,Ke3);
+
+	//Ke2=element->CreateKMatrix();
+	//for(int i=0;i<12;i++){
+	//	for(int j=0;j<12;j++){
+	//		if(Ke->values[i*12+j] != Ke2->values[i*12+j]){
+	//			Ke->Echo();
+	//			Ke2->Echo();
+	//			printf("ID = %i\n",element->Id());
+	//			_error_("S");
+	//		}
+	//	}
+	//}
 
 	/*clean-up and return*/