[ecco-support] [EXTERNAL] Calculating transports with gcmfaces in Matlab

[Wang, Ou (US 329B)]

Samar,

What are the two variables uTrans and vTrans in T=calc_transports(uTrans,vTrans,mygrid.LINES_MASKS,{'dh','dz','hfac'});?
The names suggest they are transport, but they should be UVELMASS.*mygrid.mskW and VVELMASS.*mygrid.mskS.

You may want to replace T=calc_transports(uTrans,vTrans,mygrid.LINES_MASKS,{'dh','dz','hfac'}) with the following:
uTrans = UVELMASS.*mygrid.mskW;
vTrans = VVELMASS.*mygrid.mskS;
T=calc_transports(uTrans,vTrans,mygrid.LINES_MASKS,{'dh','dz'});

Best regards,
Ou Wang

From: ecco-support <ecco-support-bounces at mit.edu> on behalf of Samar Khatiwala <samkat6 at gmail.com>
Date: Monday, June 10, 2024 at 2:59 AM
To: ECCO support list, wider membership <ecco-support at mit.edu>
Subject: [EXTERNAL] [ecco-support] Calculating transports with gcmfaces in Matlab
Hello,

I'm trying to reproduce the calculations for the overturning streamfunction described on https://ecco-v4-python-tutorial.readthedocs.io/ECCO_v4_Example_OSNAP.html<https://urldefense.us/v3/__https:/ecco-v4-python-tutorial.readthedocs.io/ECCO_v4_Example_OSNAP.html__;!!PvBDto6Hs4WbVuu7!NGtdgDhELdz3kVixjJCZ0q9D_VYms1T7Hzzy8y1T4DYa0C8LDN7tENkuljBhQzluMQEaUo_MK2G514L6FD8$> with the Matlab version of gcmfaces. The OSNAP example described in the tutorial uses the python function calc_section_stf which returns 'psi_moc'. I can't quite figure out what the equivalent Matlab function is. The documentation refers to calc_overturn.m and calc_transports.m, but the former returns transports for all (179) quasi-zonal lines, while the output of calc_transports is nowhere close to what is plotted in the tutorial. After looking through the python code the closest I've gotten is with:

lonPairs=[-44 -5;-56 -45];
latPairs=[60 56;51 60];
names={'OSNAP East','OSNAP West'};
gcmfaces_lines_transp(lonPairs,latPairs,names);
T=calc_transports(uTrans,vTrans,mygrid.LINES_MASKS,{'dh','dz','hfac'});
psi_moc=cumsum(T')/1e6;

But I'm not confident this is correct. Can someone please tell me if this is right and if not what the correct way is?

Thanks,

Samar

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mit.edu/pipermail/ecco-support/attachments/20240610/a0f3ecb9/attachment-0001.htm>