MeteoInfoLab脚本示例：风场矢量图

MeteoInfo

读取风场U/V变量数据，可以从U/V计算出风速：speed = sqrt(u*u+v*v)。quiverm函数用来绘制风场矢量图，参数中包括U/V变量，如果要绘制彩色风场还需要第三个变量，这里是风速speed，10是颜色的分级数，size=8指定了风场箭头的基准大小。

脚本程序：

1. f = addfile('D:/Temp/GrADS/model.ctl')
   
2. u = f['U'][0,'500','10:60','60:140']
   
3. v = f['V'][0,'500','10:60','60:140']
   
4. speed = sqrt(u*u+v*v)
   
5. #Plot
   
6. axesm()
   
7. geoshow('country')
   
8. layer = quiverm(u, v, speed, 10, size=8)
   
9. quiverkey(layer, 0.74, 0.18, 15, bbox={'edge':True, 'fill':True})
   
10. colorbar(layer)
    
11. title('Wind field')
    
12. yticks([20,40,60])

f = addfile('D:/Temp/GrADS/model.ctl')<br /> u = f['U'][0,'500','10:60','60:140']<br /> v = f['V'][0,'500','10:60','60:140']<br /> speed = sqrt(u*u+v*v)<br /> #Plot<br /> axesm()<br /> geoshow('country')<br /> layer = quiverm(u, v, speed, 10, size=8)<br /> quiverkey(layer, 0.74, 0.18, 15, bbox={'edge':True, 'fill':True})<br /> colorbar(layer)<br /> title('Wind field')<br /> yticks([20,40,60])

Linda.W

王老师，您好，我看到nc数据画风场矢量图，但是都是固定一个时次的，我想画多个时次求平均的季平均图，不知道该如何实现，希望王老师能够指点，谢谢老师。附上我试着编写的程序，一直有错误。@MeteoInfo
fs=addfile(r'I:\Learning\data2018\ECMWF\ws2013.nc')
uwnd=fs['u10']
vwnd=fs['v10']
us=uwnd['1:360','25:36','96:110']
vs=vwnd['1:360','25:36','96:110']
u=us.ave()
v=vs.ave()
speed=sqrt(u*u+v*v)
axesm()
layer = quiverm(u, v, speed, 10, size=8)
colorbar(layer)

MeteoInfo

unique 发表于 2020-7-9 15:55
王老师，用该脚本画风场图，与Grads的结果不一致，见附件。找了半天没找到原因，请王老师指教。
f = addfi ...

layer = quiverm(u, v, size=200, width=0.5)

piwa

U/V变量数据在哪下载呢？我想做青藏高原的风场矢量图，求大神指点啊

ldepn

菜鸟一个，弱弱地问一句什么是风场矢量图？从风场矢量图能得到一些什么信息？求大神普及一下。问度娘没有找到答案！

oubiwang

楼主，为什么我完全按您的脚本运行会报错？好像是说不支持相乘？

MeteoInfo

oubiwang 发表于 2015-12-15 18:59
楼主，为什么我完全按您的脚本运行会报错？好像是说不支持相乘？

脚本不会错，你看看MeteoInfo是不是最新版本，或者数据文件是不是有问题。

oubiwang

MeteoInfo 发表于 2015-12-15 19:17
脚本不会错，你看看MeteoInfo是不是最新版本，或者数据文件是不是有问题。

谢谢楼主，应该是版本问题

3870383

找到了，在设置里，两个锤子图标，能控制箭头密度

1铁芒萁1

NCEP/NCAR里下载

zenghaode

请问，我在程序目录下并没有U，V数据文件，怎么也能绘出图呢，是从什么地方读取的呢，谢谢！！

半颗麦芽糖

王老师，看了您的这个帖子，再结合该帖http://bbs.06climate.com/forum.php?mod=viewthread&tid=45172&extra=page%3D1
我想画一个wrfout风场矢量图，但是现在存在一个问题，wrfout里面垂直分层是采用的eta坐标，我不知道每个值与气压是如何对应的，然后我想画500hPa的图，请问脚本上能否实现这个？还有我在画图的时候发现图像上太密了，怎么调都没办法，请问王老师该怎么处理？谢谢