code: def timeseries(datasets):
      
      expanded_ds = []
      for ds in datasets:
          tmp = ds.expand_dims("time")
          tmp.coords["time"] = pd.DatetimeIndex([ds.attrs["start_time"]])
          expanded_ds.append(tmp)
  
      res = xr.concat(expanded_ds, dim="time")
      res.attrs = combine_metadata(*[x.attrs for x in expanded_ds])
      return res
docs: Expands dataset with and concats by time dimension