Grid masks

We need masked grids of some kind. Masks grids don't work with xarray anymore. Further, it was suggested to add masks to the specification rather than to the data.

Back to NumPy, we have numpy.ma for this. We should check how this integrates with pint. (Seems to be the case according to a first trial).

There is still the question about integrating it into grid specs.