replaced pyfits with astropy in my code (but not in donutlib!)