You should be able to do it with a single time point, it just needs to be Nx1, so doing
`return SourceEstimate(source_power[:, np.newaxis], vertno, fmin, tstep, subject=subject)` should work without the line using `np.repeat`.