Change SOFA by ERFA

This PR should fix #1195. I have added erfa.c and erfa.h to cextern/erfa. (I have created them using erfa-fetch).

I have changed the references to sofa and the iau_prefix in astropy.time by erfa and era_. I have updated also the tests, the documentation, and the changelog