Added a sources argument for plot(), and make sure that the output_dir=None mode works for sources with no fits