Optimize grid's `data_shape`

Grid property data_shape is called quite often, but involves considerable computation for certain grid types (e.g. calculating cell centers).

The data shape should be determined at construction, or be cached after the first call to data_shape.