在处理数据中，遇到多年按月或按季节求均值的问题，在这里做一记录！
1、加载nc数据

import xarray as xr
ds=xr.open_dataset(r'E:\Allmode\salem_russia\russia_pr\salemtrans_0.5_pr_Amon_BCC-CSM2-MR_ssp126_r1i1p1f1_gn_2020_2050.nc.nc').load()
pr=ds.pr*86400*30
lon=ds.lon
lat=ds.lat

打印一下：

2、截取研究时间段并按月分组

pr_year=pr.loc['2020-01-16':'2029-12-16']
data=pr_year.groupby('time.month')

3、对月数据进行提取，后面的1是代表时间维

mon_1=list(data)[0][1]
mon_2=list(data)[1][1]
mon_3=list(data)[2][1]
mon_4=list(data)[3][1]
mon_5=list(data)[4][1]
mon_6=list(data)[5][1]
mon_7=list(data)[6][1]
mon_8=list(data)[7][1]
mon_9=list(data)[8][1]
mon_10=list(data)[9][1]
mon_11=list(data)[10][1]
mon_12=list(data)[11][1]

4、求取DJA，春季的平均

DJA=(mon_12+mon_1+mon_2)/4#直接求平均，在这里的结果还是10年各自的春季
DJA_f10=DJA.mean(dim='time')#对上一步再操作，得到十年的春季平均