Description: Parallactic angle for a given hour angle and declination.
The source code for the ANSI C routine Hd2pa is available: