Change test function
* Use periodic BCs on all boundary
* Print solution
* Write solution to VTK-file
Change permeability field of test model
Also print solution (pressure and velocity) to VTK-file