#
# Get data from the simulation
#
?"Fetching grid data";
tri = getdata("CHARGE","charge","elements");
vtx3 = getdata("CHARGE","charge","vertices");

svtx = size(vtx3);
nvtx = svtx(1);

stri = size(tri);
ntri = stri(1);

selectXY = 0;
simOrientation = getnamed("CHARGE simulation region","dimension"); 
if (simOrientation == '2D Z-Normal') {
  selectXY = [1,2];
} else { if (simOrientation == '2D Y-Normal') {
  selectXY = [1,3];
} else { if (simOrientation == '2D X-Normal') {
  selectXY = [2,3];
} } }

vtx = vtx3(1:nvtx,selectXY);