xy2sky implemented and tested against Fortran output.  Values are the 
same up to a delta of 0.000005
Removed get_order function as nord is a parameter of xy2sky