Index: /issm/trunk/src/c/EnumDefinitions/EnumAsString.cpp
===================================================================
--- /issm/trunk/src/c/EnumDefinitions/EnumAsString.cpp	(revision 4934)
+++ /issm/trunk/src/c/EnumDefinitions/EnumAsString.cpp	(revision 4935)
@@ -197,5 +197,4 @@
 		case TemperatureEnum : return "Temperature";
 		case TemperatureOldEnum : return "TemperatureOld";
-		case TemperatureAverageEnum : return "TemperatureAverage";
 		case ThicknessEnum : return "Thickness";
 		case TypeEnum : return "Type";
Index: /issm/trunk/src/c/EnumDefinitions/EnumDefinitions.h
===================================================================
--- /issm/trunk/src/c/EnumDefinitions/EnumDefinitions.h	(revision 4934)
+++ /issm/trunk/src/c/EnumDefinitions/EnumDefinitions.h	(revision 4935)
@@ -224,5 +224,4 @@
 	TemperatureEnum,
 	TemperatureOldEnum,
-	TemperatureAverageEnum,
 	ThicknessEnum,
 	TypeEnum,
Index: /issm/trunk/src/c/EnumDefinitions/StringAsEnum.cpp
===================================================================
--- /issm/trunk/src/c/EnumDefinitions/StringAsEnum.cpp	(revision 4934)
+++ /issm/trunk/src/c/EnumDefinitions/StringAsEnum.cpp	(revision 4935)
@@ -195,5 +195,4 @@
 	else if (strcmp(name,"Temperature")==0) return TemperatureEnum;
 	else if (strcmp(name,"TemperatureOld")==0) return TemperatureOldEnum;
-	else if (strcmp(name,"TemperatureAverage")==0) return TemperatureAverageEnum;
 	else if (strcmp(name,"Thickness")==0) return ThicknessEnum;
 	else if (strcmp(name,"Type")==0) return TypeEnum;
Index: /issm/trunk/src/c/Makefile.am
===================================================================
--- /issm/trunk/src/c/Makefile.am	(revision 4934)
+++ /issm/trunk/src/c/Makefile.am	(revision 4935)
@@ -554,6 +554,4 @@
 					./modules/InputExtrudex/InputExtrudex.cpp\
 					./modules/InputExtrudex/InputExtrudex.h\
-					./modules/MaterialUpdateFromTemperaturex/MaterialUpdateFromTemperaturex.cpp\
-					./modules/MaterialUpdateFromTemperaturex/MaterialUpdateFromTemperaturex.h\
 					./modules/OutputRiftsx/OutputRiftsx.h\
 					./modules/OutputRiftsx/OutputRiftsx.cpp\
@@ -1097,6 +1095,4 @@
 					./modules/InputExtrudex/InputExtrudex.cpp\
 					./modules/InputExtrudex/InputExtrudex.h\
-					./modules/MaterialUpdateFromTemperaturex/MaterialUpdateFromTemperaturex.cpp\
-					./modules/MaterialUpdateFromTemperaturex/MaterialUpdateFromTemperaturex.h\
 					./modules/OutputRiftsx/OutputRiftsx.h\
 					./modules/OutputRiftsx/OutputRiftsx.cpp\
Index: /issm/trunk/src/c/modules/modules.h
===================================================================
--- /issm/trunk/src/c/modules/modules.h	(revision 4934)
+++ /issm/trunk/src/c/modules/modules.h	(revision 4935)
@@ -43,5 +43,4 @@
 #include "./InputUpdateFromVectorx/InputUpdateFromVectorx.h"
 #include "./MassFluxx/MassFluxx.h"
-#include "./MaterialUpdateFromTemperaturex/MaterialUpdateFromTemperaturex.h"
 #include "./MaxAbsVxx/MaxAbsVxx.h"
 #include "./MaxAbsVyx/MaxAbsVyx.h"
Index: /issm/trunk/src/c/objects/Elements/Element.h
===================================================================
--- /issm/trunk/src/c/objects/Elements/Element.h	(revision 4934)
+++ /issm/trunk/src/c/objects/Elements/Element.h	(revision 4935)
@@ -54,5 +54,4 @@
 		virtual void   UpdateGeometry(void)=0;
 		virtual void   InputToResult(int enum_type,int step,double time)=0;
-		virtual void   MaterialUpdateFromTemperature(void)=0;
 		virtual void   ProcessResultsUnits(void)=0;
 		virtual void   MinVel(double* pminvel, bool process_units)=0;
Index: /issm/trunk/src/c/objects/Elements/Penta.cpp
===================================================================
--- /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 4934)
+++ /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 4935)
@@ -1183,36 +1183,4 @@
 double Penta::MassFlux( double* segment){
 	ISSMERROR(" not supported yet!");
-}
-/*}}}*/
-/*FUNCTION Penta::MaterialUpdateFromTemperature{{{1*/
-void  Penta::MaterialUpdateFromTemperature(void){
-
-	/*Intermediaries*/
-	double Taverage;
-	double B;
-
-	/*inputs: */
-	bool   collapse;
-	Input *TempInput = NULL;
-	Input *RheologyBInput = NULL;
-
-	/*retrieve inputs :*/
-	inputs->GetParameterValue(&collapse,CollapseEnum);
-
-	if (collapse){
-		/*find temperature average*/
-		this->inputs->GetParameterAverage(&Taverage,TemperatureAverageEnum);
-	}
-	else{
-		/*find 3d temperature*/
-		this->inputs->GetParameterAverage(&Taverage,TemperatureEnum);
-	}
-
-	/*Get B from T using Paterson*/
-	B=Paterson(Taverage);
-
-	/*Add input to Matice*/
-	this->matice->inputs->AddInput(new DoubleInput(RheologyBEnum,B));
-
 }
 /*}}}*/
@@ -5342,5 +5310,4 @@
 				name==VzEnum ||
 				name==TemperatureEnum ||
-				name==TemperatureAverageEnum ||
 				name==RheologyBEnum ||
 				name==RheologyNEnum ||
Index: /issm/trunk/src/c/objects/Elements/Penta.h
===================================================================
--- /issm/trunk/src/c/objects/Elements/Penta.h	(revision 4934)
+++ /issm/trunk/src/c/objects/Elements/Penta.h	(revision 4935)
@@ -93,5 +93,4 @@
 		void   InputScale(int enum_type,double scale_factor);
 		void   InputToResult(int enum_type,int step,double time);
-		void   MaterialUpdateFromTemperature(void);
 		double MassFlux(double* segment);
 		void   MaxAbsVx(double* pmaxabsvx, bool process_units);
Index: /issm/trunk/src/c/solutions/steadystate_core.cpp
===================================================================
--- /issm/trunk/src/c/solutions/steadystate_core.cpp	(revision 4934)
+++ /issm/trunk/src/c/solutions/steadystate_core.cpp	(revision 4935)
@@ -35,10 +35,4 @@
 		thermal_core(femmodel);
 
-		if(verbose)_printf_("%s\n","computing depth average temperature");
-		InputDepthAveragex(femmodel->elements,femmodel->nodes, femmodel->vertices,femmodel->loads, femmodel->materials,femmodel->parameters,TemperatureEnum,TemperatureAverageEnum);
-	
-		if(verbose)_printf_("%s\n","update B from temperature");
-		MaterialUpdateFromTemperaturex(femmodel->elements,femmodel->nodes, femmodel->vertices,femmodel->loads, femmodel->materials,femmodel->parameters);
-
 		if(verbose)_printf_("%s\n","computing new velocity");
 		diagnostic_core(femmodel);
Index: /issm/trunk/src/c/solutions/transient3d_core.cpp
===================================================================
--- /issm/trunk/src/c/solutions/transient3d_core.cpp	(revision 4934)
+++ /issm/trunk/src/c/solutions/transient3d_core.cpp	(revision 4935)
@@ -47,10 +47,4 @@
 		thermal_core_step(femmodel,step,time);
 
-		if(verbose)_printf_("%s\n","computing depth average temperature");
-		InputDepthAveragex(femmodel->elements,femmodel->nodes, femmodel->vertices,femmodel->loads, femmodel->materials,femmodel->parameters,TemperatureEnum,TemperatureAverageEnum);
-
-		if(verbose)_printf_("%s\n","update B from temperature");
-		MaterialUpdateFromTemperaturex(femmodel->elements,femmodel->nodes, femmodel->vertices,femmodel->loads, femmodel->materials,femmodel->parameters);
-
 		if(verbose)_printf_("%s\n","computing new velocity");
 		diagnostic_core(femmodel);