Feature release

- Code reorganization
- Master switch for plot output
- Calculate velocity in operator form