Index: /issm/trunk/src/c/objects/Elements/Penta.cpp
===================================================================
--- /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 3975)
+++ /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 3976)
@@ -2116,5 +2116,4 @@
 		}
 	}
-
 
 	/*Add Ke_gg to global matrix Kgg: */
Index: /issm/trunk/src/c/objects/Elements/Tria.cpp
===================================================================
--- /issm/trunk/src/c/objects/Elements/Tria.cpp	(revision 3975)
+++ /issm/trunk/src/c/objects/Elements/Tria.cpp	(revision 3976)
@@ -599,5 +599,27 @@
 /*FUNCTION Tria::UpdateInputsFromSolutionBalancedthickness {{{1*/
 void  Tria::UpdateInputsFromSolutionBalancedthickness(double* solution, int analysis_type, int sub_analysis_type){
-	ISSMERROR(" not supported yet!");
+
+	int i;
+
+	const int    numvertices=3;
+	const int    numdofpervertex=1;
+	const int    numdof=numdofpervertex*numvertices;
+
+	int          doflist[numdof];
+	double       values[numdof];
+	double       thickness[numvertices];
+
+	int          dummy;
+
+	/*Get dof list: */
+	GetDofList(&doflist[0],&dummy);
+
+	/*Use the dof list to index into the solution vector: */
+	for(i=0;i<numdof;i++){
+		values[i]=solution[doflist[i]];
+	}
+
+	/*Add thickness as inputs to the tria element: */
+	this->inputs->AddInput(new TriaVertexInput(ThicknessEnum,values));
 }
 /*}}}*/