Index: /issm/trunk-jpl/src/c/classes/gauss/GaussTria.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/gauss/GaussTria.cpp	(revision 25440)
+++ /issm/trunk-jpl/src/c/classes/gauss/GaussTria.cpp	(revision 25441)
@@ -309,5 +309,4 @@
 	 *
 	 */
-	int         ig;
 	IssmDouble x,y;
 	IssmDouble xy_list[3][2];
@@ -325,15 +324,15 @@
 	this->weights=xNew<IssmDouble>(this->numgauss);
 
-	for(ig=0;ig<gauss1->numgauss;ig++){ // Add the first triangle gauss points
-		this->coords1[ig]=gauss1->coords1[ig];
-		this->coords2[ig]=gauss1->coords2[ig];
-		this->coords3[ig]=gauss1->coords3[ig];
-		this->weights[ig]=gauss1->weights[ig];
-	}
-	for(ig=0;ig<gauss2->numgauss;ig++){ // Add the second triangle gauss points
-		this->coords1[gauss1->numgauss+ig]=gauss2->coords1[ig];
-		this->coords2[gauss1->numgauss+ig]=gauss2->coords2[ig];
-		this->coords3[gauss1->numgauss+ig]=gauss2->coords3[ig];
-		this->weights[gauss1->numgauss+ig]=gauss2->weights[ig];
+	for(int ii=0;ii<gauss1->numgauss;ii++){ // Add the first triangle gauss points
+		this->coords1[ii]=gauss1->coords1[ii];
+		this->coords2[ii]=gauss1->coords2[ii];
+		this->coords3[ii]=gauss1->coords3[ii];
+		this->weights[ii]=gauss1->weights[ii];
+	}
+	for(int ii=0;ii<gauss2->numgauss;ii++){ // Add the second triangle gauss points
+		this->coords1[gauss1->numgauss+ii]=gauss2->coords1[ii];
+		this->coords2[gauss1->numgauss+ii]=gauss2->coords2[ii];
+		this->coords3[gauss1->numgauss+ii]=gauss2->coords3[ii];
+		this->weights[gauss1->numgauss+ii]=gauss2->weights[ii];
 	}
 
@@ -343,4 +342,5 @@
 
 	/*Initialize static fields as undefined*/
+	ig    = -1;
 	weight=UNDEF;
 	coord1=UNDEF;