 load wind                                                    
 spd = sqrt(u.*u + v.*v + w.*w);                              
 p = patch(isosurface(x,y,z,spd, 40));                        
 isonormals(x,y,z,spd, p)                                     
 set(p, 'FaceColor', 'red', 'EdgeColor', 'none');             
 p2 = patch(isocaps(x,y,z,spd, 40));                          
 set(p2, 'FaceColor', 'interp', 'EdgeColor', 'none')          
 daspect([1 1 1]);                                            
 [f verts] = reducepatch(isosurface(x,y,z,spd, 30), .2);      
 colormap(jet)                                                
 axis tight; box on                                           
 camproj perspective; camva(24);                              
 campos([165 -20 65]); camtarget([100 40 -5])                 
 camlight left; lighting gouraud   