Pass kwargs to sel (#304)

* Pass kwargs to sel

* Add test for sel with  method=nearest

* [pre-commit.ci] auto fixes from pre-commit.com hooks

for more information, see https://pre-commit.ci

* changelog

* stylistical changes

---------

Co-authored-by: pre-commit-ci[bot] <66853113+pre-commit-ci[bot]@users.noreply.github.com>
Co-authored-by: Justus Magin <[email protected]>

Author: 3 people

docs/whats-new.rst

@@ -4,7 +4,8 @@ What's new
 ==========
 0.5.1 (*unreleased*)
 --------------------
-
+- Pass ``sel`` options to the wrapped array (:pull:`304`, :issue:`303`)
+  By `Bhavin Patel <https://github.com/bpatel2107>`_.
 
 0.5 (09 Jun 2025)
 ------------------

pint_xarray/index.py

@@ -52,11 +52,11 @@ def stack(cls, variables, dim):
     def unstack(self):
         raise NotImplementedError()
 
-    def sel(self, labels):
+    def sel(self, labels, **options):
         converted_labels = conversion.convert_indexer_units(labels, self.units)
         stripped_labels = conversion.strip_indexer_units(converted_labels)
 
-        return self.index.sel(stripped_labels)
+        return self.index.sel(stripped_labels, **options)
 
     def isel(self, indexers):
         subset = self.index.isel(indexers)

pint_xarray/tests/test_index.py

@@ -104,6 +104,42 @@ def test_sel(labels, expected):
     )
 
 
+@pytest.mark.parametrize(
+    ["labels", "expected"],
+    (
+        (
+            {"x": ureg.Quantity(1.1, "m")},
+            IndexSelResult(dim_indexers={"x": np.array(0)}),
+        ),
+        (
+            {"x": ureg.Quantity(3100, "mm")},
+            IndexSelResult(dim_indexers={"x": np.array(2)}),
+        ),
+        (
+            {"x": ureg.Quantity(0.0021, "km")},
+            IndexSelResult(dim_indexers={"x": np.array(1)}),
+        ),
+        (
+            {"x": ureg.Quantity([0.0021, 0.0041], "km")},
+            IndexSelResult(dim_indexers={"x": np.array([1, 3])}),
+        ),
+    ),
+)
+def test_sel_nearest(labels, expected):
+    index = PintIndex(
+        index=PandasIndex(pd.Index([1, 2, 3, 4]), dim="x"), units={"x": ureg.Unit("m")}
+    )
+
+    actual = index.sel(labels, method="nearest")
+
+    assert isinstance(actual, IndexSelResult)
+    assert actual.dim_indexers.keys() == expected.dim_indexers.keys()
+    assert all(
+        indexer_equal(actual.dim_indexers[k], expected.dim_indexers[k])
+        for k in expected.dim_indexers.keys()
+    )
+
+
 @pytest.mark.parametrize(
     "indexers",
     ({"y": 0}, {"y": [1, 2]}, {"y": slice(0, None, 2)}, {"y": xr.Variable("y", [1])}),