global_ocean.128x64x15/diags_matlab/mit_globalmean.m

function [mdt, mdtt, mdtstd] = mit_globalmean(dt,rac3d)
% function [mdt, mdtt, mdtstd] = mit_globalmean(dt,rac3d)
  
  denom = nansum(nansum(rac3d));
  mdt    = squeeze(nansum(nansum(dt.*rac3d))./denom);
  mdtt   = squeeze(nansum(nansum(abs(dt).*rac3d))./denom);
  mdtstd = sqrt(squeeze(nansum(nansum(dt.^2.*rac3d))./denom));

  return
1	cnh	1.1	function [mdt, mdtt, mdtstd] = mit_globalmean(dt,rac3d)
2			% function [mdt, mdtt, mdtstd] = mit_globalmean(dt,rac3d)
3
4			denom = nansum(nansum(rac3d));
5			mdt = squeeze(nansum(nansum(dt.*rac3d))./denom);
6			mdtt = squeeze(nansum(nansum(abs(dt).*rac3d))./denom);
7			mdtstd = sqrt(squeeze(nansum(nansum(dt.^2.*rac3d))./denom));
8
9			return