Index: /issm/trunk/src/c/solutions/transient3d_core.cpp
===================================================================
--- /issm/trunk/src/c/solutions/transient3d_core.cpp	(revision 4746)
+++ /issm/trunk/src/c/solutions/transient3d_core.cpp	(revision 4747)
@@ -50,4 +50,7 @@
 		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/m/solutions/transient3d_core.m
===================================================================
--- /issm/trunk/src/m/solutions/transient3d_core.m	(revision 4746)
+++ /issm/trunk/src/m/solutions/transient3d_core.m	(revision 4747)
@@ -31,4 +31,7 @@
 		[femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters]=InputDepthAverage(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,TemperatureEnum,TemperatureAverageEnum);
 
+		displaystring(verbose,'\n%s',['Update material from temperature...']);
+		[femmodel.elements,femmodel.materials]=MaterialUpdateFromTemperature(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters);
+
 		displaystring(verbose,'\n%s',['computing new velocities...']);
 		femmodel=diagnostic_core(femmodel);