<div dir="ltr">Hello,<div><br></div><div>I'm trying to reproduce the calculations for the overturning streamfunction described on <a href="https://ecco-v4-python-tutorial.readthedocs.io/ECCO_v4_Example_OSNAP.html">https://ecco-v4-python-tutorial.readthedocs.io/ECCO_v4_Example_OSNAP.html</a> 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: <br><div><br></div><div>lonPairs=[-44 -5;-56 -45];<br>latPairs=[60 56;51 60];<br>names={'OSNAP East','OSNAP West'};<br>gcmfaces_lines_transp(lonPairs,latPairs,names);<br>T=calc_transports(uTrans,vTrans,mygrid.LINES_MASKS,{'dh','dz','hfac'});</div><div>psi_moc=cumsum(T')/1e6;<br></div><div><br></div><div>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?</div></div><div><br></div><div>Thanks,</div><div><br></div><div>Samar</div><div><br></div></div>