Interpolation of 4D data

So I am currently working on plotting some 4D data. I have some geophyiscal data that each have their own x,y,z position in space and a corresponding data value, c. The data is in the form:

if true
x, y, z, c
37 1 1 value1
37 2 1 value2
. . . .
. . . .
. . . .
37 1 2 value_n
37 1 2 value_n2
end

Anyways I use the reshape command to put the data in the right form so that I can use the surf command:

if true
xq = reshape(x,xCols,yCols);
yq = reshape(y,xCols,yCols);
zq = reshape(z,xCols,yCols);
cq = reshape(c,xCols,yCols);
surfc(xq,yq,zq,cq,'edgecolor','none')
end

The results of the code:

How can I interpolate my values onto a finer mesh? In other words I want to interpolate my current xq, yq, zq, and cq to a finer xq, yq, zq and cq. Also is there a way to add contour lines to each cross-section?

Thanks.