MeteoInfoLab脚本示例：风场矢量图

rceclx

半颗麦芽糖 发表于 2016-11-30 16:23
王老师，我利用nc格式资料通过脚本和界面画了图，但是不知道为什么画出的图有差异（图中气旋所在），画 ...

楼主，u = f['U'][103,[500],[20,50],[70,120]]        
高度、纬度和经度用[]表示提示越界错误，改用''没问题，是什么原因啊？

rceclx

王老师，u = f['U'][0,[500],[10,60],[60,140]]         
高度、纬度和经度用[]表示提示越界错误，改用' ' 没问题，是什么原因啊？

rceclx

王老师，u = f['U'][0,[500],[10,60],[60,140]]         
高度、纬度和经度用[]表示提示越界错误，改用' ' 没问题，是什么原因啊？

unique

王老师，用该脚本画风场图，与Grads的结果不一致，见附件。找了半天没找到原因，请王老师指教。
f = addfile('G:/MeteoInfo_2.2.6/nc/UVdif.nc', 'r')
u = f['ufinal'][:,:]
u = u[::12,::12]
v = f['vfinal'][:,:]
v= v[::12,::12]
speed = sqrt(u*u+v*v)
#Plot
axesm()
geoshow('country')
layer = quiverm(u, v, speed, 1, size=200, color=1,scale=0,width=0.5)
quiverkey(layer, 0.74, 0.18, 0.5, bbox={'edge':True, 'fill':True})
yticks([20,40,60])

MeteoInfo

unique 发表于 2020-7-9 15:55
王老师，用该脚本画风场图，与Grads的结果不一致，见附件。找了半天没找到原因，请王老师指教。
f = addfi ...

layer = quiverm(u, v, size=200, width=0.5)

guo066600

请问王老师，如何将图中的文字(标题/经纬度和图标)改为Times New Roman?

银河骑士陈先生

王老师，这个脚本可以读取分析GDAS1数据嘛？

MeteoInfo

银河骑士陈先生 发表于 2020-9-1 22:27
王老师，这个脚本可以读取分析GDAS1数据嘛？

可以

银河骑士陈先生

MeteoInfo 发表于 2020-9-2 08:44
可以

谢谢王老师！

guo066600

王老师，这里用来画风场图的数据格式没有要求吗？grib2,nc格式的都可以，对吗？