switchtolayout;
mesh;

grid = getresult('CHARGE','grid');

id = pinch(grid.ID);
tri = grid.elements;
mask = find(id == 3);

tri2 = tri(mask,1:3);

grid2 = unstructureddataset('test',grid.x,grid.y,grid.z,tri2);
grid2.addparameter('v',0);
grid2.addattribute('N',grid.N);

extractstructure(grid2);
set("name","new_silicon");