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