import numpy as np
import xarray as xr
import matplotlib.pyplot as plt
import cartopy.crs as ccrs
import cartopy.feature as cfeature
import cartopy.mpl.ticker as cticker
import datetime as dt
from cartopy.util import add_cyclic_point
from cartopy.mpl.ticker import LongitudeFormatter, LatitudeFormatter
q =r'C:\Users\86182\Desktop\lunwenshuju\q.nc'
u = r'C:\Users\86182\Desktop\lunwenshuju\u.nc'
v = r'C:\Users\86182\Desktop\lunwenshuju\v.nc'
du = xr.open_dataset(u)
dv = xr.open_dataset(v)
dq = xr.open_dataset(q)
u = du['u']
v = dv['v']
q = dq['q']
u = u.loc['2001-01-01':'2020-12-01',:,:,:]
v = v.loc['2001-01-01':'2020-12-01',:,:,:]
q = q.loc['2001-01-01':'2020-12-01',:,:,:]
q = q.loc[q.time.dt.month.isin(7)]
u
11-23
399
![](https://csdnimg.cn/release/blogv2/dist/pc/img/readCountWhite.png)
04-16
3万+
![](https://csdnimg.cn/release/blogv2/dist/pc/img/readCountWhite.png)
02-03
3062
![](https://csdnimg.cn/release/blogv2/dist/pc/img/readCountWhite.png)
“相关推荐”对你有帮助么?
-
非常没帮助
-
没帮助
-
一般
-
有帮助
-
非常有帮助
提交