Index: /issm/trunk/src/c/DataSet/DataSet.cpp
===================================================================
--- /issm/trunk/src/c/DataSet/DataSet.cpp	(revision 4063)
+++ /issm/trunk/src/c/DataSet/DataSet.cpp	(revision 4064)
@@ -685,19 +685,4 @@
 	xfree((void**)&alltruedofs);
 
-}
-/*}}}*/
-/*FUNCTION DataSet::InputExtrude(int enum_type) {{{1*/
-void  DataSet::InputExtrude(int enum_type){
-	
-	vector<Object*>::iterator object;
-	Penta* penta=NULL;
-
-	for ( object=objects.begin() ; object < objects.end(); object++ ){
-
-		if((*object)->Enum()==PentaEnum){
-			penta=(Penta*)(*object);
-			penta->InputExtrude(enum_type);
-		}
-	}
 }
 /*}}}*/
Index: /issm/trunk/src/c/DataSet/DataSet.h
===================================================================
--- /issm/trunk/src/c/DataSet/DataSet.h	(revision 4063)
+++ /issm/trunk/src/c/DataSet/DataSet.h	(revision 4064)
@@ -78,5 +78,4 @@
 		void  MeltingConstraints(int* pconverged, int* pnum_unstable_constraints);
 		DataSet* Copy(void);
-		void  InputExtrude(int enum_type);
 		int   DeleteObject(Object* object);
 		void  OutputRifts(Vec riftproperties);
Index: /issm/trunk/src/c/Makefile.am
===================================================================
--- /issm/trunk/src/c/Makefile.am	(revision 4063)
+++ /issm/trunk/src/c/Makefile.am	(revision 4064)
@@ -524,6 +524,6 @@
 					./modules/InputConvergencex/InputConvergencex.cpp\
 					./modules/InputConvergencex/InputConvergencex.h\
-					./modules/ExtrudeInputx/ExtrudeInputx.cpp\
-					./modules/ExtrudeInputx/ExtrudeInputx.h\
+					./modules/InputExtrudex/InputExtrudex.cpp\
+					./modules/InputExtrudex/InputExtrudex.h\
 					./modules/OutputRiftsx/OutputRiftsx.h\
 					./modules/OutputRiftsx/OutputRiftsx.cpp\
@@ -1038,6 +1038,6 @@
 					./modules/InputConvergencex/InputConvergencex.cpp\
 					./modules/InputConvergencex/InputConvergencex.h\
-					./modules/ExtrudeInputx/ExtrudeInputx.cpp\
-					./modules/ExtrudeInputx/ExtrudeInputx.h\
+					./modules/InputExtrudex/InputExtrudex.cpp\
+					./modules/InputExtrudex/InputExtrudex.h\
 					./modules/OutputRiftsx/OutputRiftsx.h\
 					./modules/OutputRiftsx/OutputRiftsx.cpp\
Index: /issm/trunk/src/c/modules/DepthAverageInputx/DepthAverageInputx.cpp
===================================================================
--- /issm/trunk/src/c/modules/DepthAverageInputx/DepthAverageInputx.cpp	(revision 4063)
+++ /issm/trunk/src/c/modules/DepthAverageInputx/DepthAverageInputx.cpp	(revision 4064)
@@ -4,5 +4,5 @@
 
 #include "./DepthAverageInputx.h"
-
+#include "../InputExtrudex/InputExtrudex.h"
 #include "../../shared/shared.h"
 #include "../../include/include.h"
@@ -25,5 +25,6 @@
 	}
 
+
 	/*Then extrude vertically the new inputs*/
-	elements->InputExtrude(average_enum_type);
+	InputExtrudex( elements,nodes,vertices,loads,materials,parameters,enum_type);
 }
Index: /issm/trunk/src/c/modules/InputExtrudex/InputExtrudex.cpp
===================================================================
--- /issm/trunk/src/c/modules/InputExtrudex/InputExtrudex.cpp	(revision 4064)
+++ /issm/trunk/src/c/modules/InputExtrudex/InputExtrudex.cpp	(revision 4064)
@@ -0,0 +1,27 @@
+/*!\file InputExtrudex
+ * \brief: average field throfieldh thickness
+ */
+
+#include "./InputExtrudex.h"
+#include "../../shared/shared.h"
+#include "../../include/include.h"
+#include "../../toolkits/toolkits.h"
+#include "../../EnumDefinitions/EnumDefinitions.h"
+
+void InputExtrudex( DataSet* elements,DataSet* nodes, DataSet* vertices, DataSet* loads, DataSet* materials, Parameters* parameters,int enum_type){
+
+	int i;
+
+	/*First, get elements*/
+	elements->Configure(elements,loads, nodes,vertices, materials,parameters);
+
+	/*Then extrude vertically the new inputs*/
+	for(i=0;i<elements->Size();i++){
+		Element* element=(Element*)elements->GetObjectByOffset(i);
+		if (element->Enum()==PentaEnum){
+			Penta* penta=(Penta*)element;
+			penta->InputExtrude(enum_type);
+		}
+	}
+
+}
Index: /issm/trunk/src/c/modules/InputExtrudex/InputExtrudex.h
===================================================================
--- /issm/trunk/src/c/modules/InputExtrudex/InputExtrudex.h	(revision 4064)
+++ /issm/trunk/src/c/modules/InputExtrudex/InputExtrudex.h	(revision 4064)
@@ -0,0 +1,14 @@
+/*!\file:  InputExtrudex.h
+ * \brief header file for averaging a field throfieldh thickness
+ */ 
+
+#ifndef _INPUTEXTRUDEX_H
+#define _INPUTEXTRUDEX_H
+
+#include "../../DataSet/DataSet.h"
+
+/* local prototypes: */
+void InputExtrudex( DataSet* elements,DataSet* nodes, DataSet* vertices, DataSet* loads, DataSet* materials, Parameters* parameters,int enum_type);
+
+#endif  /* _INPUTEXTRUDEX_H */
+
Index: /issm/trunk/src/c/modules/modules.h
===================================================================
--- /issm/trunk/src/c/modules/modules.h	(revision 4063)
+++ /issm/trunk/src/c/modules/modules.h	(revision 4064)
@@ -63,5 +63,5 @@
 #include "./SplitSolutionVectorx/SplitSolutionVectorx.h"
 #include "./DepthAverageInputx/DepthAverageInputx.h"
-#include "./ExtrudeInputx/ExtrudeInputx.h"
+#include "./InputExtrudex/InputExtrudex.h"
 #include "./InputToResultx/InputToResultx.h"
 #include "./GetSolutionFromInputsx/GetSolutionFromInputsx.h"
Index: /issm/trunk/src/c/solutions/balancedthickness_core.cpp
===================================================================
--- /issm/trunk/src/c/solutions/balancedthickness_core.cpp	(revision 4063)
+++ /issm/trunk/src/c/solutions/balancedthickness_core.cpp	(revision 4064)
@@ -34,5 +34,5 @@
 
 	_printf_("extrude computed thickness on all layers:\n");
-	ExtrudeInputx(femmodel->elements,femmodel->nodes, femmodel->vertices,femmodel->loads, femmodel->materials,femmodel->parameters,ThicknessEnum);
+	InputExtrudex(femmodel->elements,femmodel->nodes, femmodel->vertices,femmodel->loads, femmodel->materials,femmodel->parameters,ThicknessEnum);
 
 	if(verbose)_printf_("saving results:\n");
Index: /issm/trunk/src/c/solutions/balancedvelocities_core.cpp
===================================================================
--- /issm/trunk/src/c/solutions/balancedvelocities_core.cpp	(revision 4063)
+++ /issm/trunk/src/c/solutions/balancedvelocities_core.cpp	(revision 4064)
@@ -33,6 +33,6 @@
 
 	_printf_("extrude computed velocity on all layers:\n");
-	ExtrudeInputx(femmodel->elements,femmodel->nodes, femmodel->vertices,femmodel->loads, femmodel->materials,femmodel->parameters,VxEnum);
-	ExtrudeInputx(femmodel->elements,femmodel->nodes, femmodel->vertices,femmodel->loads, femmodel->materials,femmodel->parameters,VyEnum);
+	InputExtrudex(femmodel->elements,femmodel->nodes, femmodel->vertices,femmodel->loads, femmodel->materials,femmodel->parameters,VxEnum);
+	InputExtrudex(femmodel->elements,femmodel->nodes, femmodel->vertices,femmodel->loads, femmodel->materials,femmodel->parameters,VyEnum);
 
 	if(verbose)_printf_("saving results:\n");
Index: /issm/trunk/src/c/solutions/bedslope_core.cpp
===================================================================
--- /issm/trunk/src/c/solutions/bedslope_core.cpp	(revision 4063)
+++ /issm/trunk/src/c/solutions/bedslope_core.cpp	(revision 4064)
@@ -33,6 +33,6 @@
 	if (dim==3){
 		if(verbose)_printf_("%s\n","extruding bed in 3d...");
-		ExtrudeInputx(femmodel->elements,femmodel->nodes, femmodel->vertices,femmodel->loads, femmodel->materials,femmodel->parameters,BedSlopeXEnum);
-		ExtrudeInputx(femmodel->elements,femmodel->nodes, femmodel->vertices,femmodel->loads, femmodel->materials,femmodel->parameters,BedSlopeYEnum);
+		InputExtrudex(femmodel->elements,femmodel->nodes, femmodel->vertices,femmodel->loads, femmodel->materials,femmodel->parameters,BedSlopeXEnum);
+		InputExtrudex(femmodel->elements,femmodel->nodes, femmodel->vertices,femmodel->loads, femmodel->materials,femmodel->parameters,BedSlopeYEnum);
 	}
 	
Index: /issm/trunk/src/c/solutions/prognostic_core.cpp
===================================================================
--- /issm/trunk/src/c/solutions/prognostic_core.cpp	(revision 4063)
+++ /issm/trunk/src/c/solutions/prognostic_core.cpp	(revision 4064)
@@ -30,5 +30,5 @@
 		
 	_printf_("extrude computed thickness on all layers:\n");
-	ExtrudeInputx(femmodel->elements,femmodel->nodes, femmodel->vertices,femmodel->loads, femmodel->materials,femmodel->parameters,ThicknessEnum);
+	InputExtrudex(femmodel->elements,femmodel->nodes, femmodel->vertices,femmodel->loads, femmodel->materials,femmodel->parameters,ThicknessEnum);
 
 	if(verbose)_printf_("saving results:\n");
Index: /issm/trunk/src/c/solutions/surfaceslope_core.cpp
===================================================================
--- /issm/trunk/src/c/solutions/surfaceslope_core.cpp	(revision 4063)
+++ /issm/trunk/src/c/solutions/surfaceslope_core.cpp	(revision 4064)
@@ -33,6 +33,6 @@
 	if (dim==3){
 		if(verbose)_printf_("%s\n","extruding slope in 3d...");
-		ExtrudeInputx(femmodel->elements,femmodel->nodes, femmodel->vertices,femmodel->loads, femmodel->materials,femmodel->parameters,SurfaceSlopeXEnum);
-		ExtrudeInputx(femmodel->elements,femmodel->nodes, femmodel->vertices,femmodel->loads, femmodel->materials,femmodel->parameters,SurfaceSlopeYEnum);
+		InputExtrudex(femmodel->elements,femmodel->nodes, femmodel->vertices,femmodel->loads, femmodel->materials,femmodel->parameters,SurfaceSlopeXEnum);
+		InputExtrudex(femmodel->elements,femmodel->nodes, femmodel->vertices,femmodel->loads, femmodel->materials,femmodel->parameters,SurfaceSlopeYEnum);
 	}
 	
Index: /issm/trunk/src/c/solutions/transient2d.cpp
===================================================================
--- /issm/trunk/src/c/solutions/transient2d.cpp	(revision 4063)
+++ /issm/trunk/src/c/solutions/transient2d.cpp	(revision 4064)
@@ -61,5 +61,5 @@
 
 	_printf_("create finite element model:\n");
-	femmodel=new FemModel(fid,solution_type,analyses,7);
+	femmodel=new FemModel(fid,solution_type,analyses,2);
 
 	/*add outputfilename in parameters: */
Index: /issm/trunk/src/c/solvers/solver_diagnostic_nonlinear.cpp
===================================================================
--- /issm/trunk/src/c/solvers/solver_diagnostic_nonlinear.cpp	(revision 4063)
+++ /issm/trunk/src/c/solvers/solver_diagnostic_nonlinear.cpp	(revision 4064)
@@ -137,6 +137,6 @@
 	if (dim==3){
 		if(verbose)_printf_("%s\n","extruding velocity and pressure in 3d...");
-		ExtrudeInputx( fem->elements,fem->nodes,fem->vertices,fem->loads,fem->materials,fem->parameters,VxEnum);
-		ExtrudeInputx( fem->elements,fem->nodes,fem->vertices,fem->loads,fem->materials,fem->parameters,VyEnum);
+		InputExtrudex( fem->elements,fem->nodes,fem->vertices,fem->loads,fem->materials,fem->parameters,VxEnum);
+		InputExtrudex( fem->elements,fem->nodes,fem->vertices,fem->loads,fem->materials,fem->parameters,VyEnum);
 	}