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