Index: /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 16738)
+++ /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 16739)
@@ -158,5 +158,5 @@
 
 	_assert_(this->material);
-	this->material->inputs->AddInput(new TriaInput(input_enum,values,interpolation_enum));
+	this->material->inputs->AddInput(new PentaInput(input_enum,values,interpolation_enum));
 }
 /*}}}*/
@@ -2021,6 +2021,9 @@
 			else if (original_input->ObjectEnum()==DoubleInputEnum)
 			 total_integrated_input=new DoubleInput(average_enum_type,0.0);
-			else
+			else{
+				printf("id %i \n",this->id);
+				printf("enum %i \n",average_enum_type);
 			 _error_("object " << EnumToStringx(original_input->ObjectEnum()) << " not supported yet");
+			}
 		}
 
@@ -4853,6 +4856,6 @@
 	IssmDouble  xyz_list[NUMVERTICES][3];
 	IssmDouble  values[numdof];
-	IssmDouble  B[numdof];
-	IssmDouble  B_average,s_average;
+	IssmDouble  B[numdof],surface[numdof];
+	IssmDouble  B_average;
 	int        *doflist = NULL;
 	bool        hack    = false;
@@ -4903,9 +4906,6 @@
 				break;
 			case ArrheniusEnum:
-				surface_input->GetInputAverage(&s_average);
-				B_average=Arrhenius((values[0]+values[1]+values[2]+values[3]+values[4]+values[5])/6.0,
-							s_average-((xyz_list[0][2]+xyz_list[1][2]+xyz_list[2][2]+xyz_list[3][2]+xyz_list[4][2]+xyz_list[5][2])/6.0),
-							material->GetN());
-				for(i=0;i<numdof;i++) B[i]=B_average;
+				GetInputListOnVertices(&surface[0],SurfaceEnum);
+				for(i=0;i<numdof;i++) B[i]=Arrhenius(values[i],surface[i]-xyz_list[i][2],material->GetN());
 				this->material->inputs->AddInput(new PentaInput(MaterialsRheologyBEnum,B,P1Enum));
 				break;
Index: /issm/trunk-jpl/src/c/classes/Elements/PentaRef.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/PentaRef.cpp	(revision 16738)
+++ /issm/trunk-jpl/src/c/classes/Elements/PentaRef.cpp	(revision 16739)
@@ -16,4 +16,5 @@
 
 /*Element macros*/
+#define NUMNODESP0    1
 #define NUMNODESP1    6
 #define NUMNODESP1_2d 3
@@ -1962,4 +1963,5 @@
 
 	switch(this->element_type){
+		case P0Enum:                return NUMNODESP0;
 		case P1Enum:                return NUMNODESP1;
 		case P1bubbleEnum:          return NUMNODESP1b;