[WIP] implement arithmetic mean symmetrization

problem: consistency check fails with relatively large errors (4% field
and ~16% velocity at particle positions).
But it may be due to I/O.

further checks are needed...