# basic test, where I put one particle in LB fluid, move LB fluid with ext_force and output the position after 500 steps
#
# parameter initiation
#

import espressomd

from espressomd import lb
import numpy as np

system              = espressomd.System()
system.time_step    = 0.1
system.cell_system.skin         = 0.2
system.thermostat.set_lb(kT=0)

system.box_l = [10, 10, 10]

system.part.add(pos=(5,5,5), type=1)

lbf = espressomd.lb.LBFluid(agrid = 1, dens = 1.0, visc = 1.0, tau = 0.1, fric = 1.0 , ext_force = [0.001, 0, 0])
system.actors.add(lbf)

i = 0
while i < 500:
    system.integrator.run(steps=1)
    i += 1
tmppos = system.part[0].pos
print i,tmppos