Hi all,
As Ryo suggested there are easier way to calculate volume transport through a wall.
I.e. if you want to calculate volume transport across 10N, between 60:90E and 0:300 m, then it may be represented as below:
Vol_tr=v[x=60:90@din,y=10,z=0:300@din]
Thank you. I had something like that in mind.
It really depends on your purpose and the accuracy you pursue, but I just put some random thoughts/caveats here:
- When your grid is non-uniform, @DIN will result in some error, if the axis isn't "properly" set. List the XBOX, YBOX, or ZBOX variable to check. (If we are dealing with a netCDF file, I'm referring to the "_bnds" information.)
- If your model uses a C-grid, then @DIN is the right thing. If your model uses a B-grid, you will first need to map the velocity onto the C-grid points, like
let vm = v[gx=C_grid_xaxis@LIN]
- If you really need volume fluxes through each face of each gridbox, they can be calculated as
let vflx = v * xbox[g=v] * zbox[g=v]
provided that your non-uniform axes are "properly" set.
Hope this helps,
Ryo