code: def format_index_raw(data):
      
      tz_raw = data.columns[3]
      timezone = TZ_MAP.get(tz_raw, tz_raw)
      year = data.Year.apply(str)
      jday = data.DOY.apply(lambda x: '{:03d}'.format(x))
      time = data[tz_raw].apply(lambda x: '{:04d}'.format(x))
      index = pd.to_datetime(year + jday + time, format="%Y%j%H%M")
      data = data.set_index(index)
      data = data.tz_localize(timezone)
      return data
docs: Create DatetimeIndex for the Dataframe localized to the timezone provided
      as the label of the third column.
  
      Parameters
      ----------
      data: Dataframe
          Must contain columns 'Year' and 'DOY'. Timezone must be found as the
          label of the third (time) column.
  
      Returns
      -------
      data: Dataframe
          The data with a Datetime index localized to the provided timezone.