Add openmp parallel code fieldizing
Disable it because it is a bit slower
Tweak plots
Unify the column density functions and fix a unit bug
Deal with periodic boundary conditions