Index: /issm/trunk/src/c/objects/Elements/Penta.cpp
===================================================================
--- /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 5977)
+++ /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 5978)
@@ -2397,6 +2397,6 @@
 		matice->GetViscosity3d(&viscosity, &epsilon[0]);
 		matice->GetViscosity3d(&oldviscosity, &oldepsilon[0]);
-
 		newviscosity=viscosity+viscosity_overshoot*(viscosity-oldviscosity);
+
 		D_scalar=2*newviscosity*gauss->weight*Jdet;
 		for (i=0;i<5;i++) D[i][i]=D_scalar;
@@ -3438,5 +3438,4 @@
 	double      l1l6[6]; //for the six nodes of the penta
 	double      Jdet2d;
-	Tria*       tria=NULL;
 	GaussPenta  *gauss=NULL;
 
@@ -4658,7 +4657,6 @@
 	int          dummy;
 	double       pressure[NUMVERTICES];
-	double       surface[NUMVERTICES];
-	double       rho_ice,g;
 	double       xyz_list[NUMVERTICES][3];
+	double       stokesreconditioning;
 
 	double *vzpattyn_ptr         = NULL;
@@ -4672,4 +4670,5 @@
 	GetDofList(&doflistp,PattynApproximationEnum,GsetEnum);
 	GetDofList(&doflists,StokesApproximationEnum,GsetEnum);
+	this->parameters->FindParam(&stokesreconditioning,StokesReconditioningEnum);
 
 	/*Get node data: */
@@ -4689,5 +4688,5 @@
 		vy[i]=stokes_values[i*numdofpervertexs+1]+pattyn_values[i*numdofpervertexp+1];
 		vzstokes[i]=stokes_values[i*numdofpervertexs+2];
-		pressure[i]=stokes_values[i*numdofpervertexs+3];
+		pressure[i]=stokes_values[i*numdofpervertexs+3]*stokesreconditioning;
 	}