-Add 21cm frequency
-Use astropy unit for unit of frequency