Index: /issm/trunk-jpl/src/c/classes/Elements/Element.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Element.cpp	(revision 27448)
+++ /issm/trunk-jpl/src/c/classes/Elements/Element.cpp	(revision 27449)
@@ -3922,12 +3922,12 @@
 
    /* Start looping on the number of vertices: */
-   GaussTria gauss;
+	Gauss* gauss=this->NewGauss();
    for(int iv=0;iv<numvertices;iv++){
-      gauss.GaussVertex(iv);
+		gauss->GaussVertex(iv);
 
       /* Get variables */
-      bed_input->GetInputValue(&bed,&gauss);
-      qsg_input->GetInputValue(&qsg,&gauss);
-		TF_input->GetInputValue(&TF,&gauss);
+      bed_input->GetInputValue(&bed,gauss);
+      qsg_input->GetInputValue(&qsg,gauss);
+		TF_input->GetInputValue(&TF,gauss);
 
       if(basin_icefront_area[basinid]==0.) meltrates[iv]=0.;