tutorial_global_oce_latlon/diags_matlab/mit_zonalmean.m

function meanfield = mit_zonalmean(field,hfac,dx);
%function meanfield = mit_zonalmean(field,hfac,dx);

% $Header:  $
% $Name:  $

  [nx ny nz] = size(field);

  area = hfac.*repmat(dx,[1 1 nz]);
  meanfield = squeeze(nanmean(field.*area)./nanmean(area));
  
  knil = find(meanfield == 0);
  meanfield(knil) = NaN;
  
  return
1	function meanfield = mit_zonalmean(field,hfac,dx);
2	%function meanfield = mit_zonalmean(field,hfac,dx);
3
4	% $Header: $
5	% $Name: $
6
7	[nx ny nz] = size(field);
8
9	area = hfac.*repmat(dx,[1 1 nz]);
10	meanfield = squeeze(nanmean(field.*area)./nanmean(area));
11
12	knil = find(meanfield == 0);
13	meanfield(knil) = NaN;
14
15	return