There are 2 classes of physical of boundary conditions - those that store the value ON the boundary in the first ghost cell, and those that fill ALL the ghost cells with values extrapolated to those cell-centers.
The latter case should "Just work", e.g., symmetry and hoextraptocc, so I removed the Abort in those cases.
For the former case I previously was careful to make sure that interpolated values to the actual boundary in the BDS algorithm respected the boundary value that was stored in the ghost cell.
Address Issue #96
There are 2 classes of physical of boundary conditions - those that store the value ON the boundary in the first ghost cell, and those that fill ALL the ghost cells with values extrapolated to those cell-centers.
The latter case should "Just work", e.g., symmetry and hoextraptocc, so I removed the Abort in those cases.
For the former case I previously was careful to make sure that interpolated values to the actual boundary in the BDS algorithm respected the boundary value that was stored in the ghost cell.