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