def chi_plot (table, sid, band = 'j', harmonic=1) :
''' Does everything for one source and plots a phase plot.'''
f = chi_analyze (table, sid, band)
period = 1./f
tr_plot.plot_phase (table, sid, period*harmonic, band)
return period
