Make MS_expand work for orthorhombic case

This needs 9 constants and only fills in
another 3 but we may as well include it.

Also add a test case.
Merge branch 'master' of github.com:andreww/MSAT
Create MS_save to save Cij to a file

This does the inverse of MS_load but we
still need to add tests and support for
the ematrix format. Also remove a spurious
semi-colon from MS_load.
Allow MS_TI to handle eta

In MS_TI we Berkeley VTI paramaterisation ('panning'), but not the
more general VTI one.  Panning and Romanowicz in their tomography
models assume eta≈1, so you don't get the constants you expect if
you allow it to vary.

Add a routine for general VTI with all five parameters into MS_TI
activated with 'global'.
Add options to finess plots to MS_plot

We can now change the contours and FSW lines
to make nicer plots.
Added functionality for proper treatment of exactly parallel or perpendicular
layers in the stack of operators.

- added a function 'aggregate' which pre-processes the stack of operators, and
  amalgamates wherever possible, prior to effective splitting calculation.

- added appropriate tests into the test suite.
Merge branch 'master' of ssh://majorite/geophysics/walker/git-repos/msat
Make a start on the updated release notes.
Put units in docs for MS_effective_splitting_N

Hz, seconds and degrees.
Import files from external working dir

These are the scripts used to do the D''
topotaxy work.
Update comments for topotaxy example
Update relase notes for this branch
Allow contour bounds to be set in MS_sphere
Added files
Functional prototype
Apparently working prototype with tests
Symmetrical mode.

- Added a symmetrical misfit mode. This runs the lam2 misfit process with the two
  operators in both directions, and takes the average misfit of the two.
Initiated a splitting misfit example.
First stab at effective splitting using wavelet

Use [fast tlag] = MS_effective_splitting_N(0.1, ...
 55.0, [55.0 10.0], [1.0 1.0], 'mode', 'gaussianwavelet')
to do effective splitting calcualtion using this new
method.
Add tests, unwind output

Make the output consistent with the S&S method. Add tests
following some of those from S&S. Note that in three cases
(tests commented out) the result is a few degrees different.
Shows up more for longer periods. Is this close to a faliure
mode for S&S or is something wrong with the wavelets?