Load data from EBSD data files.
Add phase vels plot for each layer

And report some layer data. Also switch to
analytical N layer splitting calc for speed
(for now).
Fix length unit for thickness

MSAT works in units of km/s for velocity so
the thicknesses provided to MS_effective_splitting_N
need to be km. Fix this in the calculation of thickness.
Also, I don't think we need to convert from ray to
geographcal reference frame for the step by step
polarizations. We can go back to the Gaussian
Wavelet as this now gives results in a sensible
time.
Make C and rho of ice f(T)

These equations should be good to -30C at least.
added EBSD data