Use a more efficient implementation for mpdaf_sky_ref

Sort the lbda array once, and use an index array, instead of iterating on the whole pixtable for each wavelength bin. Faster than the previous implementation with only one core. (600s vs 625s on musicos1)