Everything back to CPU, test with aligning global memory to 128 byte blocks. For now only 2nd order diffusion
Added advection + CFL calculation. To-do: write routine that does reductions on the domain itself, excluding ghost/padding cells. Overhead can get quite large now....
Re-added time integration
Added boundary conditions