tutorial_global_oce_latlon/diags_matlab/mit_meridflux.m

function T = mit_meridflux(surface_flux,dx,dy)
% function T = mit_meridflux(surface_flux,dx,dy)
  
% $Header:  $
% $Name:  $

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