- step class for LB phasefield generic enough to work with 3-phase and
  N-phase models
- cahn hilliard can either be solved by LBM or by finite differences
- 3 phase model can be solved with rho phase or without

- 3 phase model works now with step sytem

- plotting
- calculation of force on boundary
- vtk outpu

- Step working for serial CPU scenarios

- scaling interface width eta instead of surface tensions tau to correct
  interface profile & surface tensions

because boundary handling does its own flag management