function eco_clim=eco_climat2(input_array,p,days)
% eco_climat - Create statistics and climatology
%
% Use as: [ecobb_stats,ecobb_clim]=eco_climat(ecobb_output_array,p,7)
%
% Inputs: input_array = Array of calibrated data and configuration
% directory as a string. p are properties an output from getprop, and days are the number of days
% processed.
%
% Output: eco_stats = an array of statistics, mean, std, counts.
%         eco_clim = an array of climatological values

% col(8) contains the calibrated data.  col(1) is the time
% Initialize the variables
%disp('Eco_clim');
carray=('rb');
sdn=input_array(:,1);
pacific_time=input_array(:,11);
ecodata=input_array(:,8);
jday=date2jul(sdn);
year=datestr(sdn,10);
date_time=datestr(sdn,0);
hours=datestr(sdn,15);
hour=(str2num(hours(:,1))*10)+str2num(hours(:,2));
minutes=(str2num(hours(:,4))*10)+str2num(hours(:,5));
jtime=jday+((hour/24) + ((minutes/60)/24));
pacific_date = sdn-(7/24);
for j=1:2,
    if j==2
        hours=datestr(pacific_time,15);
        hour=(str2num(hours(:,1))*10)+str2num(hours(:,2));
    end %if j==2
    for i=1:24,
        eval(['hr_ndx=find(hour==',num2str(i-1),');' ]);
        hr_array(i,1)=(i-1);
        hr_array(i,2)=mean(ecodata(hr_ndx));
        hr_array(i,3)=std(ecodata(hr_ndx));
        hr_array(i,4)=length(ecodata(hr_ndx));
        hr_array(i,5)=median(ecodata(hr_ndx));
        hr_array(i,6)=min(ecodata(hr_ndx));
        hr_array(i,7)=max(ecodata(hr_ndx));
        hr_array(i,8)=cov(ecodata(hr_ndx));
    end %for i=1:24
    subplot(2,1,j)
    eval(['plot(hr_array(:,1),hr_array(:,2),''',carray(:,j),'.-'');']);
    %hold;
    %errorbar(hr_array(:,1),hr_array(:,2),hr_array(:,3));
    set(gca,'XGrid','on','XTick',[1,6,12,18,24],'XLim',[1 23],'Fontsize',10);
    xlabel('Hour (Pacific Time)');
    if j==1
        xlabel('Hour (GMT)');
        title([num2str(days),' day climatology']);
    end %if j==2
    eco_clim=hr_array;
    ylab=[p.longname,' ',p.units];
    ylabel(ylab);
end %or j=1:2,
set(gcf,'paperposition',[.25,.5,5.5,6]);
eval(['print -djpeg ', p.fdir, '/ecoclim_',p.fname,num2str(days),'.jpg;']);
eval(['print -dpng ', p.fdir, '/ecoclim_',p.fname,num2str(days),'.png;']);
% Write to log file
fid=fopen([p.ldir,'/',p.log],'at+');
fprintf(fid,'%s\n',['  4. ECO_CLIMAT: Climatology hourly averages.']);
fclose(fid);