diff --git a/uxarray/grid/coordinates.py b/uxarray/grid/coordinates.py
index 1c32e67a0..1113036e7 100644
--- a/uxarray/grid/coordinates.py
+++ b/uxarray/grid/coordinates.py
@@ -833,4 +833,4 @@ def points_atleast_2d_xyz(points):
             "Points are neither Cartesian (shape N x 3) nor Spherical (shape N x 2)."
         )
 
-    return points_xyz
+    return np.ascontiguousarray(points_xyz)
diff --git a/uxarray/grid/point_in_face.py b/uxarray/grid/point_in_face.py
index a622eb8dc..f83cf5ae0 100644
--- a/uxarray/grid/point_in_face.py
+++ b/uxarray/grid/point_in_face.py
@@ -205,6 +205,7 @@ def _point_in_face_query(
         Number of valid face‐hits per point.
     """
     pts = np.asarray(points, dtype=np.float64)
+    pts = np.ascontiguousarray(pts)
     if pts.ndim == 1:
         pts = pts[np.newaxis, :]
     # Cull with k-d tree
diff --git a/uxarray/plot/matplotlib.py b/uxarray/plot/matplotlib.py
index 5e43a30af..8ae686826 100644
--- a/uxarray/plot/matplotlib.py
+++ b/uxarray/plot/matplotlib.py
@@ -148,7 +148,10 @@ def _get_points_from_axis(ax: GeoAxes, *, pixel_ratio: float = 1):
     xs = np.cos(lats) * np.cos(lons)
     ys = np.cos(lats) * np.sin(lons)
     zs = np.sin(lats)
-    pts = np.vstack((xs, ys, zs)).T
+
+    # TODO:
+    pts = np.column_stack((xs, ys, zs))
+    # pts = np.vstack((xs, ys, zs)).T
 
     return pts, valid, nx, ny