Modernize: use astropy and add docstrings