This MR adds the extra waveform arguments that are need to run with higher order modes. I haven't tested this on real data yet, so it's not ready to be merged. I do have a series of questions though, which I'll add as comments in the relevant diffs below.