%%
% Define crystal and sample symmetry
cs_ol = crystalSymmetry('Pnma',[9.915267 5.886196 4.713400])
ss = specimenSymmetry()
% Plotting convention
plotx2east
% File with Euler angles
fname = fullfile('/path/to/folder/',...
    '','euler.txt');

%%
% Load file and create the corresponding list of orientations
ori = load(fname)
n = size(ori);
n = n(1);

for i = 1:n
    phi1 = ori(i,1)*degree;
    Phi = ori(i,2)*degree;
    phi2 = ori(i,3)*degree;
    ori0(i) = orientation('Euler',phi1,Phi,phi2,cs_ol,ss);
end

%%
% Plots grains using a color scale base on ODF fitted to the grains

% Fit an ODF to the data
odf = calcODF(ori0,15*degree)


% Plot grains as IPF of the compression (if compression is in X), colored based on ODF
plotIPDF(ori0,eval(odf,ori0),[vector3d.X],'antipodal','MarkerSize',10, 'MarkerEdgeColor', 'b')
CLim(gcm,[0,10])
mtexColorbar