for i in range(10000):
    integrator.run(100)
    espressopp.tools.vmd.imd_positions(system, sock)