Index: /issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp	(revision 19636)
+++ /issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp	(revision 19637)
@@ -530,5 +530,5 @@
 
 	/* Start  looping on the number of gaussian points: */
-	Gauss* gauss=element->NewGauss(2);
+	Gauss* gauss=element->NewGauss(4);
 	for(int ig=gauss->begin();ig<gauss->end();ig++){
 		gauss->GaussPoint(ig);
@@ -648,5 +648,5 @@
 
 	/* Start  looping on the number of gaussian points: */
-	Gauss* gauss=element->NewGaussBase(2);
+	Gauss* gauss=element->NewGaussBase(4);
 	for(int ig=gauss->begin();ig<gauss->end();ig++){
 		gauss->GaussPoint(ig);
@@ -737,5 +737,5 @@
 
 	/* Start  looping on the number of gaussian points: */
-	Gauss* gauss=element->NewGauss(3);
+	Gauss* gauss=element->NewGauss(4);
 	for(int ig=gauss->begin();ig<gauss->end();ig++){
 		gauss->GaussPoint(ig);
@@ -844,6 +844,6 @@
 
 	/* Start  looping on the number of gaussian points: */
-	Gauss* gauss=element->NewGaussBase(2);
-	Gauss* gaussup=element->NewGaussTop(2);
+	Gauss* gauss=element->NewGaussBase(4);
+	Gauss* gaussup=element->NewGaussTop(4);
 	for(int ig=gauss->begin();ig<gauss->end();ig++){
 		gauss->GaussPoint(ig);
@@ -938,5 +938,5 @@
 
 	/* Start  looping on the number of gaussian points: */
-	Gauss* gauss=element->NewGaussBase(2);
+	Gauss* gauss=element->NewGaussBase(4);
 	for(int ig=gauss->begin();ig<gauss->end();ig++){
 		gauss->GaussPoint(ig);
Index: /issm/trunk-jpl/src/c/analyses/ThermalAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/ThermalAnalysis.cpp	(revision 19636)
+++ /issm/trunk-jpl/src/c/analyses/ThermalAnalysis.cpp	(revision 19637)
@@ -247,5 +247,5 @@
 
 	/* Start  looping on the number of gaussian points: */
-	Gauss* gauss=element->NewGaussBase(2);
+	Gauss* gauss=element->NewGaussBase(4);
 	for(int ig=gauss->begin();ig<gauss->end();ig++){
 		gauss->GaussPoint(ig);
@@ -313,5 +313,5 @@
 
 	/* Start  looping on the number of gaussian points: */
-	Gauss* gauss=element->NewGauss(2);
+	Gauss* gauss=element->NewGauss(4);
 	for(int ig=gauss->begin();ig<gauss->end();ig++){
 		gauss->GaussPoint(ig);
@@ -452,5 +452,5 @@
 
 	/* Start  looping on the number of gaussian points: */
-	Gauss* gauss   = element->NewGaussBase(2);
+	Gauss* gauss   = element->NewGaussBase(4);
 	for(int ig=gauss->begin();ig<gauss->end();ig++){
 		gauss->GaussPoint(ig);
@@ -511,5 +511,5 @@
 
 	/* Start  looping on the number of gaussian points: */
-	Gauss* gauss=element->NewGaussBase(2);
+	Gauss* gauss=element->NewGaussBase(4);
 	for(int ig=gauss->begin();ig<gauss->end();ig++){
 		gauss->GaussPoint(ig);
@@ -571,5 +571,5 @@
 
 	/* Start  looping on the number of gaussian points: */
-	Gauss* gauss=element->NewGauss(3);
+	Gauss* gauss=element->NewGauss(4);
 	for(int ig=gauss->begin();ig<gauss->end();ig++){
 		gauss->GaussPoint(ig);
@@ -714,5 +714,5 @@
 	int        *doflist   = NULL;
 	IssmDouble *xyz_list  = NULL;
-	bool        hack      = false;
+	bool        hack      = true;
 
 	/*Fetch number of nodes and dof for this finite element*/