load windxmin = min(x(:));xmax = max(x(:));ymin = min(y(:));ymax = max(y(:));zmin = min(z(:));daspect([2,2,1])xrange = linspace(xmin,xmax,8);yrange = linspace(ymin,ymax,8);zrange = 3:4:15;[cx cy cz] = meshgrid(xrange,yrange,zrange);hcones = coneplot(x,y,z,u,v,w,cx,cy,cz,5);set(hcones,'FaceColor','red','EdgeColor','none')hold onwind_speed = sqrt(u.^2 + v.^2 + w.^2);hsurfaces = slice(x,y,z,wind_speed,[xmin,xmax],ymax,zmin);set(hsurfaces,'FaceColor','interp','EdgeColor','none')hold offaxis tight; view(30,40); axis offcamproj perspective; camzoom(1.5)camlight right; lighting phongset(hsurfaces,'AmbientStrength',.6)set(hcones,'DiffuseStrength',.8) --------------------simwe
|
