Index: /issm/trunk/src/c/objects/Elements/Tria.cpp
===================================================================
--- /issm/trunk/src/c/objects/Elements/Tria.cpp	(revision 4493)
+++ /issm/trunk/src/c/objects/Elements/Tria.cpp	(revision 4494)
@@ -5871,5 +5871,25 @@
 /*FUNCTION Tria::InputUpdateFromSolutionBalancedthickness2 {{{1*/
 void  Tria::InputUpdateFromSolutionBalancedthickness2(double* solution){
-	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];
+
+	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));
 }
 /*}}}*/
Index: /issm/trunk/src/c/solutions/balancedthickness2_core.cpp
===================================================================
--- /issm/trunk/src/c/solutions/balancedthickness2_core.cpp	(revision 4493)
+++ /issm/trunk/src/c/solutions/balancedthickness2_core.cpp	(revision 4494)
@@ -20,5 +20,5 @@
 
 	/*activate formulation: */
-	femmodel->SetCurrentConfiguration(BalancedthicknessAnalysisEnum);
+	femmodel->SetCurrentConfiguration(Balancedthickness2AnalysisEnum);
 	
 	/*recover parameters: */
@@ -34,10 +34,9 @@
 	solver_linear(NULL,femmodel);
 
-	_printf_("Averaging over vertices:\n");
-	ISSMERROR(" not supported yet!");
-//FieldAverageOntoVerticesx(&h_g,fem_p->elements,fem_p->nodes,fem_p->vertices,fem_p->loads,fem_p->materials,fem_p->parameters);
+	//_printf_("Averaging over vertices:\n");
+	//FieldAverageOntoVerticesx(&h_g,fem_p->elements,fem_p->nodes,fem_p->vertices,fem_p->loads,fem_p->materials,fem_p->parameters);
 
-//	_printf_("extrude computed thickness on all layers:\n");
-//	FieldExtrudex( h_g, fem_p->elements,fem_p->nodes, fem_p->vertices,fem_p->loads, fem_p->materials,fem_p->parameters,"thickness",0);
+	//	_printf_("extrude computed thickness on all layers:\n");
+	//	FieldExtrudex( h_g, fem_p->elements,fem_p->nodes, fem_p->vertices,fem_p->loads, fem_p->materials,fem_p->parameters,"thickness",0);
 
 	if(verbose)_printf_("saving results:\n");
Index: /issm/trunk/src/m/solutions/balancedthickness2_core.m
===================================================================
--- /issm/trunk/src/m/solutions/balancedthickness2_core.m	(revision 4493)
+++ /issm/trunk/src/m/solutions/balancedthickness2_core.m	(revision 4494)
@@ -14,9 +14,6 @@
 
 	displaystring(verbose,'\n%s',['depth averaging velocities...']);
-	%[femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters]=InputDepthAverage(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VxEnum,VxAverageEnum);
-	%[femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters]=InputDepthAverage(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VyEnum,VyAverageEnum);
-	if dim==3,
-	%	[femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters]=InputDepthAverage(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VzEnum,VzAverageEnum);
-	end
+	[femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters]=InputDepthAverage(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VxEnum,VxAverageEnum);
+	[femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters]=InputDepthAverage(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VyEnum,VyAverageEnum);
 
 	displaystring(verbose,'\n%s',['call computational core...']);