Forcing term
test on density conservation periodic and walls, pass walls only if velocity has no vertical component
* Fix BC using numpy copy
* Correct rho calculation no sum funtion, but f0+f1+...
* Fix force calculation term F3 & F4
* Test script: forced fluid on periodic BC