added TorsionAngle function to Coordinates.f90 - to be compiled