When running the example I am getting an error when printing the data:

print(data.Geopotential_height.coords)

hght_500 = data['Geopotential_height_isobaric'].sel(time1=vtimes[0], isobaric=500)
uwnd_500 = data['u-component_of_wind_isobaric'].sel(time1=vtimes[0], isobaric=500)
vwnd_500 = data['v-component_of_wind_isobaric'].sel(time1=vtimes[0], isobaric=500)

I get below error:

KeyError Traceback (most recent call last)
Cell In[20], line 1
----> 1 hght_500 = data['Geopotential_height_isobaric'].sel(time1=vtimes[0], isobaric=500)
2 uwnd_500 = data['u-component_of_wind_isobaric'].sel(time1=vtimes[0], isobaric=500)
3 vwnd_500 = data['v-component_of_wind_isobaric'].sel(time1=vtimes[0], isobaric=500)

File ~/anaconda3/lib/python3.13/site-packages/xarray/core/dataarray.py:1717, in DataArray.sel(self, indexers, method, tolerance, drop, **indexers_kwargs)
1601 def sel(
1602 self,
1603 indexers: Mapping[Any, Any] | None = None,
(...) 1607 **indexers_kwargs: Any,
1608 ) -> Self:
1609 """Return a new DataArray whose data is given by selecting index
1610 labels along the specified dimension(s).
1611
(...) 1715 Dimensions without coordinates: points
1716 """
-> 1717 ds = self._to_temp_dataset().sel(
1718 indexers=indexers,
1719 drop=drop,
1720 method=method,
1721 tolerance=tolerance,
1722 **indexers_kwargs,
1723 )
1724 return self._from_temp_dataset(ds)

File ~/anaconda3/lib/python3.13/site-packages/xarray/core/dataset.py:3001, in Dataset.sel(self, indexers, method, tolerance, drop, **indexers_kwargs)
2933 """Returns a new dataset with each array indexed by tick labels
2934 along the specified dimension(s).
2935
(...) 2998
2999 """
3000 indexers = either_dict_or_kwargs(indexers, indexers_kwargs, "sel")
-> 3001 query_results = map_index_queries(
3002 self, indexers=indexers, method=method, tolerance=tolerance
3003 )
3005 if drop:
3006 no_scalar_variables = {}

File ~/anaconda3/lib/python3.13/site-packages/xarray/core/indexing.py:206, in map_index_queries(obj, indexers, method, tolerance, **indexers_kwargs)
203 options = {"method": method, "tolerance": tolerance}
205 indexers = either_dict_or_kwargs(indexers, indexers_kwargs, "map_index_queries")
--> 206 grouped_indexers = group_indexers_by_index(obj, indexers, options)
208 results = []
209 for index, labels in grouped_indexers:

File ~/anaconda3/lib/python3.13/site-packages/xarray/core/indexing.py:167, in group_indexers_by_index(obj, indexers, options)
165 grouped_indexers[index_id][key] = label
166 elif key not in obj.dims:
--> 167 raise KeyError(
168 f"{key!r} is not a valid dimension or coordinate for "
169 f"{obj.class.name} with dimensions {obj.dims!r}"
170 )
171 elif len(options):
172 raise ValueError(
173 f"cannot supply selection options {options!r} for dimension {key!r}"
174 "that has no associated coordinate or index"
175 )

KeyError: "'isobaric' is not a valid dimension or coordinate for Dataset with dimensions FrozenMappingWarningOnValuesAccess({'x': 614, 'y': 428, 'time1': 1, 'isobaric1': 39})"

Any hint to solve this?

THank you