function report_plots(bottle_var,ctd_var1,ctd_sta,ctd_prs,nis,smp_btl,hcolor,datatype,sensor,variable,unit)

deep_idx=find(ctd_prs>=hcolor);

if datatype(1) == 'U',
    type = 'uncal';
else
    type = 'cal';
end;


figure
plot(ctd_sta,bottle_var-ctd_var1,'k+')
grid on;
xlabel('Station number')
ylabel(['Bottle ',variable,' - CTD ', variable,' (',unit,')'])
%ylim(-0.01,0.01)
diff_median=nanmedian(bottle_var-ctd_var1);
aux=find(abs(bottle_var-ctd_var1)<0.002);
disp([type,' Ratio of points with WOCE standard of +/- 0.002'])
length(aux)/length(bottle_var-ctd_var1)
diff_stddev=nanstd(bottle_var-ctd_var1);
title([datatype,' ',sensor,' Sensor - Median Diff: ',num2str(diff_median),' Std Dev: ',num2str(diff_stddev)])

jpeg=['print -deps figures/',type,'_ctd_',variable(1:3),'_st','.eps'];
eval(jpeg)
if type(1) == 'c',
saveas(gcf,['figures/',type,'_ctd_',variable(1:3),'_st.fig'])
end;
clear diff_median diff_stddev

figure  
plot(ctd_prs,bottle_var-ctd_var1,'k+')
grid on;
xlabel('Pressure (dbar)')
ylabel(['Bottle ',variable,' - CTD ', variable,' (',unit,')'])
%ylim(-0.01,0.01)
diff_median=nanmedian(bottle_var-ctd_var1);
diff_stddev=nanstd(bottle_var-ctd_var1);
title([datatype,' ',sensor,' Sensor - Median Diff: ',num2str(diff_median),' Std Dev: ',num2str(diff_stddev)])  

jpeg=['print -deps figures/',type,'_ctd_',variable(1:3),'_pr','.eps'];
eval(jpeg)
if type(1) == 'c',
saveas(gcf,['figures/',type,'_ctd_',variable(1:3),'_pr.fig'])
end;
clear diff_median diff_stddev

%%%%%%%%%% begin plots >hcolor dbar %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

figure
plot(ctd_sta(deep_idx),bottle_var(deep_idx)-ctd_var1(deep_idx),'k+')
grid on;
xlabel('Station number')
ylabel(['Bottle ',variable,' - CTD ', variable,' (',unit,')'])
%ylim(-0.01,0.01)
diff_median=nanmedian(bottle_var(deep_idx)-ctd_var1(deep_idx));
aux=find(abs(bottle_var(deep_idx)-ctd_var1(deep_idx))<0.002);
disp([type,' Ratio of points below',num2str(hcolor), ' dbar with WOCE standard of +/- 0.002'])
length(aux)/length(bottle_var(deep_idx)-ctd_var1(deep_idx))
diff_stddev=nanstd(bottle_var(deep_idx)-ctd_var1(deep_idx));
title([datatype,' Data >',num2str(hcolor),' dbar - Median Diff: ',num2str(diff_median),' Std Dev: ',num2str(diff_stddev)])  
title([datatype,' ',sensor,' Sensor - Below',num2str(hcolor),' dbar - Median Diff: ',num2str(diff_median),' Std Dev: ',num2str(diff_stddev)])  

jpeg=['print -deps figures/',type,'_ctd_',variable(1:3),'_deep_st','.eps'];
eval(jpeg)
if type(1) == 'c',
saveas(gcf,['figures/',type,'_ctd_',variable(1:3),'_deep_st.fig'])
end;
clear diff_median diff_stddev

figure
plot(ctd_prs(deep_idx),bottle_var(deep_idx)-ctd_var1(deep_idx),'k+')
grid on;
xlabel('Pressure (dbar)')
ylabel(['Bottle ',variable,' - CTD ', variable,' (',unit,')'])
%ylim(-0.01,0.01)
diff_median=nanmedian(bottle_var(deep_idx)-ctd_var1(deep_idx));
diff_stddev=nanstd(bottle_var(deep_idx)-ctd_var1(deep_idx));
title([datatype,' ',sensor,' Sensor - Below',num2str(hcolor),' dbar - Median Diff: ',num2str(diff_median),' Std Dev: ',num2str(diff_stddev)])  

jpeg=['print -deps figures/',type,'_ctd_',variable(1:3),'_deep_pr','.eps'];
eval(jpeg)
if type(1) == 'c',
saveas(gcf,['figures/',type,'_ctd_',variable(1:3),'_deep_pr.fig'])
end;

return
