Description: Celestial-to-intermediate matrix given B-P-N.
The source code for the Fortran routine C2IBPN is available: