Index: /issm/trunk/src/m/classes/@triaelem/CreateKMatrix.m
===================================================================
--- /issm/trunk/src/m/classes/@triaelem/CreateKMatrix.m	(revision 894)
+++ /issm/trunk/src/m/classes/@triaelem/CreateKMatrix.m	(revision 895)
@@ -278,5 +278,13 @@
 	end
 	if(temperature_average_is_present) temperature_average_list(i)=temperature_average_param(dof);end;
-end
+
+end
+
+%Update material if temperature is provided.
+if temperature_average_is_present,
+	temperature_average=1/3*([1 1 1]*temperature_average_list);
+	matice.B=paterson(temperature_average);
+end
+
 if (element_debug & triaelem.id==element_debugid),
 	disp(sprintf('El id %i TriaElemnet input list before gaussian loop: \n',element_debugid)); 
@@ -381,10 +389,4 @@
 	%disp(sprintf('Epsilon: %f %f %f\n',epsilon(1),epsilon(2),epsilon(3)));
 	
-	%Update material if temperature is provided.
-	if temperature_average_is_present,
-		temperature_average=GetParameterValue(triaelem,temperature_average_list,gauss_l1l2l3);
-		matice.B=paterson(temperature_average);
-	end
-
 	%Update material if flow law is specified. This will erase the previous change 
 	%on B when temperature is provided.