This line should be:
omega_nyq_line = np.array((np.nan, nyquistfrq_plot, nyquistfrq_plot))

... so that the nyquist limit will be plotted in Hz when requested. Currently the nyquist limit always plots as rad/s, even if the Hz argument is True.