Added function VGP_calc for calculating poles from directional data