twinx non-linear mapping between shared y axes on a plot in matplotlib

Question

I have some 3d data that I am plotting with pcolormesh.

On the x-axis is time, on the y-axis is height.

The height has a potential (E) associated with it, but