% 
% plot_Hc2_vs_T_5.m
%
% plot_Hc2_vs_T_5.m - Changes from plot_Hc2_vs_T_4.m
% 1. Use SJB_OUTPUTS.dat.
%
% plot_Hc2_vs_T_4.m - Changes from plot_Hc2_vs_T_3.m
% 1. Use OUTPUTv6 - revert to the old Dp/Ds ratios! They were correct!
% 2. Delete all old loading in code
%
% plot_Hc2_vs_T_3.m - Changes from plot_Hc2_vs_T_2.m
% 1. Plot theory from SJBs MathCAD models but force D_pi/D_sigma to be
% constant and equal to Noguchi's values. [6.46 for 77nm, 14.11 for 160nm
% WRONG values - see point 3].
% 2. Remove datapoint at 450 Oe. Poor temperature control.
% 3. OUTPUTv3 - from now on use a Dp/Ds ration of Dp/Ds(c) from Noguchi
% table. Previously used wrong Ds(ab).
% 4. OUTPUTv4 use a Tc=26.32K (average of two "down" temperature sweeps)
% instead of 26.27K (average of two "down" and two "up" sweeps.
% 5. OUTPUTv5 - no differences in model, just changed name as now doing 3
% (4) individual trends as per Simons emial 10/03/14.
% 
%
% plot_Hc2_vs_T_2.m - Changes from plot_Hc2_vs_T.m
% 1. Plot theory from SJBs MathCAD models rather than simple polyfits
% 2. Remove Tc(H) plots
%
t77=load('SJB_OUTPUT77.dat');
t160=load('SJB_OUTPUT160.dat');
%
figure
hold on
%
load MgB2_77nm_Hc2T_data.mat
%title('Noguchi''s Dp/Ds ratios - varying Dp')
Ha=Ha/100;
h1=herrorbar(Tcs_50,Ha,(Tcs_50-Tcs_10),(Tcs_90-Tcs_50),'o');
set(h1(2),'linewidth',2) ;
h2=plot(t77(:,1),t77(:,2)*10,'b','linewidth',2);
%
load MgB2_160nm_Hc2T_data.mat
Ha=Ha/100;
Ha(6)=[];Tcs_10(6)=[];Tcs_50(6)=[];Tcs_90(6)=[];
h3=herrorbar(Tcs_50,Ha,(Tcs_50-Tcs_10),(Tcs_90-Tcs_50),'rs');
h4=plot(t160(:,1),t160(:,2)*10,'r','linewidth',2);
set(h3(2),'linewidth',2) ;
ylabel('H_{c2} (kOe)','fontsize',14)
xlabel('T (K)','fontsize',14)
%legend([h1(2) h2(1) h3(1) h4(1)],{'77nm film','Noguchi fixed ratio - Dp = 4.2E-4 cm^2s^{-1}','Noguchi fixed ratio - Dp = 4.45E-4 cm^2s^{-1}','Noguchi fixed ratio - Dp = 4.2E-4 cm^2s^{-1}'})
legend([h1(2) h3(2)],{'77nm film','160nm film'},'fontsize',14,'location','South')
set(gca,'Fontsize',12)
box on
set(gca,'LineWidth',1,'TickLength',[0.025 0.025]);
axis([24 36 0 10])
%