diff --git a/include/bout/mask.hxx b/include/bout/mask.hxx index fc5f814676..9033080d44 100644 --- a/include/bout/mask.hxx +++ b/include/bout/mask.hxx @@ -26,6 +26,7 @@ #include "bout/globals.hxx" #include "bout/mesh.hxx" +#include "bout/region.hxx" /** * 3D array of bools to mask Field3Ds @@ -80,4 +81,12 @@ inline std::unique_ptr> regionFromMask(const BoutMask& mask, } return std::make_unique>(indices); } + +inline BoutMask maskFromRegion(const Region& region, const Mesh* mesh) { + BoutMask mask{mesh, false}; + + BOUT_FOR(i, region) { mask[i] = true; } + return mask; +} + #endif //BOUT_MASK_H