How do I plot a contour map on the surface of a 3-d surface?

The idea is similar to using "surf(X,Y,Z,C)" but I want contour lines in addition to the coloring provided by C. The contour lines need to correspond to f(x,y,z), not x, y, or z. If there is a non-toolbox specific answer, that would be great. Thanks.

Answer 1

[Mike Garrity]

A simple tricky I use sometimes is to use a colormap with a small number of colors.

[x,y] = meshgrid(linspace(-pi,pi,150));
z = cos(x).*cos(y);
c = sin(x).*sin(y);
surf(x,y,z,c,'FaceColor','interp','EdgeColor','none')
zlabel('cosine')
colormap(lines(7))
c = colorbar;
c.Label.String = 'sine';
view([-45 60])

This makes my CData variable look like what the contourf function would draw.

[Joseph]

Not exactly what I was looking for but it meets the intent. Thanks!

[John Mickett]

I would like to put contour lines OVER a surf plot (two different variables), so that solution doesn't work. Any suggestions? I do it in 2D all the time.

[GiWon Hong]

Thanks, It is very helpful.

I used colormap(jet(31)) instead of colormap(lines(7)). It gives more contour likely view.