use astropy units