Create PA.py

Calculation of Parallactic Angles from (ra, dec) & (lat, long) + JD