Index: /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 15751)
+++ /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 15752)
@@ -406,5 +406,5 @@
 			int approximation;
 			inputs->GetInputValue(&approximation,ApproximationEnum);
-			if(approximation==HOFSApproximationEnum){
+			if(approximation==HOFSApproximationEnum || approximation==SSAFSApproximationEnum){
 				//Do nothing condensatino already done for Stokes part
 			}
@@ -581,5 +581,5 @@
 				int approximation;
 				inputs->GetInputValue(&approximation,ApproximationEnum);
-				if(approximation==HOFSApproximationEnum){
+				if(approximation==HOFSApproximationEnum || approximation==SSAFSApproximationEnum){
 					//Do nothing, condensation already done in PVectorCoupling
 				}
@@ -7406,6 +7406,11 @@
 
 	/*compute all stiffness matrices for this element*/
-	ElementMatrix* Ke1=CreateKMatrixDiagnosticSSA3d();
-	ElementMatrix* Ke2=CreateKMatrixDiagnosticFS();
+	ElementMatrix* Ke1=CreateKMatrixDiagnosticFS();
+	int indices[3]={18,19,20};
+	Ke1->StaticCondensation(3,&indices[0]);
+	int init = this->element_type;
+	this->element_type=P1Enum; //P1 needed for HO 
+	ElementMatrix* Ke2=CreateKMatrixDiagnosticSSA3d();
+	this->element_type=init;
 	ElementMatrix* Ke3=CreateKMatrixCouplingSSAFS();
 	ElementMatrix* Ke =new ElementMatrix(Ke1,Ke2,Ke3);
@@ -8427,6 +8432,15 @@
 
 	/*compute all load vectors for this element*/
+	int init = this->element_type;
+	this->element_type=P1Enum;
 	ElementVector* pe1=CreatePVectorDiagnosticSSA();
+	this->element_type=init;
 	ElementVector* pe2=CreatePVectorDiagnosticFS();
+	int indices[3]={18,19,20};
+	this->element_type=MINIcondensedEnum;
+	ElementMatrix* Ke = CreateKMatrixDiagnosticFS();
+	this->element_type=init;
+	pe2->StaticCondensation(Ke,3,&indices[0]);
+	delete Ke;
 	ElementVector* pe3=CreatePVectorCouplingSSAFS();
 	ElementVector* pe =new ElementVector(pe1,pe2,pe3);