diff --git a/Source/Background/KerrSchild.hpp b/Source/Background/KerrSchild.hpp index 514ce11..5fd30d3 100644 --- a/Source/Background/KerrSchild.hpp +++ b/Source/Background/KerrSchild.hpp @@ -244,7 +244,7 @@ class KerrSchild public: // used to decide when to excise - ie when within the horizon of the BH // note that this is not templated over data_t - bool check_if_excised(const Coordinates &coords) const + bool check_if_excised(const Coordinates &coords, const double innerH_buffer = 1.05, const double outerH_buffer = 0.9) const { // black hole params - mass M and spin a const double M = m_params.mass; @@ -270,7 +270,7 @@ class KerrSchild bool is_excised = false; // value less than 1 indicates we are within the horizon - if (outer_horizon < 0.9 || inner_horizon < 1.05) + if (outer_horizon < outerH_buffer || inner_horizon < innerH_buffer) { is_excised = true; }