global_ocean.128x64x15/diags_matlab/mit_meridflux.m

function T = mit_meridflux(surface_flux,dx,dy)
% function T = mit_meridflux(surface_flux,dx,dy)
  
  if nargin == 2
    area = dx;
  elseif nargin == 3
    area = dx.*dy;
  else
    error('need two or three input parameters')
  end
  [nx ny nt] = size(surface_flux);
  for kt=1:nt
    fdxdy(:,:,kt) = change(surface_flux(:,:,kt),'==',NaN,0).*area;
  end

  T = squeeze(sum(cumsum(fdxdy,2),1));
  
  return
1	function T = mit_meridflux(surface_flux,dx,dy)
2	% function T = mit_meridflux(surface_flux,dx,dy)
3
4	if nargin == 2
5	area = dx;
6	elseif nargin == 3
7	area = dx.*dy;
8	else
9	error('need two or three input parameters')
10	end
11	[nx ny nt] = size(surface_flux);
12	for kt=1:nt
13	fdxdy(:,:,kt) = change(surface_flux(:,:,kt),'==',NaN,0).*area;
14	end
15
16	T = squeeze(sum(cumsum(fdxdy,2),1));
17
18	return