Index: /issm/trunk/src/c/solutions/balancedthickness_core.cpp
===================================================================
--- /issm/trunk/src/c/solutions/balancedthickness_core.cpp	(revision 4901)
+++ /issm/trunk/src/c/solutions/balancedthickness_core.cpp	(revision 4902)
@@ -27,8 +27,4 @@
 	femmodel->parameters->FindParam(&solution_type,SolutionTypeEnum);
 
-	_printf_("depth averaging velocity...\n");
-	InputDepthAveragex(femmodel->elements,femmodel->nodes, femmodel->vertices,femmodel->loads, femmodel->materials,femmodel->parameters,VxEnum,VxAverageEnum);
-	InputDepthAveragex(femmodel->elements,femmodel->nodes, femmodel->vertices,femmodel->loads, femmodel->materials,femmodel->parameters,VyEnum,VyAverageEnum);
-
 	_printf_("call computational core:\n");
 	solver_linear(femmodel);
Index: /issm/trunk/src/c/solutions/balancedvelocities_core.cpp
===================================================================
--- /issm/trunk/src/c/solutions/balancedvelocities_core.cpp	(revision 4901)
+++ /issm/trunk/src/c/solutions/balancedvelocities_core.cpp	(revision 4902)
@@ -26,9 +26,4 @@
 	femmodel->parameters->FindParam(&solution_type,SolutionTypeEnum);
 
-	_printf_("depth averaging velocity...\n");
-	InputDepthAveragex(femmodel->elements,femmodel->nodes, femmodel->vertices,femmodel->loads, femmodel->materials,femmodel->parameters,VxEnum,VxAverageEnum);
-	InputDepthAveragex(femmodel->elements,femmodel->nodes, femmodel->vertices,femmodel->loads, femmodel->materials,femmodel->parameters,VyEnum,VyAverageEnum);
-	if(dim==3) InputDepthAveragex(femmodel->elements,femmodel->nodes, femmodel->vertices,femmodel->loads, femmodel->materials,femmodel->parameters,VzEnum,VzAverageEnum);
-
 	_printf_("call computational core:\n");
 	solver_linear(femmodel);
Index: /issm/trunk/src/c/solutions/prognostic_core.cpp
===================================================================
--- /issm/trunk/src/c/solutions/prognostic_core.cpp	(revision 4901)
+++ /issm/trunk/src/c/solutions/prognostic_core.cpp	(revision 4902)
@@ -24,8 +24,4 @@
 	femmodel->parameters->FindParam(&solution_type,SolutionTypeEnum);
 
-	_printf_("depth averaging velocity...\n");
-	InputDepthAveragex(femmodel->elements,femmodel->nodes, femmodel->vertices,femmodel->loads, femmodel->materials,femmodel->parameters,VxEnum,VxAverageEnum);
-	InputDepthAveragex(femmodel->elements,femmodel->nodes, femmodel->vertices,femmodel->loads, femmodel->materials,femmodel->parameters,VyEnum,VyAverageEnum);
-	
 	_printf_("call computational core:\n");
 	solver_linear(femmodel);