On-the-fly decompression of downloaded fits.Z files

To be able to do that, I needed to wrap a subprocess call to gzip. I placed this wrapper in a new `system_tools.py` file in package `utils`. Is this the proper place? Does it sound like a good name?

With this pull request, once you have downloaded an archive file, you can keep working with `astropy.io.fits`...