parts of 'suncalc'
Meeus, example 25.a, implemented
apply more tests