JDE of perigee nearest the given decimal year.
year
double perigee(double year) { final l = _La(year, 0); return _mean(l.t) + l._pc(); }