I managed to find the subsolar and sublunar longitude routines and they are the hard ones. Now I need to test the output for accuracy.

Roger