I see.  There is a bug here, that the FITS parser is not using the correct version of `_parse_unit`.  I think that probably broke when we moved from pyparsing to PLY.  I'll try to come up with a fix and test for this, probably against this PR, just because otherwise the merging of the two could get hairy.