figure
x1 = [10.5000000000000;12.5000000000000;14.5000000000000;16.5000000000000;18.5000000000000;20.5000000000000;22.5000000000000];
y0 = [0.609816000000000;0.403299000000000;0.297548000000000;0.235849000000000;0.197940000000000;0.176118000000000;0.168769000000000];
y1 = [0.687947000000000;0.495655000000000;0.385019000000000;0.311475000000000;0.261697000000000;0.232125000000000;0.223179000000000];
y2 = [0.699088000000000;0.544312000000000;0.457475000000000;0.394991000000000;0.349956000000000;0.321353000000000;0.314208000000000];
y3 = [0.701463000000000;0.562221000000000;0.497764000000000;0.452815000000000;0.421541000000000;0.401312000000000;0.398915000000000];
y4 = [0.695163000000000;0.570071000000000;0.516115000000000;0.482496000000000;0.462676000000000;0.456812000000000;0.458639000000000];
y5 = [0.825323000000000;0.663675000000000;0.608662000000000;0.574045000000000;0.553555000000000;0.547896000000000;0.561135000000000];


plot(x1,y0*100,'vk',x1,y1*100,'xk',x1,y2*100,'dk', x1,y3*100,'ok', x1,y4*100,'sk',x1,y5*100,'+k','LineWidth',1.2)
strx='$$x/h$$';
xlabel(strx,'interpreter','latex','FontSize',18);
stry='$$\psi^{u*}_{tot}(\%Q)_{y=1h}$$';
ylabel(stry,'interpreter','latex','FontSize',18);
legend({'$$Ri=0$$','$$Ri=0.21$$','$$Ri=0.5$$','$$Ri=0.7$$','$$Ri=1$$','$$Ri=1^*$$'},'interpreter','latex','FontSize',14,'Location','southwest');
legend boxoff;
set(gca,'XTickLabel',num2str(get(gca,'XTick').'));
ylim([0 100]);
set(gca, 'YScale', 'log');
ytickformat('percentage');