# Generate a list of dates for December of each year from 2000 to 2014
december_dates = pd.date_range(start="2000-12-01", end="2014-12-01", freq="YS-DEC")
# Subset the dataset to include only these December dates and
# plot the 'deficit' variable across the selected time steps
december_subset = subset.sel(time=december_dates)
december_subset.time
<xarray.DataArray 'time' (time: 15)> Size: 120B
2000-12-01 2001-12-01 2002-12-01 2003-12-01 ... 2012-12-01 2013-12-01 2014-12-01
Coordinates: (1)