integration does count redundant nodes twice?
tommy-engels opened this issue · comments
the integration of, e.g. mean flow, for instance
params_acm%mean_flow(1) = params_acm%mean_flow(1) + sum(u(g+1:Bs+g, g+1:Bs+g, 1, 1))*dx(1)*dx(2)
appears to be counting the redundant nodes twice; shouldn't the array boundary be
params_acm%mean_flow(1) = params_acm%mean_flow(1) + sum(u(g+1:Bs+g-1, g+1:Bs+g-1, 1, 1))*dx(1)*dx(2)
instead??
I remember having changed that once in the master branch (see 0a55c5d) which I pushed apparently exactly one day after you created the new physics branch and then did not change it there.