spatial_gp(scale_coords = TRUE) (and the HSGP and multi-scale
constructors) centre and scale the coordinates before the field is built,
so the fitted field lives in those units. gp_scale_coords() scales the
training coordinates and returns the centre and scale it used;
gp_scale_new_coords() puts prediction coordinates, which arrive raw, in
the same units. A fit made without scaling leaves them unchanged.