Problems with TEMPO (Windows version)
The ephemeris data in the '0833-45.par' supplied with the Windows version of TEMPO is from epoch 51001 (July 1998). Looking at the table above it can be seen that after that epoch 6 glitches have occurred, 5 of which are > 2 ppm. Although the total of the 6 glitches amounts to nearly 13 ppm, the recovery mechanism means that the error will be less than that amount.
The latest ephemeris data for Vela from the ATNF database is from epoch 51559.32 (January 2000) and a small improvement in accuracy can be gained by updating the Windows TEMPO '0833-45.par' file to those values. The latest ephemeris data the author has is from the Hartebeesthoek Radio Astronomy Observatory (HartRAO) for epoch 57068 (February 2015) and a significant improvement in accuracy is obtained by using those values. It should be mentioned that obtaining the data from HartRAO was made possible by a sympathetic staff member of HartRAO who no longer works there.
The problem with Vela is that it will continue to glitch and so even the latest HartRAO ephemeris data will become out of date. To illustrate this, if the result from the Windows TEMPO epoch 51001 data is compared to the latest HartRAO epoch 57068 data result, there is a difference of 9.9 ppm in the prediction. So, over some 17 years there is almost a 10 ppm difference. However, until another glitch occurs, using TEMPO with the '0833-45.par' file updated to the latest HartRAO values is accurate. The next Vela glitch is due every three years and so the TEMPO prediction will be inaccurate again - and increasingly inaccurate over time unless a new ephemeris data set can be obtained.
To make available a result for Vela which doesn't drift in accuracy over time, the author loaded all the F0's found for Vela over a number of decades into a spreadsheet. Then a linear fit equation was derived. This equation includes all the glitches over time because it is derived from F0's over time. This provides a prediction curve which is smooth and continuous. Because of the sawtooth nature of the real pulse frequency curve over time (due to glitches), the smooth linear fit curve will over-estimate the pulse frequency in the time just before a glitch, but under-estimate the pulse frequency just after a glitch. In the time mid-way between glitches it will be very close to the real pulsar frequency. The error between the real pulsar frequency and the linear fit prediction is expected to remain less than ±2 ppm.
As a general rule when using the Windows version of TEMPO, a check should be made to ensure that latest epoch ephemeris data is loaded in the relevant '*.par' file. For example, for Vela the ATNF database data epoch is later than the epoch in the '0833-45.par' file supplied with the Windows version of TEMPO. Conversely, the Windows TEMPO version '0329+54.par' has a later epoch than the ATNF database - resulting in a 0.33 ppm difference in predictions.
It is worth checking for the latest epoch ephemeris...