I tried this, and well...couldn't get it to work.  The problem is that iterating over `self.variables` returns all variables, including `coordinates` which breaks this 
```python
reduce_dims = [dim for dim in da.coordinates if dim in dims]
```
since a `Coordinate` doesn't have a `coordinates` attribute.