metpy 1.0版如何使用vorticiy计算？

Question

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

新结构的*是什么？
当我使用(u, v)时，其他参数似乎是可选的，我面对ValueError: Must provide dx/dy arguments or input DataArray with latitude/longitude coordinates.
当我使用(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)

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

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

metpy 1.0版如何使用vorticiy计算？

How to use vorticiy calculation in metpy version 1.0?

python

metpy