MeteoInfoLab脚本示例：风场矢量图

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



脚本程序：

f = addfile('D:/Temp/GrADS/model.ctl')
u = f['U'][0,[500],[10,60],[60,140]]
v = f['V'][0,[500],[10,60],[60,140]]
speed = sqrt(u*u+v*v)
axesm()
mlayer = shaperead('D:/Temp/map/country1.shp')
geoshow(mlayer)
layer = quiverm(u, v, speed, 10, size=8)
title('Wind field')
colorbar(layer)
show()