global_ocean.128x64x15/diags_matlab/mit_zonalmean.m

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

  [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	[nx ny nz] = size(field);
5
6	area = hfac.*repmat(dx,[1 1 nz]);
7	meanfield = squeeze(nanmean(field.*area)./nanmean(area));
8
9	knil = find(meanfield == 0);
10	meanfield(knil) = NaN;
11
12	return