checking results, do not print
allocate(c%displacements_P(1:3, 1:n_atoms, 0:2*(Lx/2+1)-1, 0:Ly-1, 0:c%local_Lz-1))