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	cnh	1.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);