The file name for each data file is composed of "Field"-"figure number"-"time". For example, az-fig1-t1875.txt is the data for A_z, Figure 1, and t=18.75.

az - The vector potential A_z (contour of A_z represents magnetic field lines)
bx, by, bz - Magnetic field B_x, B_y, and B_z
cuz - Current density J_z
ez - Electric field E_z
ni - Number density n_i
ti - Ion temperature T_i
viz - Ion fluid velocity V_iz

The data format is “x-grid”, “y-grid”, “field”. The coordinates x-grid and y-grid are not normalized. To obtain the x and y coordinates, the following normalization should be used: (“x-grid or y-grid” - 1)/40. For example, if the x-grid is 1601, the normalized x grid is (1601-1)/40=40d_i. Field values are already normalized.

For Figure 1, each data file for the 1D panels in (a) and (b) has a name “by-fig1”-“av or y coordinate”-“time”. For example, by-fig1-av-t1875.txt is the data file for the average of B_y at t=18.75, and by-fig1-y10-t2031.txt is for the B_y along y=10 at t=20.31. The data format is “x-grid”, “1025”, “field”. Please ignore the second column “1025”. To obtain the x coordinate, please use the following normalization (“x-grid” - 1)/40. 

For Figure 2, each data file for the 1D panels in (c) has a name “fg-fig2”-“x position”. For example, fg-fig2-x44 is the data file for the x range 44<x/d_i<48. Also, each data file for the 1D panels in (d) has a name “fg-fig2”-“x position”-“y position”. For example, fg-fig2-x44-y10.txt is the data for y=10. The data format is “energy bin”, “2”, “count”. Please ignore the second column “2”. The normalized energy should be calculated as “energy bin -1”*0.005.

For Figure 3 and 4, the file name for the distribution functions is f-"figure number"-"velocity plane"-"panel". For example, f-fig3-vx-vy-b.txt is the file for the distribution function in Fig. 3, in the v_x-v_y plane, and panel (b). The data format is “vx-bin”, “vy-bin”, “count”. The vx-bin and vy-bin should be normalized as (“vx-bin”-99.5)*2 and (“vy-bin”-99.5)*2, respectively.

For Figure 6 and 7, the 3D data are in gda files and hdf5 files. gda files are single precision binary (400x500x200 in d_e electron skin depth scale, where the ion skin depth d_i=10d_e) in column major order.
hdf5 files are for particles (Nx8) where positions are 0,1,2 and velocities are 4,5,6.

For Figure 8, the files xyion-fig8-a.txt and xyion-fig8-c.txt contain the x and y positions of ions in the fifth and the sixth columns, respectively. They are not normalized, and the normalized positions in x and y are obtained as (“x”-1)/40 and (“y”-1)/40. 

For Figure 9, 10, and 11, the data files are p-ion-fig9.txt, p-ion-fig10.txt, and p-ion-fig11.txt. The data format is “time bin”, “x”, “y”, “px”, “py”, “pz”, “gamma”, “particle id”, “ex”, “ey”, “ez”, “bx”, “by”, “bz”. The data are not normalized. Normalized time should be obtained as (“time bin”-1+61000)*0.0003125. The x and y positions should be normalized as (“x”-1)/40 and (“y”-1)/40, respectively. The z position should be obtained as the time integral (using the time step 0.25) of “pz”/“gamma”, normalize by 40. The normalization for the electric field is 0.70710678. The work should be obtained as the time integral (using the time step 0.25) of (electric field)*velocity, where velocity is the averaged velocity (momentum/gamma) at “time bin” and at “time bin”+1. The normalization for the work is 1600.