function norm_data = normbass(data, burst_per_avg, bass_error);
% function norm_data = normbass(data, burst_per_avg, bass_error);
%
% 	This function normalizes BASS burst accumulated sum data
%	into averages, correcting for BASS errors
%
%	data = accumulated sum in counts, axes A thru D
%			row=sample, col=axis
%	burst_per_avg = number of burst samples in sum, a vector
%	bass_error = number of errors received in each average, 
%			axes A thru D
%	norm_data = average in counts, axes A thru D

%	Written by Marinna Martini for the US Geological Survey
%	10/18/94	rewrite to be less memory intensive

if exist('burst_per_avg')~=1,
	% if the user didn't provide this info, make it up
	burst_per_avg=ones(length(data),1).*300;
end

% make sure the matrices are all in record x axis orientation
[r,c]=size(data);
if c > r, data=data'; end
[r,c]=size(bass_error);
if c > r, bass_error=bass_error'; end
[r,c]=size(data);

norm_data=zeros(r,c);
if exist('bass_error')~=1,
	% if the user didn't provide errors, assume there were none
	bass_error=zeros(r,c);
end

for i=1:c,
	norm_data(:,i) = data(:,i)./(burst_per_avg-bass_error(:,i));
end