Index: /issm/trunk-jpl/src/c/classes/objects/Elements/Penta.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/objects/Elements/Penta.cpp	(revision 13041)
+++ /issm/trunk-jpl/src/c/classes/objects/Elements/Penta.cpp	(revision 13042)
@@ -6327,18 +6327,18 @@
 	/*Intermediaries */
 	int         i,j,ig,approximation;
-	IssmDouble      Jdet;
-	IssmDouble      viscosity, oldviscosity, newviscosity, viscosity_overshoot;
-	IssmDouble      epsilon[5],oldepsilon[5]; /* epsilon=[exx,eyy,exy,exz,eyz];*/
-	IssmDouble      epsilons[6]; //6 for stokes
-	IssmDouble      xyz_list[NUMVERTICES][3];
-	IssmDouble      B[3][numdof2d];
-	IssmDouble      Bprime[3][numdof2d];
-	IssmDouble      D[3][3]={0.0};            // material matrix, simple scalar matrix.
-	IssmDouble      D_scalar;
-	IssmDouble      Ke_gg_gaussian[numdof2d][numdof2d]; //stiffness matrix evaluated at the gaussian point.
-	Tria*       tria=NULL;
-	Penta*      pentabase=NULL;
-	GaussPenta *gauss=NULL;
-	GaussTria  *gauss_tria=NULL;
+	IssmDouble  Jdet;
+	IssmDouble  viscosity , oldviscosity, newviscosity, viscosity_overshoot;
+	IssmDouble  epsilon[5],oldepsilon[5];       /* epsilon=[exx,eyy,exy,exz,eyz];*/
+	IssmDouble  epsilons[6];                    //6 for stokes
+	IssmDouble  xyz_list[NUMVERTICES][3];
+	IssmDouble  B[3][numdof2d];
+	IssmDouble  Bprime[3][numdof2d];
+	IssmDouble  D[3][3]= {0.0};                 // material matrix, simple scalar matrix.
+	IssmDouble  D_scalar;
+	IssmDouble  Ke_gg_gaussian[numdof2d][numdof2d];
+	Tria       *tria       = NULL;
+	Penta      *pentabase  = NULL;
+	GaussPenta *gauss      = NULL;
+	GaussTria  *gauss_tria = NULL;
 
 	/*Find penta on bed as this is a macayeal elements: */