Index: /issm/trunk/src/m/solutions/balancedthickness_core.m
===================================================================
--- /issm/trunk/src/m/solutions/balancedthickness_core.m	(revision 4902)
+++ /issm/trunk/src/m/solutions/balancedthickness_core.m	(revision 4903)
@@ -13,8 +13,4 @@
 	femmodel=SetCurrentConfiguration(femmodel,BalancedthicknessAnalysisEnum);
 
-	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);
-
 	displaystring(verbose,'\n%s',['call computational core...']);
 	femmodel=solver_linear(femmodel);
Index: /issm/trunk/src/m/solutions/balancedvelocities_core.m
===================================================================
--- /issm/trunk/src/m/solutions/balancedvelocities_core.m	(revision 4902)
+++ /issm/trunk/src/m/solutions/balancedvelocities_core.m	(revision 4903)
@@ -13,11 +13,4 @@
 	femmodel=SetCurrentConfiguration(femmodel,BalancedvelocitiesAnalysisEnum);
 
-	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
-
 	displaystring(verbose,'\n%s',['call computational core...']);
 	femmodel=solver_linear(femmodel);
Index: /issm/trunk/src/m/solutions/prognostic_core.m
===================================================================
--- /issm/trunk/src/m/solutions/prognostic_core.m	(revision 4902)
+++ /issm/trunk/src/m/solutions/prognostic_core.m	(revision 4903)
@@ -12,8 +12,4 @@
 	femmodel=SetCurrentConfiguration(femmodel,PrognosticAnalysisEnum);
 
-	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);
-
 	displaystring(verbose,'\n%s',['call computational core...']);
 	femmodel=solver_linear(femmodel);