#!/usr/bin/env bash
set -euo pipefail
export GMT_USERDIR="$PWD/.gmt"
mkdir -p "$GMT_USERDIR"
gmt begin VIZ-11-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-11-005
ID
gmt text -F+f7p+jML <<'LEGEND'
55 246 Colour: survival 0-100%; circle area: n
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
echo "280 235" | gmt plot -Sc3.975732839729454p -G#00549f
echo "280 227 n=10" | gmt text -F+f6p+jTC
echo "325 235" | gmt plot -Sc8.890008890013334p -G#00549f
echo "325 227 n=50" | gmt text -F+f6p+jTC
echo "370 235" | gmt plot -Sc12.572371141874243p -G#00549f
echo "370 227 n=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 variable circle diameters are centimetres; 72 bp = 1 inch.
awk '$4>0 {print $1,$2,$3,$6/72*2.54}' data/cells.tsv | gmt plot -Sc -Cdata/survival.cpt
# END EXAMPLE
gmt end