@peron
cool routine.
@guofu : note computation of Julian Date. and then some trig for the periodic part(s)