Loop over freqs
Add some docs and some optional arguments

And the framework for more.
Graph of the output data
Make it go faster

Use glacial_splitting('mode', 'GaussianWavelet') to
recover previous behaviour. But press go before
a coffee break.