Support non-trivial horizontal coordinates for single column grids

[glwagner]

This PR makes it so we can build single column grids with non-trivial horizontal coordinates and nice syntax.

Resolves #3626

@simone-silvestri

[glwagner]

The kind of thing we want to do:

from_grid = RectilinearGrid(size=(4, 4, 4), x=(1, 2), y=(2, 3), z=(0.5, 2.5), topology=(Periodic, Periodic, Bounded))
xyz = CenterField(from_grid)
set!(xyz, (x, y, z) -> x + y + z)
Oceananigans.BoundaryConditions.fill_halo_regions!(xyz)

column_grid = RectilinearGrid(size=2, x=1.3, y=2.7,  z=(0.5, 2.5), topology=(Flat, Flat, Bounded))
xyz_column = CenterField(column_grid)
Oceananigans.Fields.interpolate!(xyz_column, xyz)

which gives

julia> xyz_column[1, 1, 1:2]
2-element Vector{Float64}:
 5.0
 6.0

[glwagner]

@simone-silvestri @navidcy please review

[glwagner]

@navidcy any idea what the documenter error is?