Added routine to calculate parallactic angle