Reshaping lost when accessing FITS_rec by field name

Please consider the code below and its output. The behaviour of `FITS_rec` does not match the one of `recarray`, which feels like the correct one.

When reshaping `FITS_rec`, of initial shape `(5,)`, the new shape `(5, 1)` is properly reported. However, when accessing a field, it behaves like the reshape never happened.

Maybe a proper fix would be that reshaping a `FITS_rec` should return a `recarray`.

``` python
from os import path
from astropy.io import fits
import numpy as np

with fits.open(path.join(fits.__path__[0], 'tests', 'data', 'stddata.fits')) as hdulist:

# As FITS_rec
data = hdulist[2].data
print(type(data))
print(data.shape)
print(data['PSFFLUX'].shape)
data = data.reshape((5,1))
print(data.shape)
print(data['PSFFLUX'].shape)

# As recarray
data = np.array(hdulist[2].data.copy())
print(type(data))
print(data.shape)
print(data['PSFFLUX'].shape)
data = data.reshape((5,1))
print(data.shape)
print(data['PSFFLUX'].shape)
```

``` python

(5,)
(5, 5)
(5, 1)
(5, 5)

(5,)
(5, 5)
(5, 1)
(5, 1, 5)
```

Start by running the supplied reproduction with astropy.io.fits against stddata.fits, comparing FITS_rec and the NumPy recarray behavior after reshape. Then inspect the FITS_rec reshape and field-access paths; done means field access preserves the reshaped dimensions consistently with recarray.

Written by the indexing model from the issue text.