add function to poisson code to write output file which can be visualized with e.g. paraview

[wikfeldt]

this would be nice to have a visual view of what the poisson code does