hold on
anomaly=zeros(34,35,36);
for j=1:34;
   for jj=1:36;
      for jjj=1:16;
    anomaly(j,jj,jjj)=1000;
       end
   end
end

for radius=1:.1:10
theta=linspace(0,2*pi);
phi=linspace(0,2*pi);
[theta,phi] = meshgrid(theta,phi);
[xs,ys,zs] = sph2cart(theta,phi,radius);
surf(xs,ys,zs);
end

[x,y,z]=meshgrid(1:35,1:34,1:16);
xslice = [1  34 ]; yslice = 15; zslice = [1  16];
slice(x,y,z,anomaly,xslice,yslice,zslice)
axis off