/usr/share/doc/gmt/examples/anim01/anim_01.sh is in gmt-examples 5.2.1+dfsg-3build1.
This file is owned by root:root, with mode 0o755.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 | #!/bin/bash
# GMT ANIMATION 01
# $Id: anim_01.sh 15178 2015-11-06 10:45:03Z fwobbe $
#
# Purpose: Make web page with simple animated GIF of sine function
# GMT progs: gmt gmtset, gmt gmtmath, gmt psbasemap, gmt pstext, gmt psxy, gmt psconvert
# Unix progs: printf, mkdir, rm, mv, echo, convert, cat
# Note: Run with any argument to build movie; otherwise 1st frame is plotted only.
#
# 1. Initialization
# 1a) Assign movie parameters
. gmt_shell_functions.sh
width=4i
height=2i
dpi=125
n_frames=18
name=anim_01
ps=${name}.ps
# 1b) Do frame-independent calculations and setup
angle_step=`gmt gmtmath -Q 360 ${n_frames} DIV =`
angle_inc=`gmt gmtmath -Q ${angle_step} 10 DIV =`
gmt psbasemap -R0/360/-1.2/1.6 -JX3.5i/1.65i -P -K -X0.35i -Y0.25i \
-BWSne+glightgreen -Bxa90g90f30+u\\312 -Bya0.5f0.1g1 \
--PS_MEDIA=${width}x${height} --FONT_ANNOT_PRIMARY=9p > $$.map.ps
# 2. Main frame loop
mkdir -p $$
frame=0
while [ ${frame} -le ${n_frames} ]; do
# Create file name using a name_##.tif format
file=`gmt_set_framename ${name} ${frame}`
cp -f $$.map.ps $$.ps
angle=`gmt gmtmath -Q ${frame} ${angle_step} MUL =`
if [ ${frame} -gt 0 ]; then # First plot has no curves
# Plot smooth blue curve and dark red dots at all angle steps so far
gmt gmtmath -T0/${angle}/${angle_inc} T SIND = $$.sin.d
gmt psxy -R -J -O -K -W1p,blue $$.sin.d >> $$.ps
gmt gmtmath -T0/${angle}/${angle_step} T SIND = $$.sin.d
gmt psxy -R -J -O -K -Sc0.1i -Gdarkred $$.sin.d >> $$.ps
fi
# Plot red dot at current angle and annotate
sin=`gmt gmtmath -Q ${angle} SIND =`
gmt psxy -R -J -O -K -Sc0.1i -Gred >> $$.ps <<< "${angle} ${sin}"
printf "0 1.6 a = %03d" ${angle} | gmt pstext -R -J -F+f14p,Helvetica-Bold+jTL -O -K \
-N -Dj0.1i/0.05i >> $$.ps
gmt psxy -R -J -O -T >> $$.ps
[[ ${frame} -eq 0 ]] && cp $$.ps ${ps}
if [ $# -eq 0 ]; then
gmt_cleanup .gmt
gmt_abort "${0}: First frame plotted to ${name}.ps"
fi
# RIP to TIFF at specified dpi
gmt psconvert -E${dpi} -Tt $$.ps
mv -f $$.tif $$/${file}.tif
echo "Frame ${file} completed"
frame=`gmt_set_framenext ${frame}`
done
# 3. Create animated GIF file and HTML for web page
${GRAPHICSMAGICK-gm} convert -delay 20 -loop 0 $$/${name}_*.tif ${name}.gif
cat << END > ${name}.html
<HTML>
<TITLE>GMT Trigonometry: The sine movie</TITLE>
<BODY bgcolor="#ffffff">
<CENTER>
<H1>GMT Trigonometry: The sine movie</H1>
<IMG src="${name}.gif">
</CENTER>
<HR>
We demonstrate how the sine function <I>y = sin(a)</I> varies with <I>a</I> over
the full 360-degree interval. We plot a bright red circle at each
new angle, letting previous circles turn dark red. The underlying
sine curve is sampled at 10 times the frame sampling rate in order to reproduce
a smooth curve. Our animation uses GraphicsMagick's convert tool to make an animated GIF file
with a 0.2 second pause between frames, set to repeat forever.
<HR>
<I>${name}.sh: Created by ${USER} on `date`</I>
</BODY>
</HTML>
END
# 4. Clean up temporary files
gmt_cleanup .gmt
|