Index: /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 16973)
+++ /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 16974)
@@ -7735,5 +7735,5 @@
 
 	/* Start  looping on the number of gaussian points: */
-	gauss=new GaussPenta(0,1,2,2);
+	gauss=new GaussPenta(0,1,2,3);
 	for(int ig=gauss->begin();ig<gauss->end();ig++){
 
@@ -8841,5 +8841,5 @@
 
 	/* Start looping on the number of gauss 2d (nodes on the bedrock) */
-	gauss=new GaussPenta(0,1,2,2);
+	gauss=new GaussPenta(0,1,2,5);
 	for(int ig=gauss->begin();ig<gauss->end();ig++){
 
Index: /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 16973)
+++ /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 16974)
@@ -3600,5 +3600,5 @@
 
 	/* Start looping on the number of gauss 1d (nodes on the bedrock) */
-	GaussTria* gauss=new GaussTria(indices[0],indices[1],2);
+	GaussTria* gauss=new GaussTria(indices[0],indices[1],3);
 	for(int ig=gauss->begin();ig<gauss->end();ig++){
 
@@ -4012,5 +4012,5 @@
 
 	/* Start looping on the number of gauss 1d (nodes on the bedrock) */
-	GaussTria* gauss=new GaussTria(indices[0],indices[1],2);
+	GaussTria* gauss=new GaussTria(indices[0],indices[1],5);
 	for(int ig=gauss->begin();ig<gauss->end();ig++){