This folder contains ASCII files for the (x,y,error) data plotted in Figs. 2c, 3d, 4b, 4e, and 
   Supplementary Figs. 1a, 1a inset, 1b, 1c, 1d, 6a, 6b, 6c, 9, 10, 11, 12, and 14
and (x,y) fits plotted in Figs. 1a, 1b, 1c, 1d, 2c, 3d, 4b, 4c, 4d, 4e, 4f, and 
   Supplementary Figs. 1a, 1b, 1c, 1d, 6a, 6b, 6c, 10, 11, 12, 14, 16a, and 16b
and the intensity colour maps in Figs. 1a, 1b, 1c, 1d, 2d, 2e, 2f, 2g, 3c, 4c, 4d, 4f and 
   Supplementary Figs. 8a, 8b, 11, 12, 15a, 15b, 15c, and 15d

The folder also contains (x,y) data - BZ.xy - to plot the Brillouin zone boundaries in the hk plane and 
an RGB colour map in colmap.dat.

The supplied MATLAB code plot_data.m plots 37 figures when run in the same directory as the unzipped data files.

(x,y) data : fits and lines/calculations
++++++++++++++++++++++++++++++++++++++++++++++
Energies are given in meV for modes ordered in increasing energy, and (h,k,l) in rlu.
Angles are guiven in degrees.
Intensities are given in arbitrary units. 
|q| is given in 1/Angstroms and is used as the x coordinate in plots along wave vector paths, which is defined
as the distance travelled along a wave vector path in reciprocal space from the starting wave vector.
d-spacings are given in Angstroms.
The data in the file BZ.xy give the outlines of the Brillouin zone when plotted in Figs. 2d, 2e, 3c, and 
Supplementary Figs. 15a, 15b, 15c, and 15d.

fig_1a.xy			 8 columns  [|q| h k l Energy(mode 1) Energy(mode 2) Energy(mode 3) Energy(mode 4)]
fig_1b.xy		 	 8 columns  [|q| h k l Energy(mode 1) Energy(mode 2) Energy(mode 3) Energy(mode 4)]
fig_1c.xy		 	 6 columns  [|q| h k l Energy(mode 1) Energy(mode 2)]
fig_1d.xy		 	 6 columns  [|q| h k l Energy(mode 1) Energy(mode 2)]
fig_2c.xy		 	 3 columns  [Angle Intensity(top band) Intensity(bottom band)]
fig_3d_eta.xy		 2 columns  [Energy Intensity(XXZeta)]
fig_3d_XXZ.xy		 2 columns  [Energy Intensity(XXZ_only)]
fig_4b.xy		 	 3 columns  [Energy Intensity(both CEF peaks) Intensity(only lower CEF peak)]
fig_4c.xy		 	 6 columns  [|q| h k l Energy(mode 1) Energy(mode 2)]
fig_4d.xy		 	 6 columns  [|q| h k l Energy(mode 1) Energy(mode 2)]
fig_4e.xy		 	 3 columns  [Angle Intensity(top band) Intensity(bottom band)]		
fig_4f.xy		 	 6 columns  [|q| h k l Energy(mode 1) Energy(modes 2)]
SM_Fig_1a_data.dat	 	 2 columns  [d-spacing Intensity] Standard counting errors apply
SM_Fig_1a_fit.dat	 	 2 columns  [d-spacing Intensity]	
SM_Fig_1b_data.dat	 	 2 columns  [d-spacing Intensity] Standard counting errors apply
SM_Fig_1b_fit.dat	 	 2 columns  [d-spacing Intensity]
SM_fig_6a.xy		 3 columns  [Angle Intensity(top band) Intensity(bottom band)]
SM_fig_6b.xy		 3 columns  [Angle Intensity(top band) Intensity(bottom band)]
SM_fig_6c.xy		 3 columns  [Angle Intensity(top band) Intensity(bottom band)]		
SM_fig_10.xy		 2 columns  [Energy Intensity]
SM_fig_11.xy		 8 columns  [|q| h k l Energy(mode 1) Energy(mode 2) Energy(mode 3) Energy(mode 4)]		
SM_fig_12.xy		 12 columns [|q| h k l Energy(mode 1, phi 90) Energy(mode 1, phi 270) 
				     Energy(mode 2, phi 90) Energy(mode 2, phi 270) Energy(mode 3, phi 90)
				     Energy(mode 3, phi 270) Energy(mode 4, phi 90) Energy(mode 4, phi 270)]		
SM_fig_14.xy		 3 columns  [Energy Intensity(negative eta) Intensity(positive eta)]
SM_fig_16a_data.xy	 	 2 columns  [Angle Energy(Zero Point Energy)]
SM_fig_16a_fit.xy	 	 2 columns  [Angle Energy(Zero Point Energy)]
SM_fig_16b_data.xy	 	 2 columns  [Energy(eta) Energy(Lambda)]
SM_fig_16b_fit.xy	 	 2 columns  [Energy(eta) Energy(Lambda)]
BZ.xy			 2 columns  [(0.5x,-x) (y,0)]

(x,y,error) data : data points with errors
++++++++++++++++++++++++++++++++++++++++++++++
Energies are given in meV and (h,k,l) in rlu.
Temperatures are given in K.
Intensities and Errors are given in arbitrary units.
Moments are given in units of Bohr magnetons. 
|q| is given in 1/Angstroms and is used as the x coordinate in plots along wave vector paths, which is defined
as the distance travelled along a wave vector path in reciprocal space from the starting wave vector.

fig_2c.xye		 5 columns [Angle Intensity(top band) Error(Int. top) Intensity(bottom band) Error (Int. bottom)]
fig_3d.xye		 3 columns [Energy Intensity Error(Intensity)]
fig_4b.xye		 3 columns [Energy Intenstiy Error(Intensity)]
fig_4e.xye		 5 columns [Angle Intensity(top band) Error(Int. top) Intensity(bottom band) Error (Int. bottom)]
SM_Fig_1a_inset_data.dat 	 3 columns [Temperature Moment Error(Moment)]		
SM_fig_6a.xye		 5 columns [Angle Intensity(top band) Error(Int. top) Intensity(bottom band) Error (Int. bottom)]
SM_fig_6b.xye		 5 columns [Angle Intensity(top band) Error(Int. top) Intensity(bottom band) Error (Int. bottom)]
SM_fig_6c.xye		 5 columns [Angle Intensity(top band) Error(Int. top) Intensity(bottom band) Error (Int. bottom)]
SM_fig_9.xye		 7 columns [h k l Energy(Observed) Error(Energy Observed) ModeNumber Energy(Calculated=\tilde{\omega})]
				    (ModeNumber = 5 denotes the centre of mass of unresolved modes 3 and 4)
SM_fig_10.xye		 3 columns [Energy Intensity Error(Intensity)]
SM_fig_11.xye		 3 columns [|q| Energy Error(Energy)]
				   First row gives number of points in modes 0, 1, 2, and 5
				   (5 denotes either of nearly degenerate modes 3 and 4, 
				    0 denotes a point between modes 1 and 2 that arises due to the presence of magnetic domains)
SM_fig_12.xye		 3 columns [|q| Energy Error(Energy)]
				   First row gives number of points in modes 0, 1, 2, and 5
				   (5 denotes either of nearly degenerate modes 3 and 4)
SM_fig_14.xye		 3 columns [Energy Intensity Error(Intensity)]

Colour map figures
++++++++++++++++++++++++++++++++++++++++++++++
The first row of the file gives the dimensions of the plot arrays in the form [nRows, nCols].

The intensities are smoothed between neighboring bins by convolving a matrix with relative weights 
A=[0.1  0.2  0.1; 0.2  0.8  0.2; 0.1  0.2  0.1] to have a smoother colour variation between neighboring 
bins in the plotted intensity map. The plotted intensities in single panels with experimental data are 
produced by applying this smoothing matrix once (the calculation plots are unsmoothed), and in the 
multipanel Figures 1a and 4c the smoothing matrix used was A'=[1 1 1; 1 1 1; 1 1 1], applied once.

To obtain the plots in the paper, use the matlab code plot_data.m, running this will plot all 37 figure panels, 
can also run separate sections of the code for each individual figure. The code uses an internal function to 
reshape the data into 2D matrices, convert from bin centers to edges, and then use pcolor with appropriate axes, 
colour limits and colour map.

plot_data.m contains the subroutine 'plot_color_fig' which plots the colour data specified in the first argument,
the second argument specifies whether the data contains a column for errors (i.e. for fig_1a.dat), and the third argument
specifies whether smoothed intensities (default, flag=1), raw intensities (flag=0), or errors (flag=-1) should be plotted.

The RGB colour map in colmap.dat was used for all colour figures.

Energies are given in meV and (h,k,l) in rlu.
Intensities and Errors are given in arbitrary units. 
|q| is given in 1/Angstroms and is used as the x coordinate in plots along wave vector paths, which is defined
as the distance travelled along a wave vector path in reciprocal space from the starting wave vector.

fig_1a.dat		8 columns [|q| h k l Energy Intensity Error(Intensity) Intenstity(Smoothed)]
fig_1b.dat		6 columns [|q| h k l Energy Intensity]
fig_1c.dat		8 columns [|q| h k l Energy Intensity Error(Intensity) Intenstity(Smoothed)]
fig_1d.dat		6 columns [|q| h k l Energy Intensity]
fig_2d.dat		5 columns [(0.5x,-x) (y,0) Intensity Error(Intensity) Intenstity(Smoothed)]
fig_2e.dat		5 columns [(0.5x,-x) (y,0) Intensity Error(Intensity) Intenstity(Smoothed)]
fig_2f.dat		8 columns [|q| h k l Energy Intensity Error(Intensity) Intenstity(Smoothed)]
fig_2g.dat		8 columns [|q| h k l Energy Intensity Error(Intensity) Intenstity(Smoothed)]
fig_3c.dat		5 columns [(0.5x,-x) (y,0) Intensity Error(Intensity) Intenstity(Smoothed)]
fig_4c.dat		8 columns [|q| h k l Energy Intensity Error(Intensity) Intenstity(Smoothed)]
fig_4d.dat		6 columns [|q| h k l Energy Intensity]
fig_4f.dat		8 columns [|q| h k l Energy Intensity Error(Intensity) Intenstity(Smoothed)]
SM_fig_8a.dat		8 columns [|q| h k l Energy Intensity Error(Intensity) Intenstity(Smoothed)]
SM_fig_8b.dat		8 columns [|q| h k l Energy Intensity Error(Intensity) Intenstity(Smoothed)]
SM_fig_11.dat		8 columns [|q| h k l Energy Intensity Error(Intensity) Intenstity(Smoothed)]
SM_fig_12.dat		8 columns [|q| h k l Energy Intensity Error(Intensity) Intenstity(Smoothed)]
SM_fig_15a.dat		5 columns [(0.5x,-x) (y,0) Intensity Error(Intensity) Intenstity(Smoothed)]
SM_fig_15b.dat		5 columns [(0.5x,-x) (y,0) Intensity Error(Intensity) Intenstity(Smoothed)]
SM_fig_15c.dat		3 columns [(0.5x,-x) (y,0) Intensity]
SM_fig_15d.dat		3 columns [(0.5x,-x) (y,0) Intensity]
