台风大涡模拟三维可视化

MeteoInfo

针对WRF_LES台风大涡模拟数据，提取风场U, V, W分量，计算出三维风速场，用volumeplot函数实现风速三维体可视化。

1. fn = r'D:\Temp\LaSW\airship\WRF_LES\wrfout_d05_2014-07-20_05_00_00'
   
2. f = addfile(fn)
   
3. u = f['U'][0]
   
4. v = f['V'][0]
   
5. w = f['W'][0]
   
6. u = meteo.wrf.destagger(u, -1)
   
7. v = meteo.wrf.destagger(v, -2)
   
8. w = meteo.wrf.destagger(w, -3)
   
9. speed = sqrt(u*u + v*v + w*w)
   
10. nz, ny, nx = u.shape
    
11. height = arange(nz)
    
12. y = u.dimvalue(1)
    
13. x = u.dimvalue(2)
    
14. 
    
15. print('Plot...')
    
16. figure(facecolor='k', newfig=False)
    
17. ax = axes3d(projection=f.proj, clip_plane=False, orthographic=False, axis=False)
    
18. volumeplot(x, y, height, speed, cmap='BlAqGrYeOrRe', ray_casting='specular',
    
19.     brightness=1.5, vmin=10, vmax=50,
    
20.     opacity_nodes=[10,20,30,35,40,50], opacity_levels=[0,0.05,0.1,0.2,0.8,1])
    
21. colorbar(aspect=30, tickcolor='w')
    
22. xlim(x[0], x[-1])
    
23. ylim(y[0], y[-1])

fn = r'D:\Temp\LaSW\airship\WRF_LES\wrfout_d05_2014-07-20_05_00_00'<br /> f = addfile(fn)<br /> u = f['U'][0]<br /> v = f['V'][0]<br /> w = f['W'][0]<br /> u = meteo.wrf.destagger(u, -1)<br /> v = meteo.wrf.destagger(v, -2)<br /> w = meteo.wrf.destagger(w, -3)<br /> speed = sqrt(u*u + v*v + w*w)<br /> nz, ny, nx = u.shape<br /> height = arange(nz)<br /> y = u.dimvalue(1)<br /> x = u.dimvalue(2)<br /> <br /> print('Plot...')<br /> figure(facecolor='k', newfig=False)<br /> ax = axes3d(projection=f.proj, clip_plane=False, orthographic=False, axis=False)<br /> volumeplot(x, y, height, speed, cmap='BlAqGrYeOrRe', ray_casting='specular',<br />     brightness=1.5, vmin=10, vmax=50,<br />     opacity_nodes=[10,20,30,35,40,50], opacity_levels=[0,0.05,0.1,0.2,0.8,1])<br /> colorbar(aspect=30, tickcolor='w')<br /> xlim(x[0], x[-1])<br /> ylim(y[0], y[-1])

ztqajw143

6啊，这也太炫酷了！

MgLily

哇！这看起来好酷

taiyang

哇哦，这个是WRF专门针对台风的模式嘛？画的好cool