[jgs55]

Hi everyone! I am trying to implement a new discretization scheme for the convective term (div(phi,U)). So, in myCustomScheme.H I need to multiply the center values of the cells belonging to the stencil around the face f to certain predefined weights.

My question is how can I reconstruct every stencil, taking into account that the stencil is formed by Nx cells in the x-direction and Ny in the y-direction. I attach an image to clarify the problem (example: I need to reconstruct a 4x3 stencil --in red-- when I am working with the face f --in blue--).

I have found some OF methods used to reconstruct stencils in high-order schemes, but I do not know if they are applicable to my problem. Thanks in advance!

P.D: I am working with structured meshes.