% ***********************************************************************
% water_extinction.m
% ***********************************************************************
%
% Data description:
% The specific absorption coefficient of water in [um^-1] (380 - 1000nm) 
% with a resolution of 5nm:
%
% Prepared by: Daniel Yim on December 2011
%
% ***********************************************************************
%
% Original data source:
%
%     K.F. Palmer and D. Williams. Optical properties of water in the 
%     near infrared. J. Opt. Soc. Am., 64:1107-1110, August 1974.
%
% ***********************************************************************
M = [ ...
  380.0, 0.00000001130000000000; ...
  385.0, 0.00000000940000000000; ...
  390.0, 0.00000000850000000000; ...
  395.0, 0.00000000810000000000; ...
  400.0, 0.00000000660000000000; ...
  405.0, 0.00000000530000000000; ...
  410.0, 0.00000000470000000000; ...
  415.0, 0.00000000440000000000; ...
  420.0, 0.00000000450000000000; ...
  425.0, 0.00000000480000000000; ...
  430.0, 0.00000000490000000000; ...
  435.0, 0.00000000530000000000; ...
  440.0, 0.00000000630000000000; ...
  445.0, 0.00000000750000000000; ...
  450.0, 0.00000000920000000000; ...
  455.0, 0.00000000960000000000; ...
  460.0, 0.00000000980000000000; ...
  465.0, 0.00000001010000000000; ...
  470.0, 0.00000001060000000000; ...
  475.0, 0.00000001140000000000; ...
  480.0, 0.00000001270000000000; ...
  485.0, 0.00000001360000000000; ...
  490.0, 0.00000001500000000000; ...
  495.0, 0.00000001730000000000; ...
  500.0, 0.00000002040000000000; ...
  505.0, 0.00000002560000000000; ...
  510.0, 0.00000003250000000000; ...
  515.0, 0.00000003960000000000; ...
  520.0, 0.00000004090000000000; ...
  525.0, 0.00000004170000000000; ...
  530.0, 0.00000004340000000000; ...
  535.0, 0.00000004520000000000; ...
  540.0, 0.00000004740000000000; ...
  545.0, 0.00000005110000000000; ...
  550.0, 0.00000005650000000000; ...
  555.0, 0.00000005960000000000; ...
  560.0, 0.00000006190000000000; ...
  565.0, 0.00000006420000000000; ...
  570.0, 0.00000006950000000000; ...
  575.0, 0.00000007720000000000; ...
  580.0, 0.00000008960000000000; ...
  585.0, 0.00000011000000000000; ...
  590.0, 0.00000013510000000000; ...
  595.0, 0.00000016720000000000; ...
  600.0, 0.00000022240000000000; ...
  605.0, 0.00000025770000000000; ...
  610.0, 0.00000026440000000000; ...
  615.0, 0.00000026780000000000; ...
  620.0, 0.00000027550000000000; ...
  625.0, 0.00000028340000000000; ...
  630.0, 0.00000029160000000000; ...
  635.0, 0.00000030120000000000; ...
  640.0, 0.00000031080000000000; ...
  645.0, 0.00000032500000000000; ...
  650.0, 0.00000034000000000000; ...
  655.0, 0.00000037100000000000; ...
  660.0, 0.00000041000000000000; ...
  665.0, 0.00000042900000000000; ...
  670.0, 0.00000043900000000000; ...
  675.0, 0.00000044800000000000; ...
  680.0, 0.00000046500000000000; ...
  685.0, 0.00000048600000000000; ...
  690.0, 0.00000051600000000000; ...
  695.0, 0.00000055900000000000; ...
  700.0, 0.00000062400000000000; ...
  705.0, 0.00000070400000000000; ...
  710.0, 0.00000082700000000000; ...
  715.0, 0.00000100700000000000; ...
  720.0, 0.00000123100000000000; ...
  725.0, 0.00000148900000000000; ...
  730.0, 0.00000186000000000000; ...
  735.0, 0.00000238000000000000; ...
  740.0, 0.00000254800000000000; ...
  745.0, 0.00000271600000000000; ...
  750.0, 0.00000279000000000000; ...
  755.0, 0.00000276500000000000; ...
  760.0, 0.00000268200000000000; ...
  765.0, 0.00000258600000000000; ...
  770.0, 0.00000249200000000000; ...
  775.0, 0.00000240500000000000; ...
  780.0, 0.00000231700000000000; ...
  785.0, 0.00000223800000000000; ...
  790.0, 0.00000216200000000000; ...
  795.0, 0.00000208800000000000; ...
  800.0, 0.00000202500000000000; ...
  805.0, 0.00000196300000000000; ...
  810.0, 0.00000192700000000000; ...
  815.0, 0.00000193300000000000; ...
  820.0, 0.00000199000000000000; ...
  825.0, 0.00000240900000000000; ...
  830.0, 0.00000282800000000000; ...
  835.0, 0.00000319300000000000; ...
  840.0, 0.00000347500000000000; ...
  845.0, 0.00000375700000000000; ...
  850.0, 0.00000395200000000000; ...
  855.0, 0.00000408900000000000; ...
  860.0, 0.00000422500000000000; ...
  865.0, 0.00000443600000000000; ...
  870.0, 0.00000469600000000000; ...
  875.0, 0.00000495600000000000; ...
  880.0, 0.00000525300000000000; ...
  885.0, 0.00000557500000000000; ...
  890.0, 0.00000589700000000000; ...
  895.0, 0.00000626600000000000; ...
  900.0, 0.00000670700000000000; ...
  905.0, 0.00000714700000000000; ...
  910.0, 0.00000775900000000000; ...
  915.0, 0.00000905300000000000; ...
  920.0, 0.00001035000000000000; ...
  925.0, 0.00001164000000000000; ...
  930.0, 0.00001363000000000000; ...
  935.0, 0.00001580000000000000; ...
  940.0, 0.00001930000000000000; ...
  945.0, 0.00002380000000000000; ...
  950.0, 0.00002980000000000000; ...
  955.0, 0.00003667000000000000; ...
  960.0, 0.00004412000000000000; ...
  965.0, 0.00004827000000000000; ...
  970.0, 0.00005023000000000000; ...
  975.0, 0.00005106000000000000; ...
  980.0, 0.00005020000000000000; ...
  985.0, 0.00004855000000000000; ...
  990.0, 0.00004690000000000000; ...
  995.0, 0.00004425000000000000; ...
  1000.0, 0.00004160000000000000];

plot(M(:,1), M(:,2), 'b', 'linewidth',1.5);
set(gca,'FontSize',14);
title('water specific absorption coefficient','FontSize',14);
ylabel('specific absorption coefficient (um^{-1})','fontsize',14);
xlabel('wavelength (nm)','fontsize',14);
axis([380 1000 0 0.6*10^-4]);
set(gcf,'PaperUnits','inches','PaperPosition',[0 0 6 5]);
grid on;

pause;
print -dpng water_extinction.png -r100

close;
