There are also 2 ways to reduce eccentricity. The first way is described in EccRed_v1.md for use with the old BNSdata. The second way is described in EccRed_v2.md for DNSdata. Additional details about the 2 ways are in Ecc-Red-Notes.txt. For either way we need scripts to fit the orbital separation to the sinusoidal curve d(t)=S0+A0*t+0.5*A1*t^2-(B/omega)cos(omega*t+phi). For the 1st way we have EccRed.py (see e-based eccentricity reduction) and for the 2nd way we use EccRed_noForceBal.py (see Omega-based eccentricity reduction). All scripts can be found at https://github.com/sgridsource/EccRed.
To check that the fit was successful, one needs to plot the fit and the original data. This can be done with tgraph which is at https://github.com/wofti/tgraph.