Index: /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp	(revision 26133)
+++ /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp	(revision 26134)
@@ -2727,8 +2727,20 @@
 	if(!element->IsOnBase() || element->IsFloating()) return NULL;
 	Element* basalelement = element->SpawnBasalElement();
-	ElementMatrix* Ke = CreateKMatrixSSAFriction(basalelement);
+	ElementMatrix* Ke    = basalelement->NewElementMatrix(MLHOApproximationEnum);
+	ElementMatrix* KeSSA = CreateKMatrixSSAFriction(basalelement); //only to get K11 and K33
+	
+	/*Fetch number of nodes and dof for this finite element*/
+	int numnodes = basalelement->GetNumberOfNodes();
+
+	for(int i=0;i<numnodes;i++){
+      for(int j=0;j<numnodes;j++){
+         Ke->values[(i+0)*2*2*numnodes+j+0] = KeSSA->values[2*i*2*numnodes+2*j]; //K11
+         Ke->values[(i+6)*2*2*numnodes+j+6] = KeSSA->values[(2*i+1)*2*numnodes+2*j+1]; //K33
+      }  
+   }
 
 	/*clean-up and return*/
 	if(basalelement->IsSpawnedElement()){basalelement->DeleteMaterials(); delete basalelement;};
+	delete KeSSA;
 	return Ke;
 }/*}}}*/