% AVGAQAABS - Original script from Peter Traykovski to average alot of ABS
% files
clear
jjind=0;
if 1   
   day1=datenum(1999,12,04,16,55,00);
   dayend=datenum(2000,4,20,0,00,00);
   
   dv=datevec(day1:.5/24:dayend);
   
   for jj=1:length(dv)
      fn=['..\aquadata\IrishLapDown\',[num2str(dv(jj,[3,2,1,4]),'%2.2d'),...
               num2str(round(dv(jj,5)+dv(jj,6)/60),'%2.2d')],'*.*'];
      %fns=['ss',num2str(dv(jj,2:4),'%2.2d'),'.**'];
      
      d=dir(fn);
      if size(d)
         jjind=jjind+1;
         D(jjind,:)=d.name;
      end
   end
   save datefile D dv
end
load datefile

for ii=1:length(D)
   %for ii=1:3
   disp([num2str(ii),'/',num2str(length(D))])

fname=['..\aquadata\IrishLapDown\', D(ii,:)];

[mdata1,mdata2,mdata3,sdatenum,freq,hdr]=readaqa(fname);
dt(ii,:)=sdatenum;
a1(ii,:)=mean(mdata1');
a2(ii,:)=mean(mdata2');
a3(ii,:)=mean(mdata3');

end
%yd=datenum(dv(:,1),dv(:,2),dv(:,3),dv(:,4),dv(:,5),dv(:,6));
%yd = dt;
%try decoding D into yd
yd=datenum(str2num(D(:,5:8)),str2num(D(:,3:4)),str2num(D(:,1:2)),str2num(D(:,9:10)),...
   str2num(D(:,11:12)),str2num(D(:,13:14)));
save aquaabsavg016 yd a1 a2 a3
figure(1);clg
rdep=1.25e6./([1:129]').^2;
%rdep=ones(size(rdep));
[m,ndata]=size(a1');
mrdep=dup(rdep,ndata);
day1_2000=datenum(2000,1,1)
day1_2000=0;

subplot(311)
imagesc(yd,1:129,log10(a1'./mrdep),[-1 3])
datetick('x',6)
%imagesc(log10(a1'./mrdep),[-1 3])
title('1.0 MHz ABS Tripod A')
subplot(312)
imagesc(yd-day1_2000,1:129,log10(a2'./mrdep),[-1 3])
%imagesc(log10(a2'./mrdep),[-1 3])
datetick('x',6)
title('2.5 MHz ABS Tripod A')

subplot(313)
imagesc(yd-day1_2000,1:129,log10(a3'./mrdep),[-1 3])
%imagesc(log10(a3'./mrdep),[-1 3])
datetick('x',6)
title('5.0 MHz ABS Tripod A')

colorbar('horiz')