Index: /issm/trunk-jpl/src/c/classes/Elements/Tetra.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Tetra.cpp	(revision 17749)
+++ /issm/trunk-jpl/src/c/classes/Elements/Tetra.cpp	(revision 17750)
@@ -107,6 +107,5 @@
 
 	IssmDouble values[NUMVERTICES];
-	//int        indices[4][3] = {{0,1,2},{0,1,3},{1,2,3},{2,0,3}};
-	int        indices[4][3] = {{0,1,2},{3,1,0},{3,2,1},{3,0,2}};
+	int        indices[4][3] = {{0,1,2},{0,3,1},{1,3,2},{0,2,3}};
 
 	/*Retrieve all inputs and parameters*/
Index: /issm/trunk-jpl/src/c/classes/gauss/GaussTetra.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/gauss/GaussTetra.cpp	(revision 17749)
+++ /issm/trunk-jpl/src/c/classes/gauss/GaussTetra.cpp	(revision 17750)
@@ -53,15 +53,15 @@
 		for(int i=0;i<numgauss;i++) coords4[i]=0.;
 	}
-	else if(index1==0 && index2==1 && index3==3){
+	else if(index1==0 && index2==3 && index3==1){
 		GaussLegendreTria(&numgauss,&coords1,&coords2,&coords4,&weights,order);
 		coords3=xNew<IssmDouble>(numgauss);
 		for(int i=0;i<numgauss;i++) coords3[i]=0.;
 	}
-	else if(index1==1 && index2==2 && index3==3){
+	else if(index1==1 && index2==3 && index3==2){
 		GaussLegendreTria(&numgauss,&coords2,&coords3,&coords4,&weights,order);
 		coords1=xNew<IssmDouble>(numgauss);
 		for(int i=0;i<numgauss;i++) coords1[i]=0.;
 	}
-	else if(index1==2 && index2==0 && index3==3){
+	else if(index1==0 && index2==2 && index3==3){
 		GaussLegendreTria(&numgauss,&coords1,&coords3,&coords4,&weights,order);
 		coords2=xNew<IssmDouble>(numgauss);