如何在metpy 1.0版中使用涡度计算

Question

我想使用metpy.calc.vorticity。以前，它的结构是(u, v, dx, dy)，但改为(u, v, *, dx=None, dy=None, x_dim=- 1, y_dim=- 2)。

1. 新结构的*是什么
2. 当我使用(u, v)时，似乎其他参数是可选的，我面对ValueError: Must provide dx/dy arguments or input DataArray with latitude/longitude coordinates.
3. 当我使用(u , v, dx, dy)时，我面对的是TypeError: too many positional arguments

然而，我的代码和数据中的geostrophic_wind(height,dx,dy,latitude)没有问题，但涡度的问题是什么？我的代码的某些部分如下：

PATH_nc = '/home/python_script/'
out_fff='gfs.t12z.pgrb2.0p50.f168'
ncfile = Dataset('{}{}.nc'.format(PATH_nc,out_fff))
lon = ncfile.variables['longitude'][:]
lon[lon < 0] = 360 + lon[lon < 0]
lat = ncfile.variables['latitude'][:]
dx, dy = mpcalc.lat_lon_grid_deltas(lon, lat)
long1, lat1 = np.meshgrid(lon, lat)
hght_850  = ncfile.variables['HGT_850mb'][:]
HGT_850  = gaussian_filter((hght_850)/10, sigma=  2)
height850 = ndimage.gaussian_filter(hght_850, sigma=1.5, order=0)
geo_wind_u_850, geo_wind_v_850 = mpcalc.geostrophic_wind(np.squeeze(height850)*units.m, dx, dy, lat1)
vwnd_700 =  units('m/s') *ncfile.variables['VGRD_700mb'][:]
avor700 = mpcalc.vorticity(np.squeeze(uwnd_700), np.squeeze(vwnd_700), dx, dy)

我希望在我的代码中写下与问题相关的所有必要部分。我将感谢任何帮助。

Answer 1

所以在vorticity:的定义中

metpy.calc.vorticity(u, v, *, dx=None, dy=None, x_dim=-1, y_dim=-2)

*之后的所有内容都被定义为仅关键字参数，这意味着函数参数/参数只能通过关键字传递。因为您当前的代码不是通过关键字传递dx和dy，而是仅通过位置传递，所以您得到的错误是：TypeError: too many positional arguments。修复方法是使用关键字参数传递dx和dy，如：

avor700 = mpcalc.vorticity(np.squeeze(uwnd_700), np.squeeze(vwnd_700), dx=dx, dy=dy)