GRHD: Use WeakBoundaries() for CubicSplineInterpolation
Sometimes our grid bounds don't align exactly with the bounds from the data. This boundary treatment is more forgiving in that sense.
Sometimes our grid bounds don't align exactly with the bounds from the data. This boundary treatment is more forgiving in that sense.