Deep copy grid using xarray.DataArray.copy instead of numpy.copy

After flipping the grid, try doing the deep copy on xarray (which should return a c_contiguous array) before retrieving the data (i.e. the numpy array).

Author: weiji14

pygmt/clib/conversion.py

@@ -113,7 +113,7 @@ def dataarray_to_matrix(grid):
         inc = [abs(i) for i in inc]
         grid = grid.sortby(variables=list(grid.dims), ascending=True)
 
-    matrix = as_c_contiguous(grid[::-1].values)
+    matrix = as_c_contiguous(grid[::-1].copy(deep=True).data)
     return matrix, region, inc