Index: ../trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp
===================================================================
--- ../trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp	(revision 18664)
+++ ../trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp	(revision 18665)
@@ -594,10 +594,11 @@
 	Friction* friction=new Friction(element,3);
 
 	/* 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);
 		element->NodalFunctions(basis,gauss);
