code: def get_extra_ds(dataset):
      
      ds_collection = {}
      for ds in dataset.attrs.get('ancillary_variables', []):
          ds_collection.update(get_extra_ds(ds))
      ds_collection[dataset.attrs['name']] = dataset
  
      return ds_collection
docs: Get the extra datasets associated to *dataset*.