Btw. I think what FT does with the dics is a bit simpler. You compute the solution for just the one freq (averaging possible) specified and get a regular spatio-temporal map. Would be nice to have both, average across freqs, iterate over time as second dimension and the other way round.