read psrfits using pfits