NavList:

Hello Paul,

Thank you very much for taking the time to check my values.

I also checked my position angle function, which i got from sofa library. it is like this:

**  Given:
**     al     double     longitude of point A (e.g. RA) in radians
**     ap     double     latitude of point A (e.g. Dec) in radians
**     bl     double     longitude of point B
**     bp     double     latitude of point B

{
   double dl, x, y, pa;


   dl = bl - al;
   y = sin(dl) * cos(bp);
   x = sin(bp) * cos(ap) - cos(bp) * sin(ap) * cos(dl);
   pa = ((x != 0.0) || (y != 0.0)) ? atan2(y, x) : 0.0;

   return pa;

I am happy that you agree on other values.

There is always room to learn something new.

Have a great week.

Mehmet