计算平均值超过更多变量的问题

Question

我有一个很大的NetCDF文件，该文件每月的平均降雨量在1987-2016年期间为世界范围。我必须分别计算每个年度和每个JJA季节（每年）的平均值。我该怎么办？我尝试过

europe = europe['precip'].groupby('time.year','time.season').mean('time', keep_attrs= True)

但是它不起作用，如果我分别按年（第一）和季节（第二）分组，则季节不再是时间序列的属性。

任何建议将不胜感激。