clear all;
close all;
load('mb15.mat');
load('mc15.mat');

mcX = mb15.x(1500:1700);
mcY = mb15.y(1500:1700);
mcZ = mb15.z(1500:1700, 1500:1700);

V = [0:10:100 100:500:500];
V = -V;
fprintf('Starting contourf\n');
[C, hF] = contour3(mcX, mcY, mcZ, V);
fprintf('contour done \n');

surf(mcX, mcY, mcZ, 'EdgeColor', 'none')
daspect
colormap(haxby(256));
daspect([8,8,1])

hold;

fprintf('Starting contourf\n');
[C, hF] = contour3(mcX, mcY, mcZ, V, '-k');
fprintf('contour done \n');

radius = 300
x = -radius:50:radius;
y = sqrt(radius^2 - x.^2);
x = x + 606000;
y = y + 4073500;
z = 1:length(y);
z = -5 .* z
plot3(x, y, z, 'or', 'MarkerFaceColor', 'r');
