The code is now capable of doing some BOMD simulation with moving lattice.
Results are not compatible with the GO results, as the ions are not treated
in reduced coordinates.