Index: /issm/trunk/src/c/objects/Elements/Penta.cpp
===================================================================
--- /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 5588)
+++ /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 5589)
@@ -2917,5 +2917,4 @@
 
 	/*matrices: */
-	double     K_terms[numdof][numdof]={0.0};
 	double     Ke_temp[27][27]={0.0}; //for the six nodes and the bubble 
 	double     Ke_reduced[numdof][numdof]; //for the six nodes only
@@ -2974,4 +2973,5 @@
 	/*retrive parameters: */
 	parameters->FindParam(&analysis_type,AnalysisTypeEnum);
+	this->parameters->FindParam(&stokesreconditioning,StokesReconditioningEnum);
 
 	/*retrieve inputs :*/
@@ -2983,7 +2983,4 @@
 	/*If on water or not Stokes, skip stiffness: */
 	if(onwater || approximation!=StokesApproximationEnum) return;
-
-	/*retrieve some parameters: */
-	this->parameters->FindParam(&stokesreconditioning,StokesReconditioningEnum);
 
 	/* Get node coordinates and dof list: */
@@ -3139,12 +3136,6 @@
 	ReduceMatrixStokes(&Ke_reduced[0][0], &Ke_temp[0][0]);
 
-	for(i=0;i<numdof;i++){
-		for(j=0;j<numdof;j++){
-			K_terms[i][j]+=Ke_reduced[i][j];
-		}
-	}
-
 	/*Add Ke_gg to global matrix Kgg: */
-	MatSetValues(Kgg,numdof,doflist,numdof,doflist,(const double*)K_terms,ADD_VALUES);
+	MatSetValues(Kgg,numdof,doflist,numdof,doflist,(const double*)Ke_reduced,ADD_VALUES);
 
 	/*Free ressources:*/