Moon rise and set calculations

Using quadratic interpolation with 24 control points. Maybe not super-accurate but should be enough for many cases. Ref #7