Index: /issm/trunk/src/c/objects/Elements/Penta.cpp
===================================================================
--- /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 4653)
+++ /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 4654)
@@ -592,5 +592,5 @@
 	double surface[numgrids];
 	double xyz_list[numgrids][3];
-	double gauss[numgrids][numgrids]={{1,0,0,0},{0,1,0,0},{0,0,1,0},{1,0,0,1},{0,1,0,1},{0,0,1,1}};
+	double gauss[numgrids][4]={{1,0,0,-1},{0,1,0,-1},{0,0,1,-1},{1,0,0,1},{0,1,0,1},{0,0,1,1}};
 
 	/*inputs: */
@@ -2095,5 +2095,5 @@
 	double       rho_ice,g;
 	double       xyz_list[numvertices][3];
-	double       gauss[numvertices][numvertices]={{1,0,0,0},{0,1,0,0},{0,0,1,0},{1,0,0,1},{0,1,0,1},{0,0,1,1}};
+	double       gauss[numvertices][4]={{1,0,0,-1},{0,1,0,-1},{0,0,1,-1},{1,0,0,1},{0,1,0,1},{0,0,1,1}};
 	
 	Input*       VzInput=NULL;
@@ -2176,5 +2176,5 @@
 	double       rho_ice,g;
 	double       xyz_list[numvertices][3];
-	double       gauss[numvertices][numvertices]={{1,0,0,0},{0,1,0,0},{0,0,1,0},{1,0,0,1},{0,1,0,1},{0,0,1,1}};
+	double       gauss[numvertices][4]={{1,0,0,-1},{0,1,0,-1},{0,0,1,-1},{1,0,0,1},{0,1,0,1},{0,0,1,1}};
 	
 	Input*       VzInput=NULL;
@@ -2257,5 +2257,5 @@
 	double       rho_ice,g;
 	double       xyz_list[numvertices][3];
-	double       gauss[numvertices][numvertices]={{1,0,0,0},{0,1,0,0},{0,0,1,0},{1,0,0,1},{0,1,0,1},{0,0,1,1}};
+	double       gauss[numvertices][4]={{1,0,0,-1},{0,1,0,-1},{0,0,1,-1},{1,0,0,1},{0,1,0,1},{0,0,1,1}};
 
 	Input*       VxInput=NULL;
@@ -2814,4 +2814,10 @@
 		values[i*numdofpervertex+2]=vz;
 		values[i*numdofpervertex+3]=p/stokesreconditioning;
+		if (
+					id==1089 ||
+					id==2761 ||
+					id==1121 ||
+					id==2762 
+			) printf("element %i vertex %i p=%g  pstokes=%g\n",id,i,p,p/stokesreconditioning);
 	}
 
@@ -5869,5 +5875,5 @@
 	int       grid=-1;
 	int       i;
-	double gauss[numnodes][numnodes]={{1,0,0,0},{0,1,0,0},{0,0,1,0},{1,0,0,1},{0,1,0,1},{0,0,1,1}};
+	double gauss[numnodes][4]={{1,0,0,-1},{0,1,0,-1},{0,0,1,-1},{1,0,0,1},{0,1,0,1},{0,0,1,1}};
 
 	/*go through 3 nodes (all nodes for tria) and identify 1st and 2nd nodes: */