%Define the misorientation axis
twindir_cal = Miller(2,0,-2,1,ebsd('Calcite').CS,'UVTW')
%Define the twinning misorientation
twin_cal = orientation('axis',twindir_cal,'angle',78*degree,ebsd('Calcite').CS,ebsd('Calcite').CS)
%Check the angle between boundary misorientations and twinning
%misorientation is small
gbcal = grains.boundary('Calcite','Calcite');
istwinning_cal = angle(gbcal.misorientation,twin_cal) < 5*degree;
%Isolate twins boundary
twinboundary_cal = gbcal(istwinning_cal);
% crystal symmetry
CS = {...
'notIndexed',...
crystalSymmetry('-3m1', [4.99 4.99 17.064], 'X||a', 'Y||b*', 'Z||c', 'mineral', 'Calcite', 'color', 'light blue');