JDE of New Moon nearest the given decimal year.
year
double newMoon(double year) { final m = _Mp(year, 0); return _mean(m.t) + m._nfc(_nc) + m._a(); }