code: def _load_dataset_with_area(self, dsid, coords):
          
          file_handlers = self._get_file_handlers(dsid)
          if not file_handlers:
              return
  
          area = self._load_dataset_area(dsid, file_handlers, coords)
  
          try:
              ds = self._load_dataset_data(file_handlers, dsid)
          except (KeyError, ValueError) as err:
              logger.exception("Could not load dataset '%s': %s", dsid, str(err))
              return None
  
          if area is not None:
              ds.attrs['area'] = area
              calc_coords = (('x' not in ds.coords) or('y' not in ds.coords)) and hasattr(area, 'get_proj_vectors_dask')
              if calc_coords and hasattr(area, 'get_proj_vectors'):
                  ds['x'], ds['y'] = area.get_proj_vectors()
              elif calc_coords:
                  # older pyresample with dask-only method
                  ds['x'], ds['y'] = area.get_proj_vectors_dask(CHUNK_SIZE)
          return ds
docs: Loads *dsid* and it's area if available.