Sinusoidal phase functions R_j(t) (red lines) that fit the JBE phase residuals (blue points) for episodes from j = 1 to 10 (the short episode 5 is included in the graph showing episode 6). Time on the abscissa is in MJD, and the scale is the same for all episodes. Green dots are numbered as in Espinoza et al. (2011) at dates of published glitches with the ordinate as calculated from R_j(t). Ordinate scales are different and adjusted to different amplitudes of oscillations. Each panel shows data of the complete episode and includes the first point of the next episode. The difference between horizontal dotted lines represents a measure of the strength of the perturbation R_j(t) caused by the glitch. Dashed gray curves are the plots of ϕ_j(t)−ϕ_j−1(t), to show the difference of phases between two contiguous episodes. Gray arrows point to dates of the > 100 MeV X-ray flares detected by Fermi (Abdo et al. 2011; Ojiha et al. 2013; Bühler et al. 2012) and AGILE (Tavani et al. 2011; Striani et al. 2013).