Index: /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 15555)
+++ /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 15556)
@@ -8216,5 +8216,5 @@
 	IssmDouble zmax=xyz_list[0][2]; for(int i=1;i<6;i++) if(xyz_list[i][2]>zmax) zmax=xyz_list[i][2];
 	IssmDouble zmin=xyz_list[0][2]; for(int i=1;i<6;i++) if(xyz_list[i][2]<zmin) zmin=xyz_list[i][2];
-	if(zmax>0 && zmin<0) gauss=new GaussPenta(area_coordinates,3,10); //refined in vertical because of the sea level discontinuity
+	if(zmax>0 && zmin<0) gauss=new GaussPenta(area_coordinates,3,30); //refined in vertical because of the sea level discontinuity
 	else                 gauss=new GaussPenta(area_coordinates,3,3);
 
@@ -8234,4 +8234,5 @@
 			pe->values[4*i+1]+= pressure*Jdet*gauss->weight*normal[1]*basis[i];
 			pe->values[4*i+2]+= pressure*Jdet*gauss->weight*normal[2]*basis[i];
+			pe->values[4*i+3]+= 0;
 		}
 	}
@@ -8239,5 +8240,4 @@
 	/*Transform coordinate system*/
 	TransformLoadVectorCoord(pe,nodes,numnodes,XYZPEnum);
-
 
 	/*Clean up and return*/