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