Added relaxation parameter THETA as well as CPU arrays for staggered grid velocities