Coordinates including ghost points

[boriskaus]

Currently, 1D coordinate arrays for the axis of the grid can be retrieved as follows:

julia> using Chmy, Chmy.Architectures, Chmy.Grids
julia> using KernelAbstractions

julia> grid = UniformGrid(Arch(CPU()); origin=(-1.0, 0.0), extent=(2.0, 2.5), dims=(10,8));

julia> coords(grid, Vertex())
(LinRange{Float64}(-1.0, 1.0, 11), LinRange{Float64}(0.0, 2.5, 9))

Yet this does not include the halo or ghost points.

It would be nice to have an option to return those as well (e.g., with_halo=true).
This would allow coupling with particle advection using JustPIC or would allow specifying initial model geometries in MPI
using GeophysicalModelGenerator