Index: /issm/trunk-jpl/src/c/classes/Inputs/PentaInput.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Inputs/PentaInput.cpp	(revision 17228)
+++ /issm/trunk-jpl/src/c/classes/Inputs/PentaInput.cpp	(revision 17229)
@@ -94,5 +94,5 @@
 		IssmDouble newvalues[3]; 
 
-		/*Create arrow of indices depending on location (0=base 1=surface)*/
+		/*Create array of indices depending on location (0=base 1=surface)*/
 		int indices[3];
 		switch(location){
Index: /issm/trunk-jpl/src/c/classes/Inputs/TriaInput.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Inputs/TriaInput.cpp	(revision 17228)
+++ /issm/trunk-jpl/src/c/classes/Inputs/TriaInput.cpp	(revision 17229)
@@ -100,13 +100,19 @@
 	/*output*/
 	SegInput* outinput=NULL;
-	IssmDouble newvalues[2]; //Assume P1 interpolation only for now
-
-	/*Create arrow of indices depending on location (0=base 1=surface)*/
-
-	newvalues[0]=this->values[index1];
-	newvalues[1]=this->values[index2];
-
-	/*Create new Seg input*/
-	outinput=new SegInput(this->enum_type,&newvalues[0],P1Enum);
+
+	if(this->element_type==P0Enum){
+		outinput=new SegInput(this->enum_type,&this->values[0],P0Enum);
+	}
+	else{
+		/*Assume P1 interpolation only for now*/
+		IssmDouble newvalues[2];
+
+		/*Create array of indices depending on location (0=base 1=surface)*/
+		newvalues[0]=this->values[index1];
+		newvalues[1]=this->values[index2];
+
+		/*Create new Seg input*/
+		outinput=new SegInput(this->enum_type,&newvalues[0],P1Enum);
+	}
 
 	/*Assign output*/