Get an iterator over the space-time coordinates of the field.
This function returns a matrix of space-time coordinates of the field. Both for the whole field as well as the truncated field (without the margin)
dimension of the original field (first dimension is time; rest is space)
template of the LC coordinates
1 2 3