Index: /issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp	(revision 18664)
+++ /issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp	(revision 18665)
@@ -595,8 +595,9 @@
 
 	/* Start  looping on the number of gaussian points: */
-	GaussPenta* gauss=new GaussPenta(2,2);
-	GaussPenta* gaussup=new GaussPenta(2,2);
+	Gauss* gauss=element->NewGaussBase(2);
+	Gauss* gaussup=element->NewGaussTop(2);
 	for(int ig=gauss->begin();ig<gauss->end();ig++){
 		gauss->GaussPoint(ig);
+		gaussup->GaussPoint(ig);
 
 		element->JacobianDeterminantBase(&Jdet,xyz_list_base,gauss);