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

import xarray as xr
import matplotlib.pyplot as plt

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

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

fig, ax1 = plt.subplots()

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')

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')

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

plt.show()