matlab-grid-generator/bin/zonalmean.m

function [Tm] = zonalmean(GRID,T,varargin);
% Tm = zonalmean(GRID,T,varargin);

T( find(isnan(T)) )=0;

if nargin==3
 dTm=(mean(T,4)-mean(varargin{1},4)).*GRID.mskc;
 mapmsk=GRID.mskc(:,:,1);
elseif nargin==4
 dTm=(mean(T,4)-mean(varargin{1},4)).*GRID.mskc;
 mapmsk=varargin{2};
else
 dTm=mean(T,4).*GRID.mskc;
 mapmsk=GRID.mskc(:,:,1);
end

Tm=squeeze(sum(dTm.*repmat(mapmsk,[1 1 size(T,3)]),1));
N=squeeze(sum( GRID.mskc.*repmat(mapmsk,[1 1 size(T,3)]),1));
Tm=Tm./sq(N);
1	function [Tm] = zonalmean(GRID,T,varargin);
2	% Tm = zonalmean(GRID,T,varargin);
3
4	T( find(isnan(T)) )=0;
5
6	if nargin==3
7	dTm=(mean(T,4)-mean(varargin{1},4)).*GRID.mskc;
8	mapmsk=GRID.mskc(:,:,1);
9	elseif nargin==4
10	dTm=(mean(T,4)-mean(varargin{1},4)).*GRID.mskc;
11	mapmsk=varargin{2};
12	else
13	dTm=mean(T,4).*GRID.mskc;
14	mapmsk=GRID.mskc(:,:,1);
15	end
16
17	Tm=squeeze(sum(dTm.*repmat(mapmsk,[1 1 size(T,3)]),1));
18	N=squeeze(sum( GRID.mskc.*repmat(mapmsk,[1 1 size(T,3)]),1));
19	Tm=Tm./sq(N);