This folder contains ASCII files for the (x,y,error) data and 
for the intensity colormaps in the various figures, the filenames 
contain the figure number and panel label and in some cases additional 
labels to uniquely identify which data they contain. 

(x,y,error) or multi-column data
+++++++++++++++++++++++
fig_1e_0T.xye       3-column [Energy(meV) Intensity Error] 
fig_1e_0.21T.xye    3-column [Energy(meV) Intensity Error]
fig_1e_9T.xye       3-column [Energy(meV) Intensity Error]
fig_1f_5T.xye       3-column [Energy(meV) Intensity Error]
fig_1f_5T_fit.xye   2-column [Energy(meV) Intensity]
fig_1f_9T.xye       3-column [Energy(meV) Intensity Error]
fig_1h_0.21T.xye    3-column [Energy(meV) Intensity Error]
fig_1h_7T.xye       3-column [Energy(meV) Intensity Error]
fig_1h_7T_model.xye 2-column [Energy(meV) Intensity]
fig_1i_1m_lower.xye 4-column [AppField(T) InternalField(T) Energy(meV) 
                              Energy_error(meV)]  
fig_1i_1m_upper.xye 5-column [AppField(T) InternalField(T) Energy(meV) 
                              Energy_error(meV) NeutronEi(meV)]
fig_1i_2m_lower.xye 4-column [AppField(T) InternalField(T) Energy(meV) 
							  Energy_error(meV)] 
fig_1i_2m_upper.xye 5-column [AppField(T) InternalField(T) Energy(meV) 
			                  Energy_error(meV) NeutronEi(meV)]
fig 1j.xye          3-column [InternalField(T) Temperature(K) error_Temperature(K)]

fig_2_0.0T.xye	comma(,) separated 4-column 
            [Temperature(K) Error_Temperature(K) 
			HeatCapacity(J/K/mol_Yb) Error_HeatCapacity(J/K/mol_Yb)]
            same for 0.05T, 0.075T, 0.1T, 0.125T, 0.15T, 0.21T, 0.3T, 0.5T, 1.5T  
fig_2_gap.xye  5-column 
			[Applied_Field(T) Gap(meV) Gap_Error(meV) 
			Internal_Field(T) Internal_Field_Error(T)]

fig_S6a_field_along_001.xye, dispersions data (8-column format) for B=5,9 T ||[0 0 1] 
            [AppliedField(T) H(rlu) K(rlu) L(rlu) Energy(meV) 
			Error_energy(meV) ModeNumber(1-4) Calculated_Energy(meV)]               
fig_S6a_field_along_m110.xye, dispersions data (8-column format) for B=5 T ||[-1 1 0]
            [AppliedField(T) H(rlu) K(rlu) L(rlu) Energy(meV) 
			Error_energy(meV) ModeNumber(1-4) Calculated_Energy(meV)] 
			
fig_S13a_9T.xye  3-column format [Energy(meV) Intensity Error]
            same for 8T, 7T, 6T, 5T, 4T and 3T
fig_S13b.xye 5-column [AppliedField(T) Area ErrorArea ModelArea ModelAreaUncertainty]      

fig_S15_1T.xye 3-column format [Energy(meV) Intensity Error] 
            same for 1.5T, 2T, 3T and 9T

fig_S19_blue.xye comma(,) separated 4-column 
            [Temperature(K) Error_Temperature(K) HeatCapacity(J/K/mol Yb) 
			Error_HeatCapacity(J/K/mol Yb)]
fig_S19_red.xye  2-column [Temperature(K) HeatCapacity(mJ/K/mol Yb)]

fig_S20a_data.xye comma(,) separated 4-column 
            [Temperature(K) Error_Temperature(K) HeatCapacity(J/K/mol Yb) 
			Error_HeatCapacity(J/K/mol Yb)]
fig_S20a_SWT.xye 2-column [Temperature(K) CalculatedHeatCapacity(J/K/mol Yb)]
            same for panels 20b-e

fig_S21_1.8K.xye comma(,) separated 
            3-column [InternalField(T) Magnetization(mu_B/Yb) Error_magnetization]
fig_S21_90mK.xye 2-column [InternalField(T) Magnetization(mu_B/Yb)]
fig_S21_MFT.xye  2-column [Field(T) Calculated_magnetization(mu_B/Yb)]

Colour map figures
++++++++++++++++++++

The files to produce the intensity colour maps with neutron data have some 
header information (axes labels, limmits etc.) and a data table in a 
5-column format [x y Intensity Error SmoothedIntensity] with pixel bin 
coordinates (x,y), intensities (I) and corresponding errors (E), and 
SmoothedIntensities (last column). The intensities (I) are the measured 
raw neutron counts converted into dynamical correlations S(wavevector,energy) 
units minus the non-magnetic background and corrected for neutron absorption 
effects. The last column, SmoothedIntensities, are obtained after smoothing 
the I values between neighbouring bins by onvolving a matrix of weights, A, 
ns times, to have a smoother colour variation between neighbouring pixels in 
the plotted intensity map. For Figs. S1, S2, S3 and S12 the convolving matrix 
is A=[0.1  0.2  0.1; 0.2  0.8  0.2; 0.1  0.2  0.1] and the smoothing convolution 
is performed twice (ns=2). For Figs. S16 and S17 the matrix A is the uniform 
weight 3 x 3 matrix [1 1 1; 1 1 1; 1 1 1] and ns=1. For Figs. S14 and S18 the 
matrix A is the 5 x 5 uniform weight matrix and ns=1. For the simulation panels 
the ascii files are in a 3-column format [x y SmoothedIntensity] where the 
SmoothedIntensities are the calculated intensities smoothed the same way 
as for the data panels.     
  
To obtain the plots in the paper could use the matlab code supplied 
plot_slice_figure(filename);
with the appropriate .dat filename. This code reads the data header with 
all the figure label information, reshapes the data table into a 2D matrix 
format, converts from centre to bin edges, then uses pcolor with apropriate 
axes and intensity limits to plot the intensity map with the supplied 
rainbow colormap coltab.dat (RGB format).  
  
fig_1a.dat up to fig_1g.dat
>>plot_slice_figure('fig_1a.dat'); 		
	
fig_S1a.dat to fig_S1x.dat
fig_S2a.dat to fig_S2x.dat 
fig_S3a.dat to fig_S3x.dat
fig_S12a.dat to fig_S12h.dat 
fig_S14a.dat and fig_S14b.dat 
fig_S16a.dat to fig_S16l.dat 
fig_S17a.dat to fig_S17l.dat 
fig_S18a.dat to fig_S18c.dat 