The surface atoms of the slab can be determined by indexing the non-zero components of either of the two off-diagonal quadrants, since these only include the out-of plane connectivity. Indexed atoms in the top-half of the unit cell can then be robustly identified as top or bottom terminated; a pivotal first step for correct identification of active sites.