function plot_ecofl(ecofl_input_array,ecobb_input_array,bbdev_num,fldev_num,p,days)
% Plot calibrated ECO data.  Multi-panel figures of calibrated ecofl and
% ecobb
% Usage:
%       a = plot_eco(ecobb_input_array,ecofl_input_array,bbdev_num,fldev_num,cdir)
%        inputarray : data input array containing serial date, ref light, ref, thermistor, siam timestamp, oasis timestamp
%        dev_num=device number, cdir=conf dir
%        ecofl : WetLabs ECO fluorometer, ecobb : WetLabs ECO scattering
%        meter

% rpm MBARI 2004
% assign variables and calculate interim values
%figure('visible','off'); 
input_array=ecofl_input_array;
disp('Plotting ECO Devices ...');
ecoflcal=input_array(:,8); %calibrated fluorescence
ecoflraw=input_array(:,2); % averaged raw fluorescence counts
ecoflcaltemp=input_array(:,9); %temperature corrected calibrated fluorescence
ecofltemp=input_array(:,10); % temperature corrected value
ecofltempraw=input_array(:,5); % averaged raw thermistor counts
ecoflrefraw=input_array(:,7); % averaged raw reference light counts
sdn=input_array(:,1); % serial date number timestamp
ecofl_avg=num2str(mean(ecoflcal)); ecofl_std=num2str(std(ecoflcal)); ecoflmin=num2str(min(ecoflcal));ecoflmax=num2str(max(ecoflcal));
[row col]=size(input_array);
% Plot the data; a five panel figure with calibrated eco data and derived
% values

tlabel=['MBARI ',p.mooring,'/',p.project,' Mooring Fluorescence (',num2str(fldev_num), ')'];
last_day=datenum(sdn(row,:));
first_day=round(last_day-days)+0.1;
XLim = [first_day last_day];
interval=round(days/7);
XTick = [first_day:interval:last_day];
XTickLabel=datestr([first_day:interval:last_day],6);
%XTick = [datenum(sdn(1,:)):datenum(sdn(row,:))];
Height = 00.225; Width=0.775; n=3;

figNumber=figure(...
   'Name',['MBARI ',p.project,' ',p.mooring, ' : ECO'],...
	'Units','inches',...
	'Position',[.25 -0.00 6.5 6.5],...
    'PaperOrientation','portrait',...
    'Visible', 'OFF',...
	'PaperPosition',[0.25 .5 5 8]);
Position = [0.10 (0.2 + (n-1) * Height) Width Height];
subplot('position',Position)
figure; close
plot(sdn,ecoflcal,'g.'); hold; plot(sdn,ecoflcal,'k--');
ylabel('Chl (mg m^-^3)');
title(tlabel);
ymax=ceil(max(ecoflcal(:)));
set(gca,'XGrid','on','XTick',XTick,'XLim',XLim,'XTickLabel',XTickLabel,'Fontsize',10);
xlabel
%pause(10);
% ECOBB
input_array=ecobb_input_array;
sdn=input_array(:,1); % serial date number timestamp
ecobbraw=input_array(:,2); % averaged raw backscatter counts
ecobbtempraw=input_array(:,5); % averaged raw thermistor counts
ecobbrefraw=input_array(:,7); % averaged raw reference light counts
ecobbcal=input_array(:,8); %calibrated backscatter
ecobbvol=input_array(:,9); %volume scattering coefficients
ecobbpart=input_array(:,10); % volume scattering of particulates only
ecobb_avg=num2str(mean(ecobbcal)); ecobb_std=num2str(std(ecobbcal)); ecobbmin=num2str(min(ecobbcal));ecobbmax=num2str(max(ecobbcal));
[row,col] = size(input_array); % returns the size of matrix X in separate variables.

% Plot the data; a two panel figure with calibrated data and raw data
tlabel=['MBARI ',p.mooring,'/',p.project,' Mooring Scattering Meter (',num2str(bbdev_num), ')'];
last_day=datenum(sdn(row,:));
first_day=round(last_day-days)+0.1;
XLim = [first_day last_day];
interval=round(days/7);
XTick = [first_day:interval:last_day];
XTickLabel=datestr([first_day:interval:last_day],6);

n=n-2;Position = [0.10 (0.3 + (n-1) * Height) Width Height];
subplot('position',Position)
figure; close
plot(sdn,ecobbcal,'b.'); hold; plot(sdn,ecobbcal,'k--');
ylabel('B_b(\lambda) (m^-^1)');
title(tlabel);
set(gca,'XGrid','on','XTick',XTick,'XLim',XLim,'XTickLabel',XTickLabel,'Fontsize',10);

% Display Calibration File
h = axes('Position',[0 0 1 1],'Visible','off');
set(gcf,'CurrentAxes',h);
if and(str2num(datestr(now,5))>2,str2num(datestr(now,5))<10)
    datestring=[datestr(now), ' PDT'];
else
    datestring=[datestr(now), ' PST'];
end %if datestr
h2=text(.7,.07,datestring,'FontSize',6);
h3=text(.125,.23,'B_b(\lambda)=Total backscattering coefficients at 660. ','FontSize',8);
h2=text(.05,.125,['Where Project=',p.project,', Angle=117, Wavelength=660, and time is GMT.'],'FontSize',7);
h2=text(.05,.08,['Statistics:', datestr(first_day),' to ', datestr(last_day)],'FontSize',7);
h3=text(.05,.06,['  ecofl avgs:',ecofl_avg,' ecofl std:',ecofl_std,' ecofl min:',ecoflmin,' ecofl max:',ecoflmax],'FontSize',7);
h3=text(.05,.04,['  ecobb avgs:',ecobb_avg,' ecobb std:',ecobb_std,' ecobb min:',ecobbmin,' ecobb max:',ecobbmax],'FontSize',7);

set(gcf,'paperposition',[.25,.5,4.5,8]);
%eval(['print -djpeg ', p.fdir, '/eco',num2str(days),'.jpg;']);
%eval(['print -dpng ', p.fdir, '/eco',num2str(days),'.png;']);
print('-dpng', '-painters', '-f1', [p.fdir '/eco' num2str(days) '.png']);
print('-djpeg95', '-painters','-f1', [p.fdir '/eco' num2str(days) '.jpg']);

% Write to log file
fid=fopen([p.ldir,'/',p.log],'at+');
fprintf(fid,'%s\n',['     Statistics: ecofl avgs:',ecofl_avg,' ecofl std:',ecofl_std,' ecofl min:',ecoflmin,' ecofl max:',ecoflmax]);
fprintf(fid,'%s\n',['                 ecobb avgs:',ecobb_avg,' ecobb std:',ecobb_std,' ecobb min:',ecobbmin,' ecobb max:',ecobbmax]);
fprintf(fid,'%s\n',['  5. PLOT_ECO: plotting from ', datestr(first_day),' to ', datestr(last_day),' GMT']);
fprintf(fid,'%s\n',['     plot file: \\Tornado\ssdsdata\ecocimt\figures\eco',num2str(days),'.jpg']);
fprintf(fid,'%s\n',':::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::');
fprintf(fid,'\n');
fclose(fid);

disp(['Statistics: From ', datestr(first_day),' to ', datestr(last_day)]);
disp(['  ecofl avgs:',ecofl_avg,' ecofl std:',ecofl_std,' ecofl min:',ecoflmin,' ecofl max:',ecoflmax]);
disp(['  ecobb avgs:',ecobb_avg,' ecobb std:',ecobb_std,' ecobb min:',ecobbmin,' ecobb max:',ecobbmax]);