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	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