load wind                                       
%打開保存的數據
lims = [100.64 116.67 17.25 28.75 -0.02 6.86];         
%定義坐標范圍
[x,y,z,u,v,w] = subvolume(x,y,z,u,v,w,lims);            
%lims來定義數據子集
[sx sy sz] = meshgrid(110,20:5:30,1:5);               
%定義網格點
verts = stream3(x,y,z,u,v,w,sx,sy,sz,.5);               
%計算彩帶頂點
cav = curl(x,y,z,u,v,w);                            
%計算卷曲角速度
wind_speed = sqrt(u.^2 + v.^2 + w.^2);               
%計算流速
h = streamribbon(verts,x,y,z,cav,wind_speed,2);        
%繪制流彩帶圖
view(3)                                        
%三視圖