For each figure, numbered as in the paper, the 1st line provides the
name of the corresponding data file.  Each curve or set of data points
in the figure is then represented as a pair of (X,Y) values.  Curves
and data using common X values are grouped together.  Data
associated with the left (right) vertical axis are indicated by "Y1"
("Y2").  Labels and / or line styles used in the plot are in brackets after the data column
name.

The columns in the files, including the column headers, are delimited
by tabs.  The columns are numbered from left to right, beginning at 1.

Figure 1:
       neu_birth.txt
       X: col. 1, "R_lim" [R (m)]
       Y: col. 2, "Z_lim" [Z (m), black dashed]

       X: col. 3, "R_neu" [R (m)]
       Y: col. 4, "Z_neu" [Z (m), red circles]

       X: col. 5, "R_sep" [R (m)]
       Y: col. 6, "Z_sep" [Z (m), blue solid]

Figure 2:
       (a):
       i_dens.txt
       X: col. 1, "psi_n" [psi_n]
       Y: col. 2, "n_i_init" [initial]
       Y: col. 3, "n_i_with" [with neutrals]

       (b)
       temperature_pf95_mod_e5.txt
       X: col. 1, "psi_n" [psi_n]
       Y: col. 2, "T_e (eV)" [T_e]

       ion_temps.txt
       X: col. 1, "psi_n" [psi_n]
       Y: col. 2, "T_i_init" [initial]
       Y: col. 3, "T_i_with" [with neutrals]
       Y: col. 4, "T_i_no_neut" [no neutrals]

Figure 3:
       d2_energy.txt
       X: col. 1, "psi" [psi_n]
       Y: col. 2, "CX_src" [DEGAS2 CX]
       Y: col. 3, "Ion_src" [DEGAS2 ion]

       xgc1_energy.txt
       X: col. 1, "psi" [psi_n]
       Y: col. 2, "CX_src" [XGC1 CX]
       Y: col. 3, "Ion_src" [XGC1 ion]

Figure 4:
       turb_int.txt
       X: col. 1, "psi_n" [psi_n]
       Y: col. 2, "turb_int_with" [with neutrals]
       Y: col. 3, "turb_int_no_neut" [no neutals]

Figure 5:
       eta_i.txt
       X: col. 1, "time" [time (ms)]
       Y: col. 2, "eta_i_with(0.9)" [with neutrals, psi_n = 0.9, blue dashed]
       Y: col. 3, "eta_i_with(0.95)" [psi_n = 0.95, blue solid]

       nn_eta_i.txt
       X: col. 1, "time" [time (ms)]
       Y: col. 2, "eta_i_no_neut(0.9)" [no neutrals, psi_n = 0.9, green dashed]
       Y: col. 3, "eta_i_no_neut(0.95)" [psi_n = 0.95, green solid]

Figure 6:
       (a)
       shear.txt
       X: col. 1, "psi_n" [psi_n]
       Y: col. 2, "ExB_with" [with neutrals]
       Y: col. 3, "ExB_no_neut" [no neutrals]

       (b)
       shear_time.txt
       X: col. 1, "time" [time (ms)]
       Y: col. 2, "ExB_with_plus" [with neutrals, upper]
       Y: col. 3, "ExB_with_minus" [with neutrals, lower]

       nn_shear_time.txt
       X: col. 1, "time" [time (ms)]
       Y: col. 2, "ExB_no_neut_plus" [no neutrals, upper]
       Y: col. 3, "ExB_no_neut_minus" [no neutrals, lower]