Index: /issm/trunk-jpl/src/c/classes/objects/Elements/Tria.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/objects/Elements/Tria.cpp	(revision 14986)
+++ /issm/trunk-jpl/src/c/classes/objects/Elements/Tria.cpp	(revision 14987)
@@ -1762,9 +1762,4 @@
 			break;
 		#endif
-		#ifdef _HAVE_GIA_
-		case GiaAnalysisEnum:
-			InputUpdateFromSolutionGia(solution);
-			break;
-		#endif
 		case BedSlopeXAnalysisEnum:
 			InputUpdateFromSolutionOneDof(solution,BedSlopeXEnum);
@@ -1889,13 +1884,22 @@
 
 	switch(type){
-	#ifdef _HAVE_GIA_
-	case GiaAnalysisEnum:
-		InputUpdateFromVectorGia(vector);
-		break;
-	#endif
-	case VertexEnum: 
+	case VertexPIdEnum: 
 		/*Get values on the 3 vertices*/
 		for (int i=0;i<3;i++){
 			values[i]=vector[this->nodes[i]->GetVertexPid()];
+		}
+		/*update input*/
+		if (name==MaterialsRheologyBbarEnum || name==MaterialsRheologyBEnum || name==MaterialsRheologyZEnum || name==MaterialsRheologyZbarEnum){
+			material->inputs->AddInput(new TriaP1Input(name,values));
+		}
+		else{
+			this->inputs->AddInput(new TriaP1Input(name,values));
+		}
+		return; 
+
+	case VertexSIdEnum: 
+		/*Get values on the 3 vertices*/
+		for (int i=0;i<3;i++){
+			values[i]=vector[this->nodes[i]->GetVertexSid()];
 		}
 		/*update input*/
@@ -2072,5 +2076,5 @@
 				name==GradientEnum ||
 				name==OldGradientEnum ||
-        name==ConvergedEnum ||
+	         name==ConvergedEnum ||
 				name==BasisIntegralEnum ||
 				name==QmuVxEnum ||
@@ -3221,44 +3225,4 @@
 
 	return;
-}
-/*}}}*/
-/*FUNCTION Tria::InputUpdateFromSolutionGia{{{*/
-void  Tria::InputUpdateFromSolutionGia(IssmDouble* solution){
-
-	int vertexsidlist[NUMVERTICES];
-	IssmDouble    values[NUMVERTICES];
-
-	/*Get sid list: */
-	this->GetVertexSidList(&vertexsidlist[0]);
-
-	/*Use the dof list to index into the solution vector: */
-	for(int i=0;i<NUMVERTICES;i++){
-		values[i]=solution[vertexsidlist[i]];
-		if(xIsNan<IssmDouble>(values[i])) _error_("NaN found in solution vector");
-	}
-
-	/*Add input to the element: */
-	this->inputs->AddInput(new TriaP1Input(GiaWEnum,values));
-
-}
-/*}}}*/
-/*FUNCTION Tria::InputUpdateFromVectorGia{{{*/
-void  Tria::InputUpdateFromVectorGia(IssmDouble* vector){
-
-	int vertexsidlist[NUMVERTICES];
-	IssmDouble    values[NUMVERTICES];
-
-	/*Get sid list: */
-	this->GetVertexSidList(&vertexsidlist[0]);
-
-	/*Use the dof list to index into the solution vector: */
-	for(int i=0;i<NUMVERTICES;i++){
-		values[i]=vector[vertexsidlist[i]];
-		if(xIsNan<IssmDouble>(values[i])) _error_("NaN found in solution vector");
-	}
-
-	/*Add input to the element: */
-	this->inputs->AddInput(new TriaP1Input(GiadWdtEnum,values));
-
 }
 /*}}}*/
Index: /issm/trunk-jpl/src/c/classes/objects/Elements/Tria.h
===================================================================
--- /issm/trunk-jpl/src/c/classes/objects/Elements/Tria.h	(revision 14986)
+++ /issm/trunk-jpl/src/c/classes/objects/Elements/Tria.h	(revision 14987)
@@ -142,6 +142,4 @@
 		#ifdef _HAVE_GIA_
 		void   GiaDeflection(Vector<IssmDouble>* wg,Vector<IssmDouble>* dwgdt,IssmDouble* x,IssmDouble* y);
-		void   InputUpdateFromSolutionGia(IssmDouble* solution); 
-		void   InputUpdateFromVectorGia(IssmDouble* vector); 
 		#endif