Index: /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 17231)
+++ /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 17232)
@@ -1219,5 +1219,5 @@
 	::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
 	for(i=0;i<NUMVERTICES;i++) x_list[i]=xyz_list[i][0];
-	PentaRef::GetInputValue(&x,x_list,gauss);
+	PentaRef::GetInputValue(&x,x_list,gauss,P1Enum);
 
 	return x;
@@ -1234,5 +1234,5 @@
 	::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
 	for(i=0;i<NUMVERTICES;i++) y_list[i]=xyz_list[i][1];
-	PentaRef::GetInputValue(&y,y_list,gauss);
+	PentaRef::GetInputValue(&y,y_list,gauss,P1Enum);
 
 	return y;
@@ -1249,5 +1249,5 @@
 	::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
 	for(i=0;i<NUMVERTICES;i++) z_list[i]=xyz_list[i][2];
-	PentaRef::GetInputValue(&z,z_list,gauss);
+	PentaRef::GetInputValue(&z,z_list,gauss,P1Enum);
 
 	return z;
@@ -3102,5 +3102,5 @@
 /*FUNCTION Penta::ValueP1OnGauss{{{*/
 void Penta::ValueP1OnGauss(IssmDouble* pvalue,IssmDouble* values,Gauss* gauss){
-	PentaRef::GetInputValue(pvalue,values,gauss);
+	PentaRef::GetInputValue(pvalue,values,gauss,P1Enum);
 }
 /*}}}*/
Index: /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 17231)
+++ /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 17232)
@@ -1218,8 +1218,5 @@
 	::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
 	for(int i=0;i<NUMVERTICES;i++) y_list[i]=xyz_list[i][1];
-	int backup = this->element_type;
-	this->element_type = P1Enum;
-	TriaRef::GetInputValue(&y,y_list,gauss);
-	this->element_type = backup;
+	TriaRef::GetInputValue(&y,y_list,gauss,P1Enum);
 
 	return y;
@@ -2563,5 +2560,5 @@
 /*FUNCTION Tria::ValueP1OnGauss{{{*/
 void Tria::ValueP1OnGauss(IssmDouble* pvalue,IssmDouble* values,Gauss* gauss){
-	TriaRef::GetInputValue(pvalue,values,gauss);
+	TriaRef::GetInputValue(pvalue,values,gauss,P1Enum);
 }
 /*}}}*/
