bug fixed, simulation working. TODO boundary conditions, density field