function plotgrainpluspf_lucia(ebsd, vecm, unindexed, phasename, minangleboundary)
rot=rotation('axis', zvector, 'angle', 90*degree');
ebsd=rotate(ebsd,rot);
grains=calcGrains(ebsd,'angle',minangleboundary*degree);
vecmiller=vector3d(vecm);
Cristalsimetria = get(grains(phasename),'CS')
vecmiller=Miller(vecmiller,Cristalsimetria);
plotGrains_lucia(vecmiller, unindexed, grains,'colorcoding','tsl');