fix: specify interpolation method in multidimensional interpolation test

Author: QianyeSu

tests/test_interp_interpolation.py

@@ -1183,7 +1183,9 @@ def test_interp_multidim_edge_case_coordinates(self):
         lat_out = np.array([45])
         lon_out = np.array([90])
 
-        result = interp_multidim(data_in=data_in, lat_out=lat_out, lon_out=lon_out)
+        result = interp_multidim(
+            data_in=data_in, lat_out=lat_out, lon_out=lon_out, method="nearest"
+        )
 
         assert result.shape == (1, 1)
         assert result.values[0, 0] == 5