improve output

* compute vert/horiz velocities
* add some notes in code (see code_dump.txt)