Index: /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 27963)
+++ /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 27964)
@@ -3392,4 +3392,19 @@
 	inputs->SetPentaInput(enum_in,P1Enum,numindices,indices,values);
 
+}
+/*}}}*/
+void       Penta::SetElementInput(int enum_in,IssmDouble value,int type){/*{{{*/
+
+	if(type==P0Enum){
+		this->inputs->SetPentaInput(enum_in,P0Enum,this->lid,value);
+	}
+	else if(type==P1Enum){
+		IssmDouble values[6]; 
+		for(int i=0;i<6;i++)values[i]=value;
+		int lidlist[6];
+		this->GetVerticesLidList(&lidlist[0]);
+		this->inputs->SetPentaInput(enum_in,P1Enum,6,&lidlist[0],&values[0]);
+	}
+	else _error_("interpolation type not supported yet");
 }
 /*}}}*/
Index: /issm/trunk-jpl/src/c/classes/Elements/Penta.h
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Penta.h	(revision 27963)
+++ /issm/trunk-jpl/src/c/classes/Elements/Penta.h	(revision 27964)
@@ -176,5 +176,5 @@
 		void           ResetHooks();
 		void           SetElementInput(int enum_in,IssmDouble value);
-		void           SetElementInput(int enum_in,IssmDouble value,int type){_error_("not implemented yet");};
+		void           SetElementInput(int enum_in,IssmDouble value,int type);
 		void           SetElementInput(Inputs* inputs,int enum_in,IssmDouble value);
 		void           SetElementInput(Inputs* inputs,int enum_in,IssmDouble value,int type){_error_("not implemented yet");};
Index: /issm/trunk-jpl/src/c/classes/FemModel.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/FemModel.cpp	(revision 27963)
+++ /issm/trunk-jpl/src/c/classes/FemModel.cpp	(revision 27964)
@@ -2546,4 +2546,11 @@
 							InputDuplicatex(this,DamageDbarEnum,DamageDbarOldEnum);
 							this->ElementOperationx(&Element::ComputeNewDamage);
+						}
+						else if(output_enum==FrictionAlpha2Enum){
+							for(Object* & object : this->elements->objects){
+								Element* element=xDynamicCast<Element*>(object);
+								element->SetElementInput(FrictionAlpha2Enum,0.,P1Enum);
+							}
+							this->ElementOperationx(&Element::FrictionAlpha2CreateInput);
 						}