Calculating delta_frequency in source.py can cause inconsistent arrays
This was recently discovered. For certain sets of duration and sampling_frequency there are rounding errors that cause inconsistent shapes between frequency_array and the output modes.
Try running duration=32, sampling_frequency=4400 with an injected lalsim binary neutron star to reproduce this issue.
One solution might be replacing delta_frequency = frequency_array[1] - frequency_array[0]
with delta_frequency = np.mean(np.diff(frequency_array))