Archive contains data from figures in the paper arranged as follows. 

* Fig1: Fig_1a.dat contains the data plotted in Fig.1(a). Left column:
magnetic field; right column: magnetisation. 
This repeats for 5K, 10K, 20K, 30K, 40K, 50K, and 60K. (Therefore there are 14 columns.)

* Fig. 1b.dat contains the data plotted in Fig.1(b). It contains the
data from the phase diagram, as well as the inset figure. 
The left column is temperature; the right column is magnetisation. 
It repeats for 150mT and 220mT. (Therefore there are 4 columns).

* Fig.2: files nk.dat (where n is a number) contain the data measured at a temperature n Kelvin. 

* Fig.3;  the fitted parameters are in file fitted_data.dat. The relevant columns are numbers 4-7 (labelled in the file).
The file file knight_shift.dat contains the Knight shift data in the format x, y, y_error. 

Fig 4: A description of the data from simulations appearing in this
figure 4 is given in detail below. 

The archive contains the following directories. All calculations
performed for [111] orientation. 

SkX_0.15T_18nm_1cell   Skyrmion wavelength 18nm, 1 skyrmion unit cell
SkX_0.15T_8.5nm_4cell  Skyrmion wavelength 8.5nm, 4 skyrmion unit cell
rcone_0.15T_18nm_10cell    Cone wavelength 18nm, 10 skyrmion unit cell, Hc =0.5T
rhelix_0.15T_18nm_10cell   Helix wavelength 18nm, 10 skyrmion unit cell

(multiple  adjacent unit  cells have  been  sampled in  some cases  to
improve statistics)

Each directory has a subdirectory (4mTbins) containing histogrammed results.

Calculation output files:
cellmap3D.out			   Plot of spins within the unit cell
dipole.par	   		   Calculation settings
field.scan.abs.units.$number.out   Raw field output for stopping site $number
				   000 are `site 0' points
				   others are `site 1' points
				   (001 were histogrammed as representative)
field.scan.abs.units.total.out	   Raw field output for all stopping sites
output				   Run output

Histogram files:
Bmag.histogram.out		   Histogram output

Format notes:

cellmap3D.out can be plotted using gnuplot's 'with vector' option:
x  y  z  dx  dy  dz

field.scan.abs.units.$X.out:
lattice coords(x,y,z)	cartesian coords(x,y,z) Bfield(x,y,z) |B|

Bmag.histogram.out: 
upper end of bin	unrenormalised frequency	total number of data points
(note that the SkX point 0 and point 1 plots in FIG. 4 have been normalised wrt the 
total number of data points for all points)
