Index: /issm/trunk/src/c/objects/Matice.cpp
===================================================================
--- /issm/trunk/src/c/objects/Matice.cpp	(revision 944)
+++ /issm/trunk/src/c/objects/Matice.cpp	(revision 945)
@@ -375,4 +375,5 @@
 		}
 	}
+
 	#ifdef _DEBUG_
 	_printf_("Viscosity %lf\n",viscosity3d);
Index: /issm/trunk/src/c/objects/Penta.cpp
===================================================================
--- /issm/trunk/src/c/objects/Penta.cpp	(revision 944)
+++ /issm/trunk/src/c/objects/Penta.cpp	(revision 945)
@@ -1102,4 +1102,11 @@
 	//Update material if necessary
 	if(inputs->Recover("temperature_average",&temperature_list[0],1,dofs,6,(void**)nodes)){
+		if(matice){
+			temperature_average=(temperature_list[0]+temperature_list[1]+temperature_list[2]+temperature_list[3]+temperature_list[4]+temperature_list[5])/6.0;
+			B_average=Paterson(temperature_average);
+			matice->SetB(B_average);
+		}
+	}
+	if(inputs->Recover("temperature",&temperature_list[0],1,dofs,6,(void**)nodes)){
 		if(matice){
 			temperature_average=(temperature_list[0]+temperature_list[1]+temperature_list[2]+temperature_list[3]+temperature_list[4]+temperature_list[5])/6.0;
@@ -3755,5 +3762,4 @@
 	if(onbed && shelf){
 
-
 		tria=(Tria*)SpawnTria(0,1,2); //grids 0, 1 and 2 make the new tria.
 		tria->CreatePVectorThermalShelf(pg,inputs, analysis_type,sub_analysis_type);