如何在Python中一次绘制速度剖面和边界厚度剖面？

import matplotlib.pyplot as plt
import numpy as np

depth = np.arange(0, 1000, 10)  # 深度范围，从0到1000，步长为10
velocity = np.random.uniform(1000, 5000, len(depth))  # 随机生成速度数据
thickness = np.random.uniform(0, 100, len(depth))  # 随机生成厚度数据

fig, ax1 = plt.subplots()  # 创建一个绘图窗口和一个坐标轴对象

ax1.plot(velocity, depth, color='blue', label='Velocity')  # 绘制速度剖面曲线
ax1.set_xlabel('Velocity (m/s)')  # 设置X轴标签
ax1.set_ylabel('Depth (m)')  # 设置Y轴标签
ax1.set_ylim(max(depth), min(depth))  # 设置Y轴刻度范围（从最大深度到最小深度）

ax2 = ax1.twiny()  # 创建第二个坐标轴，共享Y轴
ax2.plot(thickness, depth, color='red', label='Thickness')  # 绘制边界厚度剖面曲线
ax2.set_xlabel('Thickness (m)')  # 设置第二个坐标轴的X轴标签
ax2.set_xlim(max(thickness), min(thickness))  # 设置第二个坐标轴的X轴刻度范围（从最大厚度到最小厚度）

ax1.legend(loc='upper right')  # 添加第一个坐标轴的图例
ax2.legend(loc='upper left')  # 添加第二个坐标轴的图例
plt.title('Velocity and Thickness Profiles')  # 设置图表标题

plt.show()  # 显示图表