clear

load_llc2160
tic
    [lonPairs,latPairs,names]=gcmfaces_lines_pairs;
    gcmfaces_lines_transp(lonPairs,latPairs,names);
toc    

nx=2160; ny=nx*13; nz=90;
pp='/nobackupp11/dmenemen/DYAMOND/c1440_llc2160/mit_output/STATS/';

%year
fin=[pp 'U/U_year.bin'];
disp(fin)
UU=readbin(fin,[nx ny nz]);
fin=[pp 'V/V_year.bin'];
disp(fin)
VV=readbin(fin,[nx ny nz]);

%fldU=UVELMASS.*mygrid.mskW; fldV=VVELMASS.*mygrid.mskS;
%replaced by:
fldU=convert2gcmfaces(UU).*mygrid.hFacW; %UVEL==>UVELMASS
fldV=convert2gcmfaces(VV).*mygrid.hFacS;
fldU=fldU.*mygrid.mskW;
fldV=fldV.*mygrid.mskS;

%compute barotropic stream function:
tic
%   fldBAR=calc_barostream(fldU,fldV);
   fldBAR=calc_barostream(fldU,fldV,0);
toc

%compute transports along transects:
tic
   fldTRANSPORTS=1e-6*calc_transports(fldU,fldV,mygrid.LINES_MASKS,{'dh','dz'});
toc

bt_strmfn=convert2gcmfaces(fldBAR);
bt_strmfn=single(bt_strmfn);
save SECTION_year bt_strmfn fldTRANSPORTS

%next monthly:
bt_strmfn=nan(nx, ny, 12);
fldTRANSPORTS=nan(23, nz, 12);

yr=2020;
for mn=3:14 %
tic
	f_mon=datestr(datenum(yr,mn,1),'yyyymm');
	fin=[pp 'U/U_' f_mon '.bin'];
	disp(fin)
	UU=readbin(fin,[nx ny nz]);
	fin=[pp 'V/V_' f_mon '.bin'];
	disp(fin)
	VV=readbin(fin,[nx ny nz]);

fldU=convert2gcmfaces(UU).*mygrid.hFacW; %UVEL==>UVELMASS
fldV=convert2gcmfaces(VV).*mygrid.hFacS;
fldU=fldU.*mygrid.mskW;
fldV=fldV.*mygrid.mskS;

%compute barotropic stream function:
%   fldBAR=calc_barostream(fldU,fldV);
   fldBAR=calc_barostream(fldU,fldV,0);

%compute transports along transects:
   fldTRANSPORTS(:,:,mn-3+1)=1e-6*calc_transports(fldU,fldV,mygrid.LINES_MASKS,{'dh','dz'});

bt_strmfn(:,:,mn-3+1)=convert2gcmfaces(fldBAR);
save SECTION_month bt_strmfn fldTRANSPORTS -v7.3

toc
end

bt_strmfn=single(bt_strmfn);
save SECTION_month bt_strmfn fldTRANSPORTS -v7.3