Support masked arrays in broadcast_to_shape.

This PR adds support for masked arrays to `iris.util.broadcast_to_shape`. The striding trick needed to be applied independently to the mask of masked arrays, otherwise the result is a plain numpy array including the masked values.