If you print the reference frequency of an FHD catalog read in with pyradiosky, the units are in Hz when they should be in MHz. I think one way to fix this would be to change reference_frequency=Quantity(source_freqs, "hertz") in read_fhd_catalog such that source_freqs is multiplied by 10**6 (ideally with an if statement checking the magnitude first).
If you print the reference frequency of an FHD catalog read in with pyradiosky, the units are in Hz when they should be in MHz. I think one way to fix this would be to change
reference_frequency=Quantity(source_freqs, "hertz")inread_fhd_catalogsuch thatsource_freqsis multiplied by10**6(ideally with an if statement checking the magnitude first).