You can provide either a position on the sky as a skyfield.starlib.Star object or an HFB source name as a string (which will be looked up in the HFB catalog).
You provide a list of unix times of length T and a list of rest frequencies of length F, and you get out an array of Doppler shifted frequencies of shape (F, T).
In case the source is in the HFB catalog, you can omit the rest frequency, in which case the rest frequency of the absorption feature as listed in the catalog will be used.
We could also split it up into a base function for Doppler shifting (to be used in a pipeline task) and a convenience function for quickly getting Doppler shifted frequencies of HFB sources with known absorption feature frequencies.
I'll work on the pipeline task for Doppler shifting and see how that interacts with this, but comments on this are welcome
This is currently set up as follows:
skyfield.starlib.Starobject or an HFB source name as a string (which will be looked up in the HFB catalog).Tand a list of rest frequencies of lengthF, and you get out an array of Doppler shifted frequencies of shape(F, T).We could also split it up into a base function for Doppler shifting (to be used in a pipeline task) and a convenience function for quickly getting Doppler shifted frequencies of HFB sources with known absorption feature frequencies.
I'll work on the pipeline task for Doppler shifting and see how that interacts with this, but comments on this are welcome