diff --git a/benchmarks/mpas_ocean.py b/benchmarks/mpas_ocean.py
index 1db46e0f9..659c19014 100644
--- a/benchmarks/mpas_ocean.py
+++ b/benchmarks/mpas_ocean.py
@@ -208,13 +208,17 @@ def setup(self, resolution):
         point = np.array([0.0, 0.0, 1.0])
         res = self.uxgrid.get_faces_containing_point(point)
 
+        self.point_xyz = np.array([self.uxgrid.face_x[0].values, self.uxgrid.face_y[0].values, self.uxgrid.face_z[0].values], dtype=np.float64)
+        self.point_lonlat = np.array([self.uxgrid.face_lon[0].values, self.uxgrid.face_lat.values[0]], dtype=np.float64)
+
     def teardown(self, resolution):
         del self.uxgrid
 
-    def time_face_search(self, resolution):
-        point_xyz = np.array([self.uxgrid.face_x[0].values, self.uxgrid.face_y[0].values, self.uxgrid.face_z[0].values], dtype=np.float64)
-        point_lonlat = np.array([self.uxgrid.face_lon[0].values, self.uxgrid.face_lat.values[0]], dtype=np.float64)
-        self.uxgrid.get_faces_containing_point(point_xyz=point_xyz, point_lonlat=point_lonlat)
+    def time_face_search_xyz(self, resolution):
+        self.uxgrid.get_faces_containing_point(self.point_xyz)
+
+    def time_face_search_lonlat(self, resolution):
+        self.uxgrid.get_faces_containing_point(self.point_lonlat)
 
 
 class ZonalAverage(DatasetBenchmark):