Index: /issm/trunk/src/c/objects/Elements/Penta.cpp
===================================================================
--- /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 4684)
+++ /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 4685)
@@ -2642,5 +2642,8 @@
 	int          doflist[numdof];
 	double       values[numdof];
-
+	double       B[numdof];
+	double       B_average;
+
+	bool         collapse;
 	int          dummy;
 
@@ -2651,8 +2654,19 @@
 	for(i=0;i<numdof;i++){
 		values[i]=solution[doflist[i]];
+		//B[i]=Paterson(values[i]);
+	}
+
+	B_average=Paterson((values[0]+values[1]+values[2]+values[3]+values[4]+values[5])/6.0);
+	for(i=0;i<numdof;i++){
+		B[i]=B_average;
 	}
 
 	/*Add thickness as inputs to the tria element: */
 	this->inputs->AddInput(new PentaVertexInput(TemperatureEnum,values));
+
+	/*Also update the rheology WARNING: ONLY FOR PATTYN and STOKES FOR NOW*/
+	inputs->GetParameterValue(&collapse,CollapseEnum);
+	if (!collapse) this->matice->inputs->AddInput(new PentaVertexInput(RheologyBEnum,B));
+
 }
 /*}}}*/
Index: /issm/trunk/src/c/objects/Materials/Matice.h
===================================================================
--- /issm/trunk/src/c/objects/Materials/Matice.h	(revision 4684)
+++ /issm/trunk/src/c/objects/Materials/Matice.h	(revision 4685)
@@ -16,7 +16,8 @@
 	private: 
 		int	   mid;
-		Inputs*  inputs;
 
 	public:
+		/*WARNIN: input should not be public but it is an easy way to update B from T (using UpdateFromSolution) from Pentas*/
+		Inputs*  inputs;
 
 		/*Matice constructors, destructors: {{{1*/