#!/usr/bin/env bash
set -euo pipefail
export GMT_SESSION_NAME=$$
mkdir -p figures
gmt begin figures/GMT-17-005 png,pdf E200,H2,Qg4,Qt4 -C
gmt set GMT_AUTO_DOWNLOAD off MAP_FRAME_TYPE plain PROJ_LENGTH_UNIT cm FONT_ANNOT_PRIMARY 9p FONT_LABEL 10p FONT_TITLE 11p PS_LINE_CAP round PS_LINE_JOIN round COLOR_NAN gray85
gmt grdmath -R-4/4/-4/4 -I0.1 X 1 SUB SQR Y 0.5 SUB SQR ADD NEG EXP 2 MUL X 1.4 ADD SQR Y 0.8 ADD SQR ADD -0.7 MUL EXP 1.5 MUL SUB = field.nc
gmt project -C-4/0.5 -E4/0.5 -G0.05 -N > track.txt
gmt grdtrack track.txt -Gfield.nc > sampled.txt
gmt plot sampled.txt -i0,3 -R-4/4/-2/2 -JX11c/6c '-W1.5p,#f6a800' -Bxa2+lx -Bya1+lAmplitude -BWSne
gmt basemap -B+t"GMT-17-005"
gmt end