#!/usr/bin/env bash
set -euo pipefail
export GMT_SESSION_NAME=$$
mkdir -p figures
gmt begin figures/GMT-13-006 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
gmt grdmath -R-8/30/30/48 -I1 X 12 SUB DUP MUL Y 40 SUB DUP MUL ADD -0.03 MUL EXP 1000 MUL = hill.nc
gmt grdimage hill.nc -JM13c -Cviridis -nc -E200 -Baf
gmt coast -Di -W0.3p,gray40
gmt grdcontour hill.nc -C100 -A200+f8p -S8 -W0.6p,white
gmt basemap -B+t"GMT-13-006"
gmt end