Report splitting parameters of each layer

Also resolve a merge conflict on the move
of the splitting parameters outside the
loop in prep for graphing splitting as
a function of pol and freq.
Seperate out the splitting calc.
Loop over source polarizations
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.
Tidy up code and output

* Report input parameters.
* Reformat some output data.
* Reduce line length of long lines (80 char maximum is useful).
* Introduce spaces aroundd binary ops (+, / etc)
  but not uniary ops (so "a = -1" or "a = 0 - 1").
* Introduce spaces between fucntion arguments,
  but not array indicies.
* Use vertical whitespace to seperate sections.
Start work on imaging models

Still lots of work to do on this, but we should have
SKS, icequake and local S and local P, I think.
Add 'quiet' to glacal_splitting ex

This suppresses blow by blow splitting parameter
output.
Make it possible to distiguish freq in plots

And add a legend
A test case for Euler to elasticity calc

This checks we get the same result as Geoff's
calculation using his (fortran?) executables.
Changed some plot options
Merge branch 'glacial-imaging-models' of git://github.com/andreww/MSAT into glacial

Conflicts:
	examples/glacial_splitting/glacial_splitting.m