Also new way of computing the lattice index from the coordinate.
Indices differ from previous version ones by 1 unit, no effect
in continuum limit/when averaging over configurations