diff --git a/xarray_eopf/filter.py b/xarray_eopf/filter.py
index 675ce36..07b305c 100644
--- a/xarray_eopf/filter.py
+++ b/xarray_eopf/filter.py
@@ -17,9 +17,8 @@ def filter_dataset(
         return dataset
     name_filter = NameFilter(includes=variables)
     names = set(map(str, dataset.variables.keys()))
+    # find all dataset variables including respective coordinates
     drop_names = names - set(name_filter.filter(names))
     if drop_names:
-        # TODO: also drop now unused coordinates + dimensions as
-        #  they remain even if no longer referenced by any data variables
         dataset = dataset.drop_vars(drop_names)
     return dataset