Description: Parallactic angle for a given hour angle and declination.
The source code for the Fortran routine HD2PA is available: