new_moon_atafter
double t = new_moon_atafter(midnight_in_china(rd));
t = new_moon_atafter(t);
while ((t_newmoon = new_moon_atafter(t_newmoon)) < t_end) {
double new_moon_atafter(double t);