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