2 |
|
|
3 |
% [PsiB,ylat]=calcBolusPsiCube(d,g,GMform,blkFile); |
% [PsiB,ylat]=calcBolusPsiCube(d,g,GMform,blkFile); |
4 |
% |
% |
5 |
% Compute bolus streamfunction from GM scheme |
% Compute eddy-induced streamfunction from Gent and McWilliams scheme |
6 |
% |
% |
7 |
% Input arguments: |
% Input arguments: |
8 |
% The incoming field data (d) and grid data (g) must be in a structured |
% The incoming field data (d) and grid data (g) must be in a structured |
9 |
% array format (which is the format that comes from rdmnc): |
% array format (which is the format that comes from rdmnc): |
10 |
% d [Field data] Kwx,Kwy |
% d [Field data] Kwx, Kwy (Skew flux form) or GM_PsiX, GM_PsiY (advective form) |
11 |
% g [Grid data ] drF,rA,dxC,dyC,dxG,dyG,HFacW,HFacS |
% g [Grid data ] drF,rA,dxC,dyC,dxG,dyG,HFacW,HFacS |
12 |
% GMform [string] GM form 'Skew' or 'Advc' |
% GMform [string] GM form: 'Skew' or 'Advc' |
13 |
% blkFile [file name] Broken line file |
% blkFile [file name] Broken line file |
14 |
|
% mask [structure] Optional: Mask field for computation per basin, it assumes that |
15 |
|
% maskW and maskS are provided in a structure |
16 |
% Output arguments: |
% Output arguments: |
17 |
% PsiB : bolus streamfunction at interface level (in Sv) |
% PsiB : bolus streamfunction at interface level (in Sv) |
18 |
% ylat : meridional coordinate of PsiB |
% ylat : meridional coordinate of PsiB |