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