def _doy_to_datetimeindex(doy, epoch_year=<NUM_LIT>):
doy = np.atleast_1d(doy).astype('<STR_LIT:float>')<EOL>epoch = pd.Timestamp('<STR_LIT>'.format(epoch_year - <NUM_LIT:1>))<EOL>timestamps = [epoch + dt.timedelta(days=adoy) for adoy in doy]<EOL>return pd.DatetimeIndex(timestamps)<EOL>
Convert a day of year scalar or array to a pd.DatetimeIndex.

Parameters
----------
doy : numeric
    Contains days of the year

Returns
-------
pd.DatetimeIndex