% ***********************************************************************
% kaolinite_k.m
% ***********************************************************************
%
% Data description:
% refractive index (imaginary part) of kaolinite (185-2600nm), 
%
% Prepared by: Bradley W. Kimmel on July 2009
%
% ***********************************************************************
%
% Original data source:
% 
% Chapter V: Optical Complex Index of Refraction Between
% 0.185 and 2.5 um, Table V-15, page 105
% in
%
% Egan, W.G. and Hilgeman, T.W., "Optical Properties of Inhomogeneous
% Materials:  Applications to Geology, Astronomy, Chemistry, and
% Engineering", Academic Press, 1979
%
% ***********************************************************************
clear all;

wave = [185, 190, 200, 210, 215, 220, 225, 233, 240, 260, 280, 300, ...
		326, 360, 370, 400, 433, 466, 500, 533, 566, 600, 633, 666, ...
		700, 817, 907, 1000, 1105, 1200, 1303, 1400, 1500, 1600, 1700, ...
		1800, 1900, 2000, 2100, 2200, 2300, 2400, 2500, 2600];

data = 10 .^ [...
		-3.02, -2.98, -2.92, -2.85, -2.90, -2.86, -2.91, -2.92, -2.91, ...
		-2.93, -2.93, -2.97, -3.09, -3.40, -3.56, -3.69, -3.84, -3.91, ...
		-4.02, -4.28, -4.41, -4.42, -4.38, -4.03, -4.00, -3.89, -3.91, ...
		-3.80, -3.87, -3.80, -3.69, -3.46, -3.51, -3.50, -3.39, -3.34, ...
		-3.25, -3.25, -3.16, -2.92, -2.75, -2.80, -2.48, -2.37];

fid = fopen('ecKao.txt', 'wb');
for i = 1:length(data)
	fprintf(fid, '%d  %3.4g \n', wave(i), data(i));
end;
fclose(fid);

plot(wave, data, 'b-', 'linewidth', 3);
axis([185 2600 0 5e-3]);
set(gca, 'FontSize', 14);
grid on;
title('Kaolinite', 'FontSize', 14);
ylabel('extinction coefficient', 'FontSize', 14);
xlabel('wavelength (nm)', 'FontSize', 14);
set(gcf,'PaperUnits','inches','PaperPosition',[0 0 6 5]);
pause;
print('-dpng', '-r100', 'kaolinite_k.png');
close;

