function b = read_ts ( a,b );
% ADV_HDR - Assign readable names to ADV .ts records
% ts = a ( a, b );

% Chris Sherwood, USGS
%last modified 2/6/01
% J. Lacy

[nr,nc]=size(a);
nprof = nr;     % total number of profiles
b.bnum  = a(:,1);  % burst number
b.num_in_bur = a(:,2);  %sample number within burst
if b.coord_sys == 'XYZ'   
    b.vx = a(:,3);  %Vx or Veast, cm/s
    b.vy = a(:,4);  %Vy or Vnorth, cm/s
 end
 if b.coord_sys =='ENU'
     b.ve = a(:,3);  %Vx or Veast, cm/s
     b.vn = a(:,4);  %Vy or Vnorth, cm/s
end
b.vz = a(:,5);  %Vz or Vup, cm/s
b.amp1 = a(:,6);  %signal strength, receiver 1, counts
b.amp2 = a(:,7);  %signal strength, receiver 2, counts
b.amp3 = a(:,8);  % signal strength, receiver 3, counts
b.co1 = a(:,9);  %correlation, receiver 1, per cent
b.co2 = a(:,10);  % correlation, receiver 2, per cent
b.co3 = a(:,11);  % correlation, receiver 3, per cent
b.hdg  = a(:,12); % heading (degrees magnetic)
b.pitch= a(:,13); % rotation about Y axis
b.roll = a(:,14); % rotation about X axis
b.t    = a(:,15); % temperature, deg C 
%b.p    = a(:,16); % mean pressure, dBars
if nc > 16
b.obs1 = a(:,17); %voltage from obs1
b.obs2 = a(:,18); %voltage from obs2
if nc > 18
b.p = a(:,19);  %pressure from ext pressure sensor, dbar
end
end