matlab_class/gcmfaces_calc/calc_transports.m

function [FLD]=calc_transports(fldU,fldV,SECTIONS_MASKS);
%object:    compute transports through pre-defined sections
%inputs:    fldU and fldV are the fields of grid point transport
%           SECTIONS_MASKS is the set of sections along
%               which transports will integrated (SECTIONS_MASKS should 
%               have been produced by line_greatC_TUV_mask.m)
%output:    FLD is the array of transport profiles

global mygrid;

%initialize output:
n3=max(size(fldU.f1,3),1); n4=max(size(fldV.f1,4),1);
FLD=NaN*squeeze(zeros(length(SECTIONS_MASKS),n3,n4));

%prepare fldU/fldV:
fldU(isnan(fldU))=0; fldV(isnan(fldV))=0;

dxg=mk3D(mygrid.DXG,fldU); dyg=mk3D(mygrid.DYG,fldU); drf=mk3D(mygrid.DRF,fldU);
for k4=1:n4;
fldU(:,:,:,k4)=fldU(:,:,:,k4).*dyg.*drf*1e-6;
fldV(:,:,:,k4)=fldV(:,:,:,k4).*dxg.*drf*1e-6;
end;

%use array format to speed up computation below:
fldU=convert2array(fldU); fldV=convert2array(fldV);
n1=size(fldU,1); n2=size(fldU,2);
fldU=reshape(fldU,n1*n2,n3*n4); fldV=reshape(fldV,n1*n2,n3*n4); 

for iy=1:length(SECTIONS_MASKS); 

   %get list ofpoints that form a zonal band:
   mmu=convert2array(SECTIONS_MASKS(iy).mmuIn);
   nnu=find(~isnan(mmu)); mmu=mmu(nnu)*ones(1,n3*n4);
   mmv=convert2array(SECTIONS_MASKS(iy).mmvIn);
   nnv=find(~isnan(mmv)); mmv=mmv(nnv)*ones(1,n3*n4);

   %do the area weighed average along this band:
   tmpu=sum(fldU(nnu,:).*mmu,1);
   tmpv=sum(fldV(nnv,:).*mmv,1);

   %store:
   FLD(iy,:)=reshape(tmpu+tmpv,n3,n4);

end; 


1	gforget	1.1	function [FLD]=calc_transports(fldU,fldV,SECTIONS_MASKS);
2	gforget	1.2	%object: compute transports through pre-defined sections
3			%inputs: fldU and fldV are the fields of grid point transport
4			% SECTIONS_MASKS is the set of sections along
5			% which transports will integrated (SECTIONS_MASKS should
6			% have been produced by line_greatC_TUV_mask.m)
7			%output: FLD is the array of transport profiles
8	gforget	1.1
9			global mygrid;
10
11			%initialize output:
12			n3=max(size(fldU.f1,3),1); n4=max(size(fldV.f1,4),1);
13			FLD=NaN*squeeze(zeros(length(SECTIONS_MASKS),n3,n4));
14
15			%prepare fldU/fldV:
16			fldU(isnan(fldU))=0; fldV(isnan(fldV))=0;
17
18			dxg=mk3D(mygrid.DXG,fldU); dyg=mk3D(mygrid.DYG,fldU); drf=mk3D(mygrid.DRF,fldU);
19			for k4=1:n4;
20	gforget	1.3	fldU(:,:,:,k4)=fldU(:,:,:,k4).dyg.drf*1e-6;
21			fldV(:,:,:,k4)=fldV(:,:,:,k4).dxg.drf*1e-6;
22	gforget	1.1	end;
23
24			%use array format to speed up computation below:
25			fldU=convert2array(fldU); fldV=convert2array(fldV);
26			n1=size(fldU,1); n2=size(fldU,2);
27			fldU=reshape(fldU,n1n2,n3n4); fldV=reshape(fldV,n1n2,n3n4);
28
29			for iy=1:length(SECTIONS_MASKS);
30
31			%get list ofpoints that form a zonal band:
32			mmu=convert2array(SECTIONS_MASKS(iy).mmuIn);
33			nnu=find(~isnan(mmu)); mmu=mmu(nnu)ones(1,n3n4);
34			mmv=convert2array(SECTIONS_MASKS(iy).mmvIn);
35			nnv=find(~isnan(mmv)); mmv=mmv(nnv)ones(1,n3n4);
36
37			%do the area weighed average along this band:
38			tmpu=sum(fldU(nnu,:).*mmu,1);
39			tmpv=sum(fldV(nnv,:).*mmv,1);
40
41			%store:
42			FLD(iy,:)=reshape(tmpu+tmpv,n3,n4);
43
44			end;
45
46