Index: /issm/trunk/src/m/classes/@pentaelem/CreatePVector.m
===================================================================
--- /issm/trunk/src/m/classes/@pentaelem/CreatePVector.m	(revision 945)
+++ /issm/trunk/src/m/classes/@pentaelem/CreatePVector.m	(revision 946)
@@ -657,4 +657,10 @@
 	vxvy_list= vxvyvz_list(:,1:2);
 
+	%Update material if temperature is provided
+	if temperature_is_present,
+		temperature=1/6*[1 1 1 1 1 1]*temperature_list;
+		matice.B=paterson(temperature);
+	end
+ 
 	% Get gaussian points and weights. Penta is an extrusion of a Tria, we therefore 
 	%get tria gaussian points as well as segment gaussian points. For tria gaussian 
@@ -673,8 +679,7 @@
 			gauss_coord=[first_area_gauss_coord(igarea) second_area_gauss_coord(igarea) third_area_gauss_coord(igarea) vert_gauss_coord(igvert)];
 		   
-			%Update material if temperature is provided
+			%Update temperature at gaussian point
 			if temperature_is_present,
 				temperature=GetParameterValue(pentaelem,temperature_list,gauss_coord);
-				matice.B=paterson(temperature);
 			end