Just realized that the contour matrix cc now seems to contain x and y location data,
pcolor(x, y, z);% Plot the data
holdon;
shadingflat
[cc, hh] = contour(x, y, z, [0 10^5],'k');% Overlay contour line
holdon;
cc(cc<1) = NaN;
h = plot(cc(1,:), cc(2, :),'w.');