Replace Stumpf BJD conversion with sofa.
Fix `tsreduce bjd` epoch and increase output precision.