The hard-coded number 2.3504 * 10**(-5) should be compute-able with the units package: it's just the conversion between MJy/steradians and Jy/arcsec.
I was wonderng if it would be better to store some of the other numbers in some kind of 'calibration database' so that the code doesn't need to be updated -- but I think astropy discourages this?
The hard-coded number 2.3504 * 10**(-5) should be compute-able with the units package: it's just the conversion between MJy/steradians and Jy/arcsec.
I was wonderng if it would be better to store some of the other numbers in some kind of 'calibration database' so that the code doesn't need to be updated -- but I think astropy discourages this?