global_ocean.128x64x15/diags_matlab/nansum.m

function y = nansum(x)
%NANSUM Sum ignoring NaNs.
%   NANSUM(X) returns the sum treating NaNs as missing values.  
%   For vectors, NANSUM(X) is the sum of the non-NaN elements in
%   X. For matrices, NANSUM(X) is a row vector containing the sum 
%   of the non-NaN elements in each column of X. 
%
%    See also NANMEDIAN, NANSTD, NANMIN, NANMAX, NANMEAN.

%   Copyright 1993-2000 The MathWorks, Inc. 
%   $Revision: 1.1 $  $Date: 2002/10/22 13:30:40 $

% Replace NaNs with zeros.
nans = isnan(x);
i = find(nans);
x(i) = zeros(size(i));

% Protect against an entire column of NaNs
y = sum(x);
i = find(all(nans));
y(i) = i + NaN;

1	cnh	1.1	function y = nansum(x)
2			%NANSUM Sum ignoring NaNs.
3			% NANSUM(X) returns the sum treating NaNs as missing values.
4			% For vectors, NANSUM(X) is the sum of the non-NaN elements in
5			% X. For matrices, NANSUM(X) is a row vector containing the sum
6			% of the non-NaN elements in each column of X.
7			%
8			% See also NANMEDIAN, NANSTD, NANMIN, NANMAX, NANMEAN.
9
10			% Copyright 1993-2000 The MathWorks, Inc.
11			% $Revision: 1.1 $ $Date: 2002/10/22 13:30:40 $
12
13			% Replace NaNs with zeros.
14			nans = isnan(x);
15			i = find(nans);
16			x(i) = zeros(size(i));
17
18			% Protect against an entire column of NaNs
19			y = sum(x);
20			i = find(all(nans));
21			y(i) = i + NaN;
22