for i=1:length(allarholur)
    if allarholur{i}~[]
        toppur=allarholur{i};
        topx(i)=toppur(1,1);
        topy(i)=toppur(1,2);
        topz(i)=toppur(1,3);
    end
end

delta   = 100; 
rangeX  = -4000:delta:4000; 
rangeY  = -4000:delta:3000; 

[X,Y]   = meshgrid(rangeX,rangeY); 
Z       = griddata(topx',topy',topz',X,Y); 
mesh(X,Y,Z), hold
plot3(topx,topy,topz,'o'), hold off
% shading interp 
% view(-45,60) 
% zlim([0 6]) 
% set(gca, 'CLim', [min(zlim), max(zlim)]); 
% colorbar