| 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 |