PO2的长度本身就是和温度长度不相同的，也需要遍历PO2

非常抱歉，我理解您的意思了。下面是根据您的要求生成的代码，考虑了PO2和T变量的长度不相同的情况：

fig, ax1 = plt.subplots()

ax2 = ax1.twinx()

colors = ['r', 'g', 'b']

for i, temp in enumerate(T):
    var = sol.x
    var_2 = var[2][i]  # 获取对应温度的var[2]
    var_0 = var[0][i]  # 获取对应温度的var[0]
    PO2_temp = PO2[:min(len(PO2), len(var_2))]  # 确保PO2和var[2]的长度相同
    ax1.plot(PO2_temp, var_2, color=colors[i], label=f'Temperature {temp}')
    ax2.plot(PO2_temp, var_0, color=colors[i], linestyle='--')

ax1.set_xlabel('PO2')
ax1.set_ylabel('var[2]')
ax2.set_ylabel('var[0]')

ax1.legend(loc='upper left')
ax2.legend(loc='upper right')

plt.title('Plot with Dual Y-axis')
plt.show()

在这段代码中，我们使用了切片操作 `PO2[:min(len(PO2), len(var_2))]` 来确保 `PO2` 和 `var[2]` 的长度相同。这样可以避免在绘制图形时出现长度不匹配的问题。

请将此代码添加到您的程序中，并根据需要进行调整。运行代码后，将生成带有双坐标轴的图形，其中左侧纵坐标为 `var[2]`，右侧纵坐标为 `var[0]`，横坐标为 `PO2`。每个温度 `T` 对应一条曲线，曲线向左看对应左侧纵坐标，曲线向右看对应右侧纵坐标。同时，我们确保了 `PO2` 和 `var[2]` 的长度相同。

如果您有任何其他问题，请随时提问。