Reshaping Nodal Spatial Vectors into 3D Matrices for Countourslice

I have a problem of temperature field across a cubic domain.

My nodes locations are stored in 3 column vectors of x=[x1 x2 ... xn]', y=[y1 y2 ... yn]', z=[z1 z2 ... zn]', while the temperature is saved also in a column vector of Temp=[T1 T2 ... Tn]'.

However, in order to use the countourslice function in Matlab, I found that the input to this function must be 3D matrices.

So I would be grateful, if I know how to reshape my vectors for being used in this function.