clf
directory0 = ...
'C:\Users\Alex\Dropbox\main\OUTPUT_LDN\fig_1_2_3_4_5_6_M0_-93_fullspectrum\';

name='TI_PDM_DATA';

limits=[-.4 .4];
k=[0];
length=size(k,2);

for i=1:1:length

filename = [directory0,name,num2str(k(i)),'.csv'];
M2 = csvread(filename);
filename = [directory0,'energies_',name,num2str(k(i)),'.csv'];
M1 = csvread(filename);
filename = [directory0,'bandstructure_',name,num2str(k(i)),'.csv'];
M = transpose(csvread(filename));
filename = [directory0,'Jy_',name,num2str(k(i)),'.csv'];
M3 = csvread(filename);

% % %here we check for artificial zeros%%%%%

for m=1:size(M,1)
for l=1:size(M,2)
   if M(m,l)==0
      M(m,l) = 99999; 
   end
end
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%PART 1 FOR STATE VECTOR

subplot(3,length,i+length);
% if k(i)<13
%     lim = size(M1,2)/2;
%     diff=25;
% elseif k(i) <21
%     lim = 2;
%     diff=25;
% else
%     lim = 3;
%     diff=25;
% end
% lim=size(M2,2)/4;
lim=size(M1,2)/2;
for j=lim:1:lim

    hold on
%if(M1(j)>0.2198 && M1(j)<0.022)
% hFig2 = plot(0:1:100,M2(:,4*j-3) + M2(:,4*j-2)...
%     +M2(:,4*j-1)+M2(:,4*j) + M1(j));
% end
hFig2 = plot(0:1:100,M2(:,4*j-3),'-s'); 
    plot(0:1:100,M2(:,4*j-2),'-o');
    plot(0:1:100,M2(:,4*j-1),'-x');
    plot(0:1:100,M2(:,4*j),'-.');
%end
end

% ylim(limits);
% ylim([0 0.06]);
xlim([0 100]);
set(gca,'XTick', [0 20 40 60 80 100]);
set(gca, 'XTickLabel', ...
    {'0','20','40','60','80','100'})
% % if (k(i)>0)
%     line([20 20],[-10 10],'LineStyle',':','Color','k');

%     line([50.5 50.5],[-1 1],'LineStyle',':','Color','k');
%     line([51 51],[-1 1],'LineStyle',':','Color','k');
% end
% title(['$M_0 = ',num2str(0.28 - k(i) * 0.56/30.0,'%.3f'),'$'],...
%     'Interpreter', 'LaTex');
% title(['B-field $= ',num2str(0 + k(i) * 30,'%.f'),'$ T'],...
%     'Interpreter', 'LaTex');
% title('B-field $= 3000$ T',...
%     'Interpreter', 'LaTex');
% title(['k $= ',num2str(0 + k(i) * 0.05/30.0,'%.4f'),'\mathrm{A}^-1 $'],...
%     'Interpreter', 'LaTex');
ylabel('$|\Psi|^2$','Interpreter', 'LaTex');
xlabel('\AA','Interpreter', 'LaTex');
set(gca,'TickLabelInterpreter','LaTex');
set(gca, 'FontSize', 16);
hold off

% legend({['$B=', num2str(2*2.26 - k(1)*2*2.25/30.0,'%.3f'),'$'],...
%     ['$B=', num2str(2*2.26 - k(2)*2*2.25/30.0,'%.3f'),'$'],...
%     ['$B=', num2str(2*2.26 - k(3)*2*2.25/30.0,'%.3f'),'$'],...
%     ['$B=', num2str(2*2.26 - k(4)*2*2.25/30.0,'%.3f'),'$']},...
%     'Interpreter','latex','Location','northwest');
    legend({'$|P1^+_-,\frac{1}{2}\rangle$',...
        '$|P2^-_+,\frac{1}{2}\rangle$',...
        '$|P1^+_-,-\frac{1}{2}\rangle$',...
        '$|P2^-_+,-\frac{1}{2}\rangle$'},'Interpreter','latex',...
        'Location','north');
    %,'Orientation','Horizontal');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%PART 2 FOR BANDSTRUCTURE
% r=i+1;
subplot(3,length,i);
hFig1 = plot(0:1:100,M,'-wo','MarkerEdgeColor','k');
line([0 100],[0 0],'LineStyle',':','Color','k');
% line([62.5 62.5],[-10 10],'LineStyle','-','Color','r');
ylim(limits);
% xlim([0 100]);
set(gca,'XTick', [0 20 40 60 80 100]);
set(gca, 'XTickLabel', ...
    {'-0.8','-0.48','-0.16','0.16','0.48','0.8'})
ylabel('Energy/eV','Interpreter', 'LaTex');
xlabel('$k_y/\mathrm{nm}^{-1}$','Interpreter', 'LaTex');
set(gca,'TickLabelInterpreter','LaTex');
set(gca, 'FontSize', 16);
hAxes=gca;
% title(['$B = ',num2str(k(i)),'$~T'],...
%     'Interpreter', 'LaTex');
%%%%%


%%PART 3 SPIN%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
subplot(3,length,i+2*length);
% if k(i)<13
%     lim = size(M1,2)/2;
%     diff=25;
% elseif k(i) <21
%     lim = 2;
%     diff=25;
% else
%     lim = 3;
%     diff=25;
% end
% lim=size(M2,2)/4;
% diff=25;
for j=lim:1:lim
    hold on
% hFig2 = plot(0:1:100,M2(:,4*j-3) + M2(:,4*j-2)...
%     +M2(:,4*j-1)+M2(:,4*j));
%if(M1(j)>0.2198 && M1(j)<0.022)
xs=transpose(1:1:size(M,1));
ys=zeros(size(M,1),1);%(M2(:,4*j)+M2(:,4*j-1)+M2(:,4*j-2)+M2(:,4*j-3));
us=zeros(size(M3,1),1);

vs=M3(:,j);
% if(M1(j)>-0.19 && M1(j)<-0.105)
quiver(xs,ys,us,vs,'Marker','.','AutoScaleFactor',1e-2);
%end
hold off
end
ylim([-0.1 0.1]);
xlim([0 100]);
set(gca,'XTick', [0 20 40 60 80 100]);
set(gca, 'XTickLabel', ...
    {'0','20','40','60','80','100'})
ylabel('$\langle \hat{\sigma}_z \rangle$',...
    'Interpreter', 'LaTex');
xlabel('\AA','Interpreter', 'LaTex');
set(gca,'TickLabelInterpreter','LaTex');
set(gca, 'FontSize', 16);
%%%%%%%%%%
% 
% % r=0;
% % saveas(gcf, [directory0,name,num2str(k),'.png']);
end
% subplot(3,1,2);
%     line([81 81],[-10 10],'LineStyle',':','Color','k');