html_url,issue_url,id,node_id,user,created_at,updated_at,author_association,body,reactions,performed_via_github_app,issue
https://github.com/pydata/xarray/issues/7042#issuecomment-1249353885,https://api.github.com/repos/pydata/xarray/issues/7042,1249353885,IC_kwDOAMm_X85Kd6Cd,4160723,2022-09-16T13:16:43Z,2022-09-16T13:16:43Z,MEMBER,You'd need to adjust the latitude / longitude values passed to `sel` so that your ROI contains (or is nearest to) more than one gird point in latitude or longitude. I'm afraid there's not automatic way in Xarray to do that.,"{""total_count"": 0, ""+1"": 0, ""-1"": 0, ""laugh"": 0, ""hooray"": 0, ""confused"": 0, ""heart"": 0, ""rocket"": 0, ""eyes"": 0}",,1374911473
https://github.com/pydata/xarray/issues/7042#issuecomment-1249295769,https://api.github.com/repos/pydata/xarray/issues/7042,1249295769,IC_kwDOAMm_X85Kdr2Z,4160723,2022-09-16T12:19:48Z,2022-09-16T12:19:48Z,MEMBER,"Hi @rogomichael,

It is expected in your example that `sel` with your given values and `method='nearest'` returns coordinates with duplicate longitude & latitude values.

The resolution of the grid is too coarse compared to your ROI, which is why the two provided values for latitude (longitude) both select the same grid point.

If you want to use `sel` for extracting a range over the latitude and longitude dimensions, it is better to use slices. It better illustrates that your ROI is smaller than the grid resolution:

```python
ds.sel(
    longitude=slice(49.294, 49.376),
    latitude=slice(10.971, 11.087)
)
# <xarray.Dataset>
# Dimensions:    (time: 365, longitude: 1, latitude: 0)
# Coordinates:
#   * time       (time) datetime64[ns] 2019-01-01 2019-01-02 ... 2019-12-31
#   * longitude  (longitude) float64 49.38
#   * latitude   (latitude) float64 
# Data variables:
#     tg         (time, latitude, longitude) float32 ...
``` ","{""total_count"": 0, ""+1"": 0, ""-1"": 0, ""laugh"": 0, ""hooray"": 0, ""confused"": 0, ""heart"": 0, ""rocket"": 0, ""eyes"": 0}",,1374911473
https://github.com/pydata/xarray/issues/7042#issuecomment-1248529133,https://api.github.com/repos/pydata/xarray/issues/7042,1248529133,IC_kwDOAMm_X85Kawrt,35968931,2022-09-15T19:36:46Z,2022-09-15T19:36:46Z,MEMBER,"Hi @rogomichael, thanks for raising your issue.

Could you possibly turn your problem into a [minimal reproducible example](https://github.com/pydata/xarray/discussions/5404)? That would make it much easier for us to help you.

You've included a lot of code, but from looking at it it seems to me that the only lines which matter are

```python
nc_file = xr.open_mfdataset('year2019.nc')

nc_file[""tg""].sel(
    longitude=[state_lon[0], state_lon[1]],
    latitude=[state_lat[0], state_lat[1]],
    method=""nearest""
)
```

This suggests to me that your original file has duplicate values, as all xarray has done here is open it.

(However I still am not going to be able to reproduce this error without access to the file `year2019.nc`.)

","{""total_count"": 0, ""+1"": 0, ""-1"": 0, ""laugh"": 0, ""hooray"": 0, ""confused"": 0, ""heart"": 0, ""rocket"": 0, ""eyes"": 0}",,1374911473