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

Simple program to read Wilson lines generated using IP-Glasma