% plots averaged spectra from hisdata_sim
% makes 4 plots to get idea of what data is doing
%
dum = input('Which index first?  ');

figure(1);
plot(freq,ones(index(dum),:),'r');
axis([890 910 0 150]);
title('Single Spectra');
xlabel('Wavenumber');
ylabel('Radiance');
hold;
plot(freq,ones(dum+1,:),'k');
plot(freq,ones(dum+2,:),'b');
plot(freq,ones(dum+3,:),'g');
plot(freq,ones(dum+4,:),'m');
plot(freq,ones(dum+5,:),'c');
hold;

figure(2);
plot(freq,twos((dum*5)+1,:),'r');
axis([890 910 0 150]);
title('Two Spectra Averages');
xlabel('Wavenumber');
ylabel('Radiance');
hold;
plot(freq, twos(dum*5+2,:),'k');
plot(freq, twos(dum*5+3,:),'b');
plot(freq, twos(dum*5+4,:),'g');
plot(freq, twos(dum*5+5,:),'m');
hold;

figure(3);
plot(freq,threes(dum*4+1,:),'r');
axis([890 910 0 150]);
title('Three Spectra Averages');
xlabel('Wavenumber');
ylabel('Radiance');
hold;
plot(freq,threes(dum*4+2,:),'k');
plot(freq,threes(dum*4+3,:),'b');
plot(freq,threes(dum*4+4,:),'g');
hold;

figure(4);
plot(freq,fours(dum*3+1,:),'r');
axis([890 910 0 150]);
title('Four, Five, and Six Spectra Averages');
xlabel('Wavenumber');
ylabel('Radiance');
hold;
plot(freq,fours(dum*3+2,:),'k');
plot(freq,fours(dum*3+3,:),'b');
plot(freq,fives(dum*2+1,:),'g');
plot(freq,fives(dum*2+2,:),'m');
plot(freq,sixes(dum,:),'c');
hold;

% this prints all the single spectra and averages (2,3,4,5,6) in different
% frames