Changeset 18665

Timestamp:

10/21/14 13:43:41 (10 years ago)

Author:

seroussi

Message:

BUG: fixed problem with gauss points in enthalpy Pvector ice sheet

File:

• issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp (modified) (1 diff)

issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp

@@ -595 +595 @@
 
         /* 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);