tgc function

Calculate contour geodesic torsion