NewRad will extrapolate into any boundary region where at least one face is marked as a physical output boundary.
This makes it impossible to use with Cartoon2D (and Cartoon2D’s tests fail if one tries) since Cartoon2D runs have only a single interior point that could be used as a source for extrapolation.
changes this so that NewRad does not interpolate into any region where at least one face is a symmetry boundary.
- symmetry boundaries are applied afterwards (likely since NewRad is triggerd by user code only in
CCTK_INITIALand not a boundary condition)
- the symmetry boundary conditions can fill in edges and corners where physical and symmetry boundaries overlap (true for the ones in the ET)