计算涡度脚本报错求助！

无尽的兰博

用王老师那个计算涡度的脚本画了图，用的是grib2的数据。我想要200hPa那层的图，编程如下
f = addfile('F:/data/study/fnl/fnl_20160622_18_00.grib2')
u = f['u-component_of_wind_pressure_difference_layer'][0,200,:,:]
v = f['v-component_of_wind_pressure_difference_layer'][0,200,:,:]
vort = hcurl(u, v)
axesm()
mlayer = shaperead('F:/MeteoInfo/map/country1.shp')
geoshow(mlayer)
layer = contourfm(vort, 20)
title('Vorticity')
colorbar(layer, orientation='horizontal', aspect=50)


报错如下：

Index out of range!
Index out of range!
Traceback (most recent call last):
  File "<iostream>", line 8, in <module>
  File "F:\MeteoInfo\pylib\mipylib\miplot.py", line 3984, in contourfm
    ls = __getlegendscheme(args, gdata.min(), gdata.max(), **kwargs)
AttributeError: 'NoneType' object has no attribute 'min'



不知道是哪里出了问题，需要怎么修改呢？谢谢各位了！

无尽的兰博

求助

无尽的兰博

有没有哪位来帮我一下~~

MeteoInfo

无尽的兰博 发表于 2016-10-12 16:08
有没有哪位来帮我一下~~

已经在另一个帖子里回复你了

wuyula

haonanaaaaaaaaaaaaaa