Center of buoyancy function added and documented.

Check out the example plot comparison in ts.center.buoyancy. Compares thermoD to center.buoyancy. This is a metric that Dale Robertson and others worked with at CWR in the 90's.