关于python绘制风场图

如何选择特定气压值下的风场图
import xarray as xr
import numpy as np
import datetime as dt
import cartopy.crs as ccrs
import cartopy.feature as cfeature
import matplotlib.pyplot as plt
from cartopy.mpl.gridliner import LONGITUDE_FORMATTER, LATITUDE_FORMATTER
import matplotlib.ticker as mticker
filename1 =r'F:/python/uwnd.2021.nc'
filename2 =r'F:/python/vwnd.2021.nc'

ds1=xr.open_dataset(filename1)
ds2=xr.open_dataset(filename2)

lon=ds1['lon']
lat=ds1['lat']

lons,lats=np.meshgrid(lon,lat)
u=ds1['uwnd'][124[2][][]
v=ds2['vwnd'][124][2][][]

proj = ccrs.PlateCarree(central_longitude=0)
plt.rcParams['font.sans-serif'] = ['SimHei']
plt.rcParams['axes.unicode_minus'] = False
fig = plt.figure(figsize=(10,9),dpi=500) # 创建画布
#ax=fig.subplots(1,1,subplot_kw={'projection':proj}) # 创建子图
ax=plt.axes([0,0,1,1],projection=ccrs.PlateCarree(central_longitude=0))

ax.barbs(lons[:2,::2],lats[:2,::2],u[:2,::2],v[:2,:2],
barb_increments={'half':2,'full':4,'flag':20},
length=5,zorder=5)
plt.title("2021年2月中国区域风场图",size=20)
region=[72,136,0,55]#投影的经纬度
ax.set_extent([72,136,0,55],crs=proj)

ax.coastlines('50m',lw=0.5)

extent=[72,136,0,55]
gl=ax.gridlines(crs=ccrs.PlateCarree(),draw_labels=True,
linewidth=0.2,color='k',alpha=0.5,linestyle='--')
gl.xlabels_top=False
gl.ylabels_right=False
gl.xformatter=LONGITUDE_FORMATTER
gl.yformatter=LATITUDE_FORMATTER
gl.xlocator=mticker.FixedLocator(np.arange(extent[0],extent[1],15))
gl.ylocator=mticker.FixedLocator(np.arange(extent[2],extent[3],10))
plt.show（）

u=ds1
v=ds2
这串代码中124和2分别代表什么？

怎么能够选择想要的气压值和时间下的风场图

grads 也无法读取nc文件