#!/usr/bin/env bash

# 1. Create files needed in the loop
cat << 'EOF' > pre.sh
gmt begin pre
    gmt set FONT_ANNOT_PRIMARY=15p,Times-Roman,black
    gmt set FONT_LABEL=15p,Times-Roman,black
    gmt set MAP_GRID_CROSS_SIZE_PRIMARY=5p
    gmt set MAP_FRAME_PEN=thinnest,black
    gmt set MAP_TICK_LENGTH_PRIMARY=4p/2p
    gmt grd2cpt data/gauss_model.nc -Clapaz -R0/100/0/100 -Z -D
    gmt grdimage data/gauss_model.nc -R0/100/0/100 -Bxag+l"x (m)" -Byag+l"y (m)" -JX15c/15c -X4.5c -Y1.5c 
gmt end
EOF

# 2. Set up the main frame script
cat << 'EOF' > main.sh
gmt begin
#	Plot smooth blue curve and dark red dots at all steps so far
    gmt convert data/lgd_trace.txt -qi0:${MOVIE_FRAME} > data.txt
    gmt plot data.txt -W0.05p,white -R0/100/0/100 -JX15c/15c -X4.5c -Y1.5c -i1,2
    gmt plot data.txt -Sc0.05i -Gred -i1,2
    gmt plot data.txt -W0.05p,yellow -R0/100/0/100 -JX15c/15c -i3,4
    gmt plot data.txt -Sc0.05i -Gblack -i3,4
gmt end
EOF

# 3. Run the movie
gmt movie main.sh -Sbpre.sh -Cxga -Tdata/lgd_trace.txt -Vi -D5 -Zs -Nlgd_anim -Fmp4