JDE of First Quarter nearest the given decimal year.
year
double first(double year) { final m = _Mp(year, 0.25); return _mean(m.t) + m._flc() + m._w() + m._a(); }