+ distributes across numpy arrays so you can write this as:

````
return np.unique(np.concatenate([base_maps + i * facet_offset for i in range(fs._mesh.layers - 1)]))
````
Note that unique flattens arrays anyway.