set colors classic
  set encoding iso_8859_1

  set view equal xyz
  set ticslevel 0
  set palette model HSV functions 5./6.*gray, 1, 1
  set auto cb
  unset colorbox

  set xrange [0:20]
  set yrange [0:20]
  set zrange [0:1]
  set xlabel 'x'
  set ylabel 'y'
  set zlabel 'z'
  set border 1+2+4+8 + 16+32+64+128 + 256+512+1024+2048

  f = 2.
  splot \
    'velocitygrid.txt' u ($1-0.5):($2-0.5):($3-0.5):(0.3+0.5*log10($7)):(sqrt($4**2+$5**2+$6**2)) w p pt 7 ps variable lc palette, \
    'velocitygrid.txt' u ($1-0.5-$4*f):($2-0.5-$5*f):($3-0.5-$6*f):($4*2*f):($5*2*f):($6*2*f):(sqrt($4**2+$5**2+$6**2)) w vec lc palette

  pause -1 'Press Enter to continue: '