Volume average in avg3d

Open mathrack opened this issue 4 years ago • 0 comments

Hello,

In the subroutine avg3d, when estimating the average of a 3D field, the cell at i=nx is skipped when nclx1=1 :

https://github.com/xcompact3d/Incompact3d/blob/e090f06e0d26809a3f2d730d223e31da903baff6/src/tools.f90#L841

It seems odd that skipping the cell at i=nx depends on the boundary condition at x=0 (nclx1) and not on the boundary condition at x=Lx (nclxn) ?

Oct 21 '21 11:10 mathrack