Index: /issm/trunk/src/m/solutions/jpl/balancedthickness_core.m
===================================================================
--- /issm/trunk/src/m/solutions/jpl/balancedthickness_core.m	(revision 3974)
+++ /issm/trunk/src/m/solutions/jpl/balancedthickness_core.m	(revision 3975)
@@ -6,13 +6,22 @@
 
 	%get FE model
-	m=models.bt;
+	verbose=models.bt.parameters.Verbose;
 	results.time=0;
 	results.step=1;
 
-	displaystring(m.parameters.Verbose,'\n%s',['call computational core:']);
-	results.h_g=diagnostic_core_linear(m,analysis_type,sub_analysis_type);
+	displaystring(verbose,'\n%s',['depth averaging velocity...']);
+	[models.bt.elements,models.bt.nodes,models.bt.vertices,models.bt.loads,models.bt.materials,models.bt.parameters]=DepthAverageInput(models.bt.elements,models.bt.nodes,models.bt.vertices,models.bt.loads,models.bt.materials,models.bt.parameters,VxEnum);
+	[models.bt.elements,models.bt.nodes,models.bt.vertices,models.bt.loads,models.bt.materials,models.bt.parameters]=DepthAverageInput(models.bt.elements,models.bt.nodes,models.bt.vertices,models.bt.loads,models.bt.materials,models.bt.parameters,VyEnum);
 
-	displaystring(m.parameters.Verbose,'\n%s',['extrude computed thickness on all layers:']);
-	results.h_g=FieldExtrude(m.elements,m.nodes,m.vertices,m.loads,m.materials,m.parameters,results.h_g,'thickness',0);
+	displaystring(verbose,'\n%s',['call computational core:']);
+	h_g=diagnostic_core_linear(models.bt,analysis_type,sub_analysis_type);
 
+	%Update
+	models.bt.elements=UpdateInputsFromSolution(models.bt.elements,models.bt.nodes,models.bt.vertices,models.bt.loads,models.bt.materials,models.bt.parameters,h_g,BalancedthicknessAnalysisEnum,NoneAnalysisEnum);
+
+	displaystring(verbose,'\n%s',['extrude computed thickness on all layers:']);
+	models.bt.elements=ExtrudeInput(models.bt.elements,models.bt.nodes,models.bt.vertices,models.bt.loads,models.bt.materials,models.bt.parameters,ThicknessEnum);
+
+	%NEED TO BE CLEANED
+	results.h_g=h_g;
 end %end function