During discussions about unstructured grid support for Parcels v4, we realised that it is failrly trivial to combine the particle index positions (particle.xi, particle.yi, particle.zi and particle.ti) into one particle.grid_index(?), which can then be unraveled with modulo and division operators based on xdim, ydim, zdim and tdim.
The advantage would be that the memory footprint of the ParticleSet would be less, because it would save three Variables of N particles.