Meshblock identification for domain edge

[kyriakosTapinou]

I'm trying to identify which meshblocks are on the domain edge (specifically the low-z-boundary). I am trying to figure out how the boundary condition marking (BoundaryFlag etc.) work. If someone could point me towards the exact meshblock metadata that can be accessed publicly then that would do.

Previously I was using the LogicalLocation data structure e.g. loc.lx3==0 but I think this is wrong when refinement is used.

[kyriakosTapinou]

I figured out the issue and closed the discussion.

[tomidakn]

Just in case. You can use loc.lx3==0, but if you want to identify the right hand side, compare loc.lx3 with (pmesh->nrbx3 << (loc.level - pmesh->root_level) - 1). Alternatively, if (pmb->pbval->nblevel[0][1][1] == -1) means it is touching a physical boundary.