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 in scripts.
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.