deep=find(oxy_match.ctdprs>500);
diff_median=nanmedian(oxy_match.btloxy-oxy_match.btlomml2mmk1)
diff_stddev=nanstd(oxy_match.btloxy-oxy_match.btlomml2mmk1)
diff_median_deep=nanmedian(oxy_match.btloxy(deep)-oxy_match.btlomml2mmk1(deep))
diff_stddev_deep=nanstd(oxy_match.btloxy(deep)-oxy_match.btlomml2mmk1(deep))

%%%% Begin plotting
subplot(2,1,1)
plot(oxy_match.ctdprs,oxy_match.btloxy-oxy_match.btlomml2mmk1,'g*')
hold on;
plot(oxy_match.ctdprs(deep),oxy_match.btloxy(deep)-oxy_match.btlomml2mmk1(deep),'*')
xlabel('Station (dbar) ','fonts',16)
%ylabel(['Oxygen (Draw T) - Oxygen (CTD-converted) (umol/kg)'],'fonts',16)
title(['Median Diff (Prs > 500 db): ',num2str(diff_median_deep),' Std Dev: ',num2str(diff_stddev_deep) ],'fonts',16)  
legend('Shallow ','Deep (> 500 db) ','Location','Southeast')

subplot(2,1,2)
plot(oxy_match.stnnbr,oxy_match.btloxy-oxy_match.btlomml2mmk1,'g*')
hold on;
plot(oxy_match.stnnbr(deep),oxy_match.btloxy(deep)-oxy_match.btlomml2mmk1(deep),'*')
xlabel('Pressure (dbar) ','fonts',16)
%ylabel(['Oxygen (Draw T) - Oxygen (CTD-converted) (umol/kg)'],'fonts',16)
title(['Median Diff: ',num2str(diff_median),' Std Dev: ',num2str(diff_stddev) ],'fonts',16)  
legend('Shallow ','Deep (> 500 db) ','Location','Southeast')
t=text(-5,1,['O2 (w/Draw T) - O2 (umol/l->umol/kg w/CTD)'],'fonts',16,'rotation',90);
saveas(gcf,['figures/O2_titration_compare.fig'])

figure
plot(oxy_match.btloxy,oxy_match.ctdprs,'k*')
hold on;
plot(oxy_match.btlomml2mmk1,oxy_match.ctdprs,'r*')
ylabel('Pressure')
xlabel('umol/kg')
legend('O2 w/draw T','O2 w/CTD','Location','southwest')
revyaxis
saveas(gcf,['figures/O2_titration_profiles.fig'])
