提取grib2数据中各气压层的风速有问题，如何修改？

guo066600

MI读取grib2数据，代码如下：
infn ='D:/Temp/fnl_20190401_00_00.grib2'
f = addfile(infn)

文件信息中变量的格式为（time,pressure_difference_layer1,lat,lon)
想读取某经纬度范围内850hPa这一层对应的风速u和v，并计算speed,代码如下
u = f['u-component_of_wind_pressure_difference_layer'][0,'85000','15:50','100:140']
v = f['v-component_of_wind_pressure_difference_layer'][0,'85000','15:50','100:140']
speed= sqrt(u*u+v*v)

虽然没有报错，但从画的风场图看，这样表示不对的，不知王老师和各位高手能不能帮忙改正一下？