MeteoInfoLab脚本示例：风场矢量图

2024-08-30 03:23:49

读取风场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()

最新文章

热门文章