#!/usr/bin/env bash
set -euo pipefail
export GMT_SESSION_NAME=$$
mkdir -p figures
gmt begin figures/GMT-16-002 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
printf -- '-3.70 -3.70 -0.000103\n-1.90 -3.70 -0.003495\n-0.10 -3.70 -0.001275\n1.70 -3.70 -0.000005\n3.50 -3.70 -0.000000\n-3.70 -1.90 -0.015851\n-1.78 -1.74 -0.730415\n0.02 -1.90 -0.154352\n1.82 -1.74 0.006190\n3.50 -1.90 0.000012\n-3.70 -0.10 -0.026238\n-1.90 0.06 -0.749954\n-0.10 -0.10 0.089985\n1.70 0.06 1.008522\n3.50 -0.10 0.002694\n-3.70 1.70 -0.000465\n-1.78 1.86 -0.009437\n0.02 1.70 0.176760\n1.82 1.86 0.160590\n3.50 1.70 0.000915\n-3.70 3.50 -0.000000\n-1.90 3.50 -0.000003\n-0.10 3.50 0.000073\n1.70 3.50 0.000151\n3.50 3.50 0.000000\n' > samples.txt
gmt nearneighbor samples.txt -R-4/4/-4/4 -I0.1 -S1.2 -N4+m2 -Ggridded.nc
gmt makecpt '-C#00549f,white,#f6a800' -T-2/2/0.1 -Z -H > field.cpt
gmt grdimage gridded.nc -R-4/4/-4/4 -JX11c -Cfield.cpt -E200 -nc -Bxa2+lx -Bya2+ly -BWSne
gmt plot samples.txt -i0,1 -Sc0.12c '-G#132b36' -W0.3p,white
gmt colorbar -Cfield.cpt -DJBC+w8c/0.3c+h+o0/0.9c -Bxa1+lAmplitude
gmt basemap -B+t"GMT-16-002"
gmt end