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