#!/usr/bin/env bash
set -euo pipefail
export GMT_USERDIR="$PWD/.gmt"
mkdir -p "$GMT_USERDIR"
gmt begin VIZ-10-005 pdf
gmt set \
  FONT_ANNOT_PRIMARY 7p \
  FONT_LABEL 8p \
  FONT_TITLE 8p \
  MAP_TITLE_OFFSET 4p \
  MAP_FRAME_TYPE plain \
  MAP_TICK_LENGTH_PRIMARY 3p
gmt plot -R0/420/0/280 -JX420p/280p -Gwhite -W0p <<'CANVAS'
0 0
420 0
420 280
0 280
CANVAS
gmt text -F+f8p+jTR <<'ID'
414 272 VIZ-10-005
ID
gmt text -F+f7p+jML <<'LEGEND'
55 246 Survived: 0% (light) to 100% (blue); grey: n=0
LEGEND
echo "58 234" | gmt plot -Sr6p/6p -G#e2edf8
echo "64 234" | gmt plot -Sr6p/6p -G#d7e5f3
echo "70 234" | gmt plot -Sr6p/6p -G#cbddef
echo "76 234" | gmt plot -Sr6p/6p -G#bfd6ea
echo "82 234" | gmt plot -Sr6p/6p -G#b4cee6
echo "88 234" | gmt plot -Sr6p/6p -G#a8c6e1
echo "94 234" | gmt plot -Sr6p/6p -G#9dbedc
echo "100 234" | gmt plot -Sr6p/6p -G#91b6d8
echo "106 234" | gmt plot -Sr6p/6p -G#85aed3
echo "112 234" | gmt plot -Sr6p/6p -G#7aa6cf
echo "118 234" | gmt plot -Sr6p/6p -G#6e9fca
echo "124 234" | gmt plot -Sr6p/6p -G#6397c6
echo "130 234" | gmt plot -Sr6p/6p -G#578fc1
echo "136 234" | gmt plot -Sr6p/6p -G#4b87bd
echo "142 234" | gmt plot -Sr6p/6p -G#407fb8
echo "148 234" | gmt plot -Sr6p/6p -G#3477b3
echo "154 234" | gmt plot -Sr6p/6p -G#296faf
echo "160 234" | gmt plot -Sr6p/6p -G#1d68aa
echo "166 234" | gmt plot -Sr6p/6p -G#1160a6
echo "172 234" | gmt plot -Sr6p/6p -G#0658a1
echo "55 229 0" | gmt text -F+f6p+jTC
echo "115 229 50" | gmt text -F+f6p+jTC
echo "175 229 100%" | gmt text -F+f6p+jTC
gmt basemap -R0/75/0.5/4.5 -JX345p/175p -X55p -Y45p -Bpxa15+l"Age band (years)" \
  -Bpycdata/classes-y.txt+l"Group" -BWS
# BEGIN EXAMPLE
gmt plot data/cells.tsv -i0,1 -Sr13.8p/43.75p -G#ececed
awk '$4>0 {print $1,$2,$3}' data/cells.tsv |
  gmt plot -Sr13.8p/43.75p -Cdata/survival.cpt
# END EXAMPLE
gmt end