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