#!/usr/bin/env bash
set -euo pipefail
export GMT_SESSION_NAME=$$
mkdir -p figures
gmt begin figures/GMT-18-001 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 makecpt '-C#00549f,white,#f6a800' -T-2/2/0.1 -Z -H > field.cpt
gmt grdview field.nc -R-4/4/-4/4/-2/2 -JX10c -JZ4c -Cfield.cpt -Qs '-N-2+g#e8f1fa' -p145/30 -Bxa2+lx -Bya2+ly -Bza1+lAmplitude -BWSneZ+tGMT-18-001
gmt end