Cross-section plotting with terrain as background

Question

In Atmospheric Research, the cross-section plots are often used to represent the average meteorology vertical profiles.

Using Python, I can plot the cross-section figure based on netcdf file(Simulation result).
I use this post for reference

But how to represent the terrain in specific cutting line. Example like this:

Answer 1

I've used

plt.fill_between(axis, terrain, 0, facecolor='black')