I’ve been working on a program, vfplot, which creates EPS files of curved arrows (representing a 2-d vector field). The program plots only the arrows, GMT can do the rest of the work needed to make an attractive plot.
Created with
#!/bin/sh
set -e
EXT='0.03125'
RNG="-R-$EXT/$EXT/-$EXT/$EXT"
DIM='4i/4i'
PRJ="-JX$DIM"
EPS='compose.eps'
# background scaler
GRD='/tmp/scalar-fine.grd'
grdsample grd/scalar.grd -G$GRD $RNG -I1024+n
grdimage $GRD $PRJ $RNG -Cscalar.cpt -K > $EPS
rm -f $GRD
# central cylinder
psxy $PRJ $RNG -Sc1.27i -G255 -W1p -K -O >> $EPS <<EOF
0 0
EOF
# vfplot output
psimage vfplot.eps -W$DIM -Dx0/0 -O >> $EPS
# results
psconvert -A -Tf $EPS
psconvert -A -E150 -Tg $EPS
# cleanup
rm -r $EPS
Other examples here, follow the GitLab links to inspect the GMT scripts compose.sh