yaqiang

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()

 

分类:

技术点:

相关文章: