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