function [mdata1,mdata2,mdata3,sdatenum,freq,hdr]=readaqa(fname);
% READAQA - Read Aquatec ABS data file
% [mdata1,mdata2,mdata3,sdatenum,freq,hdr]=readaqa(fname);

% Written by Peter Traykovski, WHOI
% Modified by Amit Bohara (this is his v1)
% Comments added by Chris Sherwood
% Last revised Sept. 5, 2003
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);
% next three lines were changed to 128 from 129 by Amit Bohara
mdata1=mdata(1:128,:);
mdata2=mdata(128+(1:128),:);
mdata3=mdata(2*128+(1:128),:);
%The remaining last three rows does contain some data.
%I am not sure what it is, but it doesnt show up in the
%aquatec data set.