#!/usr/bin/env bash

# =========================
# GMT settings
# =========================
gmt set FONT_ANNOT_PRIMARY   12p,Helvetica
gmt set FONT_ANNOT_SECONDARY 12p,Helvetica
gmt set FONT_LABEL           12p,Helvetica
gmt set FONT_TITLE           12p,Helvetica

# =========================
# Plot
# =========================
gmt begin gmt_map pdf

    # CPT
    gmt makecpt -Cjet -T0/100/1 -Z -D

    # Map
    gmt basemap -R120/150/20/50 -JM12c -Bafg -BWSNE
    gmt coast   -Df -W0.25p

    # Main focal mechanisms
    gmt meca _cmt_2011_gmt_5.0_sort \
        -Sa0.5c \
        -C \
        -Fa0.05c

    # Colorbar
    gmt colorbar \
        -DJBR+jBL+o2c/0c+w8c/0.4c \
        -Baf+l"Depth [km]" \
        --FONT_ANNOT=16p \
        --FONT_LABEL=16p

    # Legend focal mechanisms
    gmt meca -Sa0.5c -C -Fa0.05c -N << EOF
155.5 48 80 226 47 136 9.0
155.5 46 80 226 47 136 8.0
155.5 44 80 226 47 136 7.0
155.5 42 80 226 47 136 6.0
155.5 40 80 226 47 136 5.0
EOF

    # Legend labels
    gmt text -N << EOF
157.5 48 M9
157.5 46 M8
157.5 44 M7
157.5 42 M6
157.5 40 M5
EOF

gmt end

xdg-open gmt_map.pdf