

% Wavelengths for the Satlantic HyperPro HSE187c stored in text file
satwl = load('wl-HSE187c.txt')'

% Satlantic LOG files (converted with SatCon)
files = {'blue1-HSE187c.dat', 400;
		'blue2-HSE187c.dat', 440;
		'blue3-HSE187c.dat', 488;
		'green1-HSE187c.dat', 510;
		'green2-HSE187c.dat', 532;
		'green3-HSE187c.dat', 595;
		'red1-HSE187c.dat', 676;
		'red2-HSE187c.dat', 715;
		'red3-HSE187c.dat', 660,
		'ecoblue-HSE187c.dat', 440,
		'ecored-HSE187c.dat', 670} 

es_all = [];

for kfil = 1:size(files,1)
	
	% Skip header line, Es data is in columns 4:140
	es = textread(files{kfil,1}, '', 'headerlines', 1);
	es = es(:,4:140);
	es = es - repmat(min(es,[],2),1,size(es,2))
	es = es ./ repmat(max(es,[],2),1,size(es,2));
	
	es_all = [es_all; median(es)];
end


% Set default plotting parameters
set(0, 'DefaultLineLineWidth', 0.75)
set(0, 'DefaultLineMarkerSize', 10)
set(0, 'DefaultAxesLineWidth', 0.75)
set(0, 'DefaultTextFontSize', 10)
set(0, 'DefaultAxesFontSize', 10)

% ------ Plot for bb-9
colors_bb9 = jet(9);

figure, hold on
set(gcf, 'paperunits', 'inches', 'papersize', [6.5 4.0])
set(gcf, 'paperposition', [0.25 0.25 6.0 3.5])

	plot(satwl, es_all(1,:), '-', 'Color', colors_bb9(1,:))
	plot(satwl, es_all(2,:), '-', 'Color', colors_bb9(2,:))
	plot(satwl, es_all(3,:), '-', 'Color', colors_bb9(3,:))
	plot(satwl, es_all(4,:), '-', 'Color', colors_bb9(4,:))
	plot(satwl, es_all(5,:), '-', 'Color', colors_bb9(5,:))
	plot(satwl, es_all(6,:), '-', 'Color', colors_bb9(6,:))
	plot(satwl, es_all(9,:), '-', 'Color', colors_bb9(7,:))
	plot(satwl, es_all(7,:), '-', 'Color', colors_bb9(8,:))
	plot(satwl, es_all(8,:), '-', 'Color', colors_bb9(9,:))

	set(gca, 'Box', 'on')
	set(gca, 'xlim', [350 800], 'ylim', [0 1])
	xlabel('Wavelength [nm]')
	ylabel('LED Relative Response')
%save LED_spectra es_all satwl
    
%load Wetlabs
%	plot(wl(1,:), es_all(1,:), '--', 'Color', colors_bb9(1,:))
%	plot(wl(2,:), es_all(2,:), '--', 'Color', colors_bb9(2,:))
%	plot(wl(3,:), es_all(3,:), '--', 'Color', colors_bb9(3,:))
%	plot(wl(4,:), es_all(4,:), '--', 'Color', colors_bb9(4,:))
%	plot(wl(5,:), es_all(5,:), '--', 'Color', colors_bb9(5,:))
%	plot(wl(6,:), es_all(6,:), '--', 'Color', colors_bb9(6,:))
%	plot(wl(7,:), es_all(7,:), '--', 'Color', colors_bb9(7,:))
%	plot(wl(8,:), es_all(8,:), '--', 'Color', colors_bb9(8,:))
%	plot(wl(9,:), es_all(9,:), '--', 'Color', colors_bb9(9,:))

load Wetlabs
	plot(wl(1,:), fl_all(1,:), ':', 'Color', colors_bb9(1,:))
	plot(wl(2,:), fl_all(2,:), ':', 'Color', colors_bb9(2,:))
	plot(wl(3,:), fl_all(3,:), ':', 'Color', colors_bb9(3,:))
	plot(wl(4,:), fl_all(4,:), ':', 'Color', colors_bb9(4,:))
	plot(wl(5,:), fl_all(5,:), ':', 'Color', colors_bb9(5,:))
	plot(wl(6,:), fl_all(6,:), ':', 'Color', colors_bb9(6,:))
	plot(wl(7,:), fl_all(7,:), ':', 'Color', colors_bb9(7,:))
	plot(wl(8,:), fl_all(8,:), ':', 'Color', colors_bb9(8,:))
	plot(wl(9,:), fl_all(9,:), ':', 'Color', colors_bb9(9,:))

plot(wl(1,:),photodiode(wl(1,:)),'--')    
print -djpeg90 -r300 ledcurves-bb9.jpg 


% ------ Plot for ECOVSFs
%figure, hold on
%set(gcf, 'paperunits', 'inches', 'papersize', [6.5 4.0])
%set(gcf, 'paperposition', [0.25 0.25 6.0 3.5])

	%plot(satwl, es_all(10,:), '-', 'Color', [0.0 0.0 1.0])
	%plot(satwl, es_all(11,:), '-', 'Color', [1.0 0.0 0.0])

	%set(gca, 'Box', 'on')
	%set(gca, 'xlim', [350 800], 'ylim', [0 1])
	%xlabel('Wavelength, [nm]')
	%ylabel('LED Relative Response')

%print -djpeg90 -r300 ledcurves-ecovsf.jpg 


