Introduce descriptor variables to Grid that stores the indices of faces that contain pole points:
Grid.north_pole_face_indices
Grid.north_pole_face_indices
Can use uxarray.grid.geometry._pole_point_inside_polygon() to classify the faces, which is already implemented and has been rigorously tested.