This seems to be a very basic question, but after reading the help function and searching on the internet, I still cannot find a solution. Excuse me if I missed something obvious here.

This seems to be a very basic question, but after reading the help function and searching on the internet, I still cannot find a solution. Excuse me if I missed something obvious here.

This seems to be a very basic question, but after reading the help function and searching on the internet, I still cannot find a solution. Excuse me if I missed something obvious here.

This seems to be a very basic question, but after reading the help function and searching on the internet, I still cannot find a solution. Excuse me if I missed something obvious here.

This seems to be a very basic question, but after reading the help function and searching on the internet, I still cannot find a solution. Excuse me if I missed something obvious here.

This seems to be a very basic question, but after reading the help function and searching on the internet, I still cannot find a solution. Excuse me if I missed something obvious here.

I am trying to plot energy surfaces of triatomic molecules, A-B-C, in Python/Matplotlib. I have tabulated data of bond length (B1 and B2) combinations and computed molecular energies (MP3s) which I want to convert to a smooth 3D surface plot. A similar Excel Macro function output the following:
Excel surface output

I want to visualize a catenoid in a 3D plot with matplotlib, but instead of coloring it by height, which is the default behaviour for the built in colormaps, I want to color it by gaussian curvature – similar to what was asked in this, but in python instead.