Added units to WCS coordinate column

WCS coordinate units are determined automatically within astropy also in case of absent CUNIT. This may fail with non-standard FITS headers but has the significant advantage that we do not need to set the units manually going through all possible allowed CTYPEs.