A detailed description of the physical meaning of the plotted data can
be found in the figure captions of the arxiv version
https://arxiv.org/abs/1511.03064 of the peer-reviewed paper entitled
"Grand-canonical condensate fluctuations in weakly interacting Bose-Einstein condensates of light"
(authors Christoph Weiss and Jacques Tempere)

Fig1: generated by 
 gnuplot < GrandGanonicalWeissTempere2016Fig1.gsc

uses data files "run0ww100N2000Nmax.dat", "GCE2b.dat", "GCE3b.dat",
         "GCE1b.dat", "GCE3bnew.dat", "GCE1bnew.dat", "GCE2bnew.dat"
In the file names, 1/2/3 refers to the three versions of the LTI
         model. "new" is the weaker interaction.

In the figure:
the black line is the 8th column of
run0ww100N2000Nmax.dat as a function of the first column
the green curve is the analytic equation (x*0.5*0.312*10**8)**(0.5) 
The three thick magenta curves lie on top of each other, shown is the fifth column
of GCE1b.dat, GCE2b.dat GCE3b.dat as a function of the first column.
The three thin magenta curves also lie on top of each other, shown is the fifth column
of GCE1bnew.dat, GCE2bnew.dat GCE3bnew.dat as a function of the first column.

Fig2: generated by
gnuplot < GrandGanonicalWeissTempere2016Fig2.gsc 
uses data files "BDLfig1-1.dat" and "BDLfig1-2.dat". Shown is the
second column as a function of the first.


Fig 3: generated by 
gnuplot < GrandGanonicalWeissTempere2016Fig3.gsc 
uses the datafiles "line0.dat", "line1.dat", "line2.dat", "line3.dat", "line4.dat", "line5.dat". Shown is the
second column as a function of the first.

Fig 4: generated by
 gnuplot < GrandGanonicalWeissTempere2016Fig4.gsc 
uses the datafiles "panel1.dat", "panel2.dat", "BDLfig1-1.dat",
 "panel3.dat".

in panel a, 
the black dashed curve is the second column of panel1.dat
 as a function ot the first.
the black dotted curve is the fourth column of panel1.dat as a
 function of the first.
the blue curve is the third column of panel1.dat as a
 function of the first.

in panel b,
the black curve is the  second column of panel2.dat
 as a function ot the first.
the blue solid curve is the third column of panel2.dat
 as a function ot the first.
the blue dashed curve is the fourth column of panel2.dat
 as a function ot the first.

in panel c,
the black curve is the fifth column of panel3.dat as a function of the
first.
the blue solid line is the sixth column of panel3.dat as a function of the
first.
 the dash-dotted blue line is the fourth column of panel3.dat as a function of the
first.

Fig 5: generated by
gnuplot < GrandGanonicalWeissTempere2016Fig5.gsc 

data files:
"run0ww100N2000Nmax.dat", "GCE2b.dat", "GCE3b.dat", "paperrunww100N.dat",
         "GCE1b.dat", "GCE3bnew.dat", "GCE1bnew.dat", "GCE2bnew.dat", "GCE3.dat"

the black curve shows the third column of run0ww100N2000Nmax.dat as a
function of the first.
the brown curve shows the fifth column of paperrunww100N.dat as a
function of the first.
the thin light blue curves are the fourth columns of
GCE1bnew.dat/GCE2bnew.dat/GCE3bnew.dat as a function of the first
column.
the thick light blue curves are the  fourth columns of
GCE1b.dat/GCE2b.dat/GCE3b.dat as a function of the first
column.


Fig 6: generated by:
gnuplot < GrandGanonicalWeissTempere2016Fig6.gsc  
uses data files:
"BDLfig3-1.dat" "BDLfig3-2.dat" "BDLfig3-2.dat", shown is the second
column as a function of the first.

Fig 7 generated by:
gnuplot < GrandGanonicalWeissTempere2016Fig7.gsc 
uses files:
"GCE1b.dat", GCE2b.dat",  "GCE3b.dat",  "GCE1bnew.dat",
"GCE2bnew.dat", "GCE3bnew.dat"
shown is the sixth column as a function of the first.
