Transformation from J2000.0 equatorial coordinates to IAU 1958 Galactic coordinates. (double precision)
dr,dd double J2000.0 RA,Dec
*dl,*db double Galactic longitude and latitude l2,b2 (all arguments are radians)
slaDcs2c, slaDmxv, slaDcc2s, slaDranrm, slaDrange
The equatorial coordinates are J2000.0. Use the routine slaEg50 if conversion from B1950.0 'FK4' coordinates is required.
Blaauw et al, Mon.Not.R.astron.Soc.,121,123 (1960) P.T.Wallace Starlink 31 October 1993