load wind
figure(1)
zmax = max(z(:)); zmin = min(z(:));
h=streamslice(x,y,z,u,v,w,[],[],(zmax-zmin)/2);
set(h,'Color','black','LineWidth',2);
figure(2)
[sx,sy,sz] = meshgrid(80,20:10:50,0:5:15);
plot3(sx(:),sy(:),sz(:),'*r');
axis(volumebounds(x,y,z,u,v,w))
grid on; box; daspect([2 2 1])
h2=streamline(x,y,z,u,v,w,sx(:),sy(:),sz(:));
set(h2,'Color','black','LineWidth',2);














