Support cyclic grids for LAGRANTO interpolation
The problem is that contrary to the trilinear interpolation, the LAGRANTO interpolation first interpolates the index to access the fields at consecutive time steps. This float-type index was not allowed to go past the index in the grid, stopping computations to pass the cyclic grid boundary.
Closes #228 (closed)
Edited by Christoph Fischer