Description: Celestial-to-terrestrial matrix given CIP.
The source code for the Fortran routine C2TXY is available: