Adding Code and scripts to produce rad_vs_freq plots for a fixed wavenumber (in this case, 1).