lunar_phase_atafter
t = lunar_phase_atafter(180, t);
t_event = lunar_phase_atafter(phi_events[i], t_event);
double lunar_phase_atafter(double phi, double t);