1. Implemented reduction
   - doesn't fully work on the gpu code [few errors in final_val.dat]
   - works fine on serial_mod
2. accelerate_flow() comes at the end of the kernel. One acc. before
   the timing loop and one deceleration after the timing loop.