Index: /issm/trunk/src/c/objects/Elements/Penta.cpp
===================================================================
--- /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 4273)
+++ /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 4274)
@@ -56,15 +56,15 @@
 
 	if isnan(iomodel->upperelements[index]){
-		penta_elements_ids[0]=this->id; //upper penta is the same penta
+		penta_elements_ids[1]=this->id; //upper penta is the same penta
 	}
 	else{
-		penta_elements_ids[0]=(int)(iomodel->upperelements[index]);
+		penta_elements_ids[1]=(int)(iomodel->upperelements[index]);
 	}
 	
 	if isnan(iomodel->lowerelements[index]){
-		penta_elements_ids[1]=this->id; //lower penta is the same penta
+		penta_elements_ids[0]=this->id; //lower penta is the same penta
 	}
 	else{
-		penta_elements_ids[1]=(int)(iomodel->lowerelements[index]);
+		penta_elements_ids[0]=(int)(iomodel->lowerelements[index]);
 	}
 	this->InitHookNodes(nummodels);this->nodes=NULL;
@@ -4141,6 +4141,7 @@
 	Input* original_input=NULL;
 
-	/*Make a copy of the original input: */
+	/*Extrude the original input: */
 	original_input=(Input*)this->inputs->GetInput(enum_type);
+	original_input->Extrude();
 
 	/*recover parameters: */
Index: /issm/trunk/src/c/objects/Inputs/BeamVertexInput.h
===================================================================
--- /issm/trunk/src/c/objects/Inputs/BeamVertexInput.h	(revision 4273)
+++ /issm/trunk/src/c/objects/Inputs/BeamVertexInput.h	(revision 4274)
@@ -69,4 +69,5 @@
 		void AXPY(Input* xinput,double scalar);
 		void Constrain(double cm_min, double cm_max);
+		void Extrude(void){ISSMERROR("not supported yet");};
 		void GetVectorFromInputs(Vec vector,int* doflist);
 		void GetValuesPtr(double** pvalues,int* pnum_values);
Index: /issm/trunk/src/c/objects/Inputs/BoolInput.h
===================================================================
--- /issm/trunk/src/c/objects/Inputs/BoolInput.h	(revision 4273)
+++ /issm/trunk/src/c/objects/Inputs/BoolInput.h	(revision 4274)
@@ -69,4 +69,5 @@
 		void AXPY(Input* xinput,double scalar);
 		void Constrain(double cm_min, double cm_max);
+		void Extrude(void){ISSMERROR("not supported yet");};
 		void GetVectorFromInputs(Vec vector,int* doflist);
 		void GetValuesPtr(double** pvalues,int* pnum_values);
Index: /issm/trunk/src/c/objects/Inputs/DoubleInput.h
===================================================================
--- /issm/trunk/src/c/objects/Inputs/DoubleInput.h	(revision 4273)
+++ /issm/trunk/src/c/objects/Inputs/DoubleInput.h	(revision 4274)
@@ -68,4 +68,5 @@
 		void AXPY(Input* xinput,double scalar);
 		void Constrain(double cm_min, double cm_max);
+		void Extrude(void){ISSMERROR("not supported yet");};
 		void GetVectorFromInputs(Vec vector,int* doflist);
 		void GetValuesPtr(double** pvalues,int* pnum_values);
Index: /issm/trunk/src/c/objects/Inputs/Input.h
===================================================================
--- /issm/trunk/src/c/objects/Inputs/Input.h	(revision 4273)
+++ /issm/trunk/src/c/objects/Inputs/Input.h	(revision 4274)
@@ -52,4 +52,5 @@
 		virtual void AXPY(Input* xinput,double scalar)=0;
 		virtual void Constrain(double cm_min, double cm_max)=0;
+		virtual void Extrude()=0;
 		virtual void GetVectorFromInputs(Vec vector,int* doflist)=0;
 		virtual void GetValuesPtr(double** pvalues,int* pnum_values)=0;
Index: /issm/trunk/src/c/objects/Inputs/IntInput.h
===================================================================
--- /issm/trunk/src/c/objects/Inputs/IntInput.h	(revision 4273)
+++ /issm/trunk/src/c/objects/Inputs/IntInput.h	(revision 4274)
@@ -69,4 +69,5 @@
 		void AXPY(Input* xinput,double scalar);
 		void Constrain(double cm_min, double cm_max);
+		void Extrude(void){ISSMERROR("not supported yet");};
 		void GetVectorFromInputs(Vec vector,int* doflist);
 		void GetValuesPtr(double** pvalues,int* pnum_values);
Index: /issm/trunk/src/c/objects/Inputs/PentaVertexInput.cpp
===================================================================
--- /issm/trunk/src/c/objects/Inputs/PentaVertexInput.cpp	(revision 4273)
+++ /issm/trunk/src/c/objects/Inputs/PentaVertexInput.cpp	(revision 4274)
@@ -944,4 +944,14 @@
 }
 /*}}}*/
+/*FUNCTION PentaVertexInput::Extrude(double cm_min, double cm_max){{{1*/
+void PentaVertexInput::Extrude(void){
+
+	int i;
+
+	/*First 3 values copied on 3 last values*/
+	for(i=0;i<3;i++) this->values[3+i]=this->values[i];
+
+}
+/*}}}*/
 /*FUNCTION PentaVertexInput::GetVectorFromInputs(Vec vector,int* doflist){{{1*/
 void PentaVertexInput::GetVectorFromInputs(Vec vector,int* doflist){
Index: /issm/trunk/src/c/objects/Inputs/PentaVertexInput.h
===================================================================
--- /issm/trunk/src/c/objects/Inputs/PentaVertexInput.h	(revision 4273)
+++ /issm/trunk/src/c/objects/Inputs/PentaVertexInput.h	(revision 4274)
@@ -78,4 +78,5 @@
 		void AXPY(Input* xinput,double scalar);
 		void Constrain(double cm_min, double cm_max);
+		void Extrude(void);
 		void GetVectorFromInputs(Vec vector,int* doflist);
 		void GetValuesPtr(double** pvalues,int* pnum_values);
Index: /issm/trunk/src/c/objects/Inputs/SingVertexInput.h
===================================================================
--- /issm/trunk/src/c/objects/Inputs/SingVertexInput.h	(revision 4273)
+++ /issm/trunk/src/c/objects/Inputs/SingVertexInput.h	(revision 4274)
@@ -68,4 +68,5 @@
 		void AXPY(Input* xinput,double scalar);
 		void Constrain(double cm_min, double cm_max);
+		void Extrude(void){ISSMERROR("not supported yet");};
 		void GetVectorFromInputs(Vec vector,int* doflist);
 		void GetValuesPtr(double** pvalues,int* pnum_values);
Index: /issm/trunk/src/c/objects/Inputs/TriaVertexInput.h
===================================================================
--- /issm/trunk/src/c/objects/Inputs/TriaVertexInput.h	(revision 4273)
+++ /issm/trunk/src/c/objects/Inputs/TriaVertexInput.h	(revision 4274)
@@ -76,4 +76,5 @@
 		void AXPY(Input* xinput,double scalar);
 		void Constrain(double cm_min, double cm_max);
+		void Extrude(void){ISSMERROR("not supported yet");};
 		void GetVectorFromInputs(Vec vector,int* doflist);
 		void GetValuesPtr(double** pvalues,int* pnum_values);