Each element, consists of the data output from a single run of the model (for input parameters defined in the field ips) and contains the following fields:

pn: [2x17102 double]: Refined resolution mesh node data - see Matlab PDE toolbox documentation.

en: [7x889 double]: Refined resolution mesh edge segment data - see Matlab PDE toolbox documentation.

tn: [4x33854 double]: Refined resolution mesh triangle data - see Matlab PDE toolbox documentation.

p0: [2x2121 double]: Coarse mesh node data - see Matlab PDE toolbox documentation.

e0: [7x385 double]: Coarse mesh edge segment data - see Matlab PDE toolbox documentation.

t0: [4x4098 double]: Coarse mesh triangle data - see Matlab PDE toolbox documentation.

L2: [2121x20 double]: Shear Alfven Wave eigenvalues calculated on the coarse mesh.

ips: [1x1 struct]: Model input parameters - contact A. W. Degeling for more information.

rho: [2121x1 double]: equatorial plasma mass density, defined on the coarse mesh.

u: [342040x1 double]: ULF wave electric field raw output data, defined on the refined mesh. This data consists of the complex amplitudes of cartesian components (in the equatorial plane) for 10 spectral method basis functions (in the field-aligned direction) which are used to define the wave electric field on the equatorial mesh. The equatorial electric field components Ex and Ey are given by simply summing over these 10 amplitudes. For example, Ex and Ey (in mV/m) are obtained by the following Matlab code:

N = 10; nn = 1:N; no_ps = length(pn(1,:)); mVpm = ips.const.mVpm;

uu = reshape(u,no_ps,2.*N);

Ex = sum(uu(:,2.*nn-1),2).*mVpm;

Ey = sum(uu(:,2.*nn),2).*mVpm;

The Matlab function "make_PDP_figs.m" is used to read the data, perform the necessary post-processing operations and output the article figures. To produce all six figures, simply run the file without any input arguments.