Comment (by Roland Haas):

Thank you. I think with those things should be good to go. With the addition of adding this parfile as a test case to Carpet. Shiftout is ok since without the shiftout the first boundary point is exactly on the boundary (eg at -1.0) and there are not boundary_size points outside of the nominal domain (but boundary_size - 1). This is actually intentional since it helps with convergence to have the boundary point stay fixed.

Weren’t there plans at one point for SphericalNR to support / use mesh refinement on top of the spherical coordinates? In that case one would have to update this code it would seem.

--
Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2601/add-four-more-parameters-for-carpet