function [mdata1,mdata2,mdata3,sdatenum,freq,hdr]=readaqa_old(fname);
% READAQA_OLD - Version I used to use for processing before Sept. 2003 - CRS
fid = fopen(fname,'r');
command=fread(fid,1,'uint8');
junk=fread(fid,1,'uint8');
rec_size=fread(fid,1,'uint16');
chksum=fread(fid,1,'uint16');
sdata=fread(fid,rec_size/2,'uint16');


command=fread(fid,1,'uint8');
junk=fread(fid,1,'uint8');
rec_size=fread(fid,1,'uint16');
chksum=fread(fid,1,'uint16');
sdata=fread(fid,rec_size,'uchar');

hdr=setstr(sdata(3:70)');
year=str2num([dec2hex(sdata(82)) dec2hex(sdata(81) )]);
month=str2num([dec2hex(sdata(79))]);
day=str2num([dec2hex(sdata(80))]);
hour=str2num([dec2hex(sdata(78))]);
minute=str2num([dec2hex(sdata(77))]);
second=str2num([dec2hex(sdata(76))]);
sdatenum=datenum(year,month,day,hour,minute,second);

command=fread(fid,1,'uint8');
junk=fread(fid,1,'uint8');
rec_size=fread(fid,1,'uint16');
chksum=fread(fid,1,'uint16');
sdata=fread(fid,rec_size/2,'uint16')';
samp_rate=sdata(1);avg_over=sdata(2);
freq=samp_rate./avg_over;%in hz
nbins=sdata(11);

command=fread(fid,1,'uint8');
junk=fread(fid,1,'uint8');
rec_size=fread(fid,1,'uint16');
chksum=fread(fid,1,'uint16');
sdata=fread(fid,'uint16')';
fclose(fid);

ndata=ceil(length(sdata)/387);
%extr_pts=length(sdata)-387*ndata;
%fseek(fid,-384*2,1)
zdata=zeros(1,ndata*387);
zdata(1:length(sdata))=sdata;
mdata=reshape(zdata,387,ndata);
mdata1=mdata(1:129,:);
mdata2=mdata(129+(1:129),:);
mdata3=mdata(2*129+(1:129),:);