CMB / Lowermost Mantle Joint Tomographic Model

Output model for Muir, Tanaka and Tkalčić

Description of included data files: 

corrmat.dat - correlation matrix between slowness & radius perturbation coefficients; order of coefficients is (0,0), (1,-1), (1,0)....(8,8) for slowness, and then the same for radius
drpowers.dat - power per degree l for radius perturbation, column 1 = l, column 2 = mean, column 3 = 5%ile, column 4 = 95%ile
vppowers.dat - power per degree l for Vp perturbation, column 1 = l, column 2 = mean, column 3 = 5%ile, column 4 = 95%ile (using 13.61 km/s reference velocity)
m_dr.dat - summary statistics for radius perturbation coefficients, column 1 = l, column 2 = m, column 3 = mean, column 4 = 5%ile, column 5 = 95%ile
m_ds.dat - summary statistics for slowness perturbation coefficients, column 1 = l, column 2 = m, column 3 = mean, column 4 = 5%ile, column 5 = 95%ile
spatialcorr.dat - spatial correlation between slowness and radius, column 1 = latitude, column 2 = longitude, column 3 = correlation
tomodata.dat - summary statistics for Vp, column 1 = latitude, column 2 = longitude, column 3 = mean, column 4 = std dev
topodata.dat - summary statistics for radius, column 1 = latitude, column 2 = longitude, column 3 = mean, column 4 = std dev

Spherical harmonics are given by 
Y^m_l(phi, theta) = sqrt(2)*sqrt(((2l+1)(l-m)!)/(4pi(l+m)!)) cos(m phi) P^m_l(cos(theta)); m>0
Y^m_l(phi, theta) = sqrt(2)*sqrt(((2l+1)(l-m)!)/(4pi(l+m)!)) sin(-m phi) P^(-m)_l(cos(theta)); m<0
Y^m_l(phi, theta) = sqrt(((2l+1)(l-m)!)/(4pi(l+m)!)) P^(-m)_l(cos(theta)); m=0

where P^m_l is the associated legendre function including the Condon-Shortley phase