function eco_clim=eco_climat(cdir)
% eco_climat - Create statistics and climatology
%
% Use as: eco_clim=eco_climat('\\tornado\ssdsdata\ecocimt\configuration\ecofl2.properties')
%
% Inputs: cdir properties directory
%
% Output: eco_clim = an array of statistics, mean, std, counts.
%         

% col(8) contains the calibrated data.  col(1) is the time
% Initialize the variables
%disp('Eco_clim');
% Get properties
%cdir = ['\\tornado\ssdsdata\ecocimt\configuration\ecofl2.properties'];
prop = readProperties(cdir);
otype='plot';
p=getprops(prop,otype);
% Write to log file
fid=fopen([p.ldir,'/',p.log],'at+');
fprintf(fid,'%s\n',['ECO ',p.ecotype,' device number: ', p.dev, ', serial number ', p.serial]);
fprintf(fid,'%s\n',['  1. ECO_CLIMAT: ' DATESTR(NOW)]);
fclose(fid);
% open calibrated data file which contains all current deployment data
disp(['Opening file: \\Tornado\ssdsdata\ecocimt\data\',p.ecotype,'.csv']);
input_array=load(['\\Tornado\ssdsdata\ecocimt\data\',p.ecotype,'.csv']);
carray=('rb'); %color array for plots
sdn=input_array(:,3);
pacific_time=input_array(:,14);
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);
minday=datestr(min(sdn),2);maxday=datestr(max(sdn),2);
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);   % hour
        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',[0 23],'Fontsize',10);
    xlabel('Hour (Pacific Time)');
    if j==1
        xlabel('Hour (GMT)');
        title(['Climatology: ',minday,' to ',maxday]);
        eco_clim=hr_array;
    end %if j==1
    ylab=[p.longname,' ',p.units];
    ylabel(ylab);
end %for j=1:2,
set(gcf,'paperposition',[.25,.5,5.5,6]);
eval(['print -djpeg ', p.fdir, '/ecoclim_',p.fname,'.jpg;']);
eval(['print -dpng ', p.fdir, '/ecoclim_',p.fname,'.png;']);

% Write to log file
fid=fopen([p.ldir,'/',p.log],'at+');
fprintf(fid,'%s\n',['  2. ECO_CLIMAT: Climatology hourly averages from ',minday,' to ',maxday,'.']);
fclose(fid);