In order to incorporate a new term, I am dsicretizing a fourth order term for which I need to access nodes which are second adjacent ( for example, north of north node) and so on.
Like done for north node as, i_left = equation_stencil%c(-1,0,0) , I tried to make similar changes and declared a variable i_left_left=equation_stencil%c(-2,0,0) and so on.
There was no problem in compilation.
However, while running the code, there was an error:
Subscript array C has value -2 which is less than lower bound of -1.
I searched in the files and found that a similar declaration has already been used ( for example, variable i_front2 and so on).
May you please help in this regard.
( The changes are being made in the file " discretize_cell_transport_equation.f90"