Index: /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 16745)
+++ /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 16746)
@@ -4851,4 +4851,5 @@
 	IssmDouble  values[numdof];
 	IssmDouble  B[numdof],surface[numdof];
+	IssmDouble  B_average;
 	int        *doflist = NULL;
 	bool        hack    = false;
@@ -4969,9 +4970,9 @@
 				break;
 			case PatersonEnum:
-				for(i=0;i<numdof;i++) B[i]=Paterson(values[i]);
+				for(i=0;i<numdof;i++) B[i]=Paterson(temperatures[i]);
 				this->material->inputs->AddInput(new PentaInput(MaterialsRheologyBEnum,B,P1Enum));
 				break;
 			case ArrheniusEnum:
-				for(i=0;i<numdof;i++) B[i]=Arrhenius(values[i],surface[i]-xyz_list[i][2],material->GetN());
+				for(i=0;i<numdof;i++) B[i]=Arrhenius(temperatures[i],surface[i]-xyz_list[i][2],material->GetN());
 				this->material->inputs->AddInput(new PentaInput(MaterialsRheologyBEnum,B,P1Enum));
 				break;