function m = meanRefl(data)

% Calculates the mean of a 180x288 reflectivity matrix and plots the
% result.
% Input: data - 180x288 reflectivity matrix
% Output: m - 180x1 vector of longitudal mean values for each latitude

% Initialize sizes
m = zeros(180,1);

for i=1:180
  lonSum = 0;
  count = 0;
  % Summarize and count values which are not NaN for each latitude
  for j=1:288
    if ~isnan(data(i,j))
      lonSum = lonSum + data(i,j);
      count = count+1;
    end
    % Loop next longitude
  end
  % Calculate mean if any values
  if count > 0
    m(i) = lonSum/count;
  else
    m(i) = NaN;
  end
  % Loop next latitude
end

% Plot the result
plot(-89.5:89.5,m);
axis([-89.5 89.5 0 100]);
title('Mean Reflectivity');
xlabel('Latitude [deg]');
ylabel('Reflectivity [%]');