Velocity field in a magnetohydrodynamics simulation illustrating the suppression of turbulence by shearing, data kindly provided by A. P. L. Newton (Sheffield University).

With the CSV input data converted to GMT grd files, we generate the arrows layer with

vfplot \
    --verbose \
    --width 7in \
    --margin 0.5/3/0 \
    --aspect 4 \
    --scale 5e-4 \
    --iterations 400/5 \
    --pen 0.1mm \
    --fill 255 \
    --length 0/1in \
    --cache 513 \
    --symbol triangle \
    --output vfplot.eps \
    grd/u.grd grd/v.grd

Background of divergence is calculated from the vector field using GMT's grdmath, composited with the arrows layer using GMT.

Download the plot as pdf or inspect the GitLab sources.