使用xarray绘制多年来以x轴表示的双y轴的月平均值

，可以通过以下步骤完成：

首先，导入需要的库和模块：

import xarray as xr
import matplotlib.pyplot as plt

加载数据集：使用xarray的open_dataset方法打开数据集文件。

data = xr.open_dataset('data.nc')

数据处理：对数据进行筛选和处理，计算月平均值。

data_monthly = data.resample(time='1M').mean()

创建图表对象：使用matplotlib创建一个图表对象。

fig, ax1 = plt.subplots()

绘制第一个y轴的数据：使用ax1对象绘制第一个y轴的数据。

ax1.plot(data_monthly['x'], data_monthly['y1'], color='blue', label='y1')
ax1.set_xlabel('x')
ax1.set_ylabel('y1', color='blue')
ax1.tick_params('y', colors='blue')

绘制第二个y轴的数据：使用ax1对象的twinx()方法创建一个共享x轴的新y轴，并绘制第二个y轴的数据。

ax2 = ax1.twinx()
ax2.plot(data_monthly['x'], data_monthly['y2'], color='red', label='y2')
ax2.set_ylabel('y2', color='red')
ax2.tick_params('y', colors='red')

添加图例：使用ax1对象和ax2对象的legend()方法添加图例。

lines, labels = ax1.get_legend_handles_labels()
lines2, labels2 = ax2.get_legend_handles_labels()
ax1.legend(lines + lines2, labels + labels2, loc='upper right')