function flag = bbopkq(data, depth, dww, max1stDiff, maxStd)
% BBOPKQ    - Quantify quality of incident illumination (Es)
%
% Use as: flag = bbopkq(data, depth, dww, max1stDiff, maxStd)
%
% Inputs: data  = prr data in engineering units (preferably only give downcast data) (vector only)
%         depth = a vector of depths corresponding to data
%         dww   = Depth window width (in meters)
%         m1    = Max mean 1st difference
%         mStd  = Max standard deviation
%
% Output: flag = vector of qc flags, 1 = good, 0 = bad

% Brian Schlining
% 27 Jan 2000

[r c] = size(data);
flag  = ones(r,1); % 1 is a good, 0 is bad

warning off
for i = 1:r
   top = depth(i) - dww/2;
   bottom = depth(i) + dww/2;
   good = find(depth > top & depth <= bottom);
   x = data(good);
   mean1stDiff = sum(abs(diff(x)))/(length(good)-1);
   xStd = std(x);
   if mean1stDiff >= max1stDiff
      flag(i) = 0;
   end
   if xStd >= maxStd
      flag(i) = 0;
   end
   %fprintf(1, 'n = %04i, 1stDiff = %10.6f, std = %10.6f, Depth = %7.2f, FLAG = %1i\n', ...
   %   i, mean1stDiff, xStd, depth(i), flag);
end
warning on