function [] = map_3D(data)

clf

nul=find(data(:,1)==0);
nz=length(nul)
sz=length(data(:,1));

patch(data(1:(nul(1)-1),3),data(1:(nul(1)-1),1),data(1:(nul(1)-1),2),'y');
hold on;

for i=1:(nz-1),
  patch(data((nul(i)+1):(nul(i+1)-1),3),data((nul(i)+1):(nul(i+1)-1),1), ...
     data((nul(i)+1):(nul(i+1)-1),2),'y');
end

patch(data((nul(nz)+1):sz,3),data((nul(nz)+1):sz,1),data((nul(nz)+1):sz,2), ...
   'y');
hold off;

set(gca,'AspectRatio',[1 1]);