Data are generated by the C++-language based VPIC simuation, post-processed by Fortran or IDL codes.

pic_data....zip contains PIC simulation data in the binary format.
Files of pic_data_x120-180_y-30-0_**.zip are for the sub-domain of x=120~180 di, y=-30~0 di.
Files of pic_data_x120-180_y30-60_**.zip are for x=120~180 di, y=30~60 di.
For each file, the prefix (e.g., bx, ey, ne, etc.) indicates the quantity. The number in the suffix indicates the time step in units of omega_pe^-1, where 3733 omega_pe^-1 corresponds to 0.5 omega_ci^-1. 
Each file contains float data arrays with a size of 3136x1568x1, corresponding to x-y dimensions. Data can be read by softwares like IDL, python, Matlab, etc, using the standard data reading method.

The following is an example of an IDL program to read a data file.

------------------
  
datafile='Az_119456.gda'
openr, lun1,datafile,/get_lun
nx=3136L
ny=1568L
nz=1L
struct = {data:fltarr(nx,ny,nz)}
field = assoc(lun1,struct)
struct=field[0]
fulldata=reform(temporary(struct.data))
  
------------------

fermi_dat.zip contains ascii data for the Fermi decomposition calculations. Each file is for one event. The 3 columns are for time (in unit of omega_ci^-1), Fermi_xy, and Fermi_z terms. Fermi_xy and Fermi_z are for the unit area, and the total Fermi term would be equal to Fermi_xy+Fermi_z.

je_te_dat.zip contains ascii data for Te and je.E decompositions. Each file is for one event, where the region is shown in the filename. The columns are:
time (omega_ci^-1)
number of X-lines in the region,
total je.E (integrated over the area),
jepara.Epara,
jeperp.Eperp,
sum of the perpendicular jeperp.Eperp decompositions (sum of the following 4 terms),
Fermi term,
Betatron term,
magnetization term,
demagnetized term (jeperp.Eperp values in regions with Kappa<1),
average Te over the area,
delta_Te/miVA^2 where delta_Te is the average Te over the area subtracting the minimum Te in the area, and VA is based on the average |B| and n in the region
delta_Te/miVA^2, where delta_Te is the same as above, and VA is based on the magnetic field amplitude in the x-y plane
average Te over 2dix2di surrounding X-lines
delta_Te/miVA^2, where detal_Te is the average Te over 2dix2di subtracting the minimum Te in the region, and VA is based on average |B| and n in the region
average Te over 2dex2de surrounding X-lines
delta_Te/miVA^2, where delta_Te is the average Te over 2dex2de subtracting the minimum Te in the region, and VA is based on the magnetic field amplitude in the x-xy plane
delta_Te/miVA^2, where delta_Te is the average Te over 2dex2de subtracting the minimum Te in the region, and VA is based on the average |B| and n in the region.
average |B| in the region
average magnetic field amplitude in the x-y plane
average electron density n
minimum Te in the region


je_te_norx_dat.zip contains ascii data of Te and je.E for non-reconnecting current sheets. The format is the same with je_te_dat.zip

thickness_data_plot.zip contains images for individual current sheets at the time with a minimum thickness. The filename contains information about the current sheet. For example, cs_te_twci15.5_x45.292_y8.184_dte-1.07_d3.2_dve6.1_jm0.23.png
means it is at the time t*omega_ci=15.5, the x-line location is at x=45.292 di, y=8,184di, delta_Te/mi_VA^2 is -1.07%, where we take a cut along N across the X-line, delta_Te is the difference between the average Te in the current sheet and that just outside of the current sheet, and VA is based on the inflow parameters outside of the current sheets on both sides of the current sheet. The thickness of the current sheet is 3.2 de. The VeL shear flow at the edges of the current sheet is 6.1 VA, and the maximum |jz| is 0.23.