#!/bin/bash
gmt begin charyn3.19_epicenter_mislocation pdf,ps
proj=-JX12
reg=-R-30/30/-30/30
gmt basemap $proj $reg -Bxa25f5+l'km' -Bya25f5+l'km' -BWeSn+t'Epicentral Mislocation charyn3.19'

# Circle with radius of median vector length
gmt plot $proj $reg -SE -Wthick,red << EOF
     0.000     0.000     0.000    12.284    12.284
EOF

# Mislocation vectors, with event numbers for the longer ones
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.554     2.221
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.663    12.528
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -2.663    12.528 2
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -8.398     5.453
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -6.725     9.604
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.992    28.309
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -4.992    28.309 5
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.230     5.823
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -9.306    21.923
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -9.306    21.923 7
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.781     2.123
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.199     0.112
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -6.396     7.358
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
   -13.011    16.067
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
   -13.011    16.067 11
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.905     8.615
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.983    -3.024
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.421     3.668
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.784     5.022
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.915     7.594
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.955     6.394
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.828     9.466
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -8.209    -6.186
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -7.033    11.705
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -7.033    11.705 20
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
   -12.118    10.168
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
   -12.118    10.168 21
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.972     8.633
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.724     8.272
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.258     7.375
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
   -12.625     3.148
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
   -12.625     3.148 25
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.848     0.710
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.473    11.652
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -4.473    11.652 27
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -5.236    22.679
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -5.236    22.679 28
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -8.250    17.692
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -8.250    17.692 29
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
   -12.477    20.765
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
   -12.477    20.765 30
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
   -10.224    18.445
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
   -10.224    18.445 31
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -7.773    18.311
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -7.773    18.311 32
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
   -10.054    19.733
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
   -10.054    19.733 33
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -9.243     0.971
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -5.689    14.082
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -5.689    14.082 35
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.983     7.816
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.520    15.913
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -2.520    15.913 37
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.254     7.368
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.281     5.353
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.017    25.365
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -4.017    25.365 40
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.409     5.260
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -5.921     4.626
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -7.621    16.343
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -7.621    16.343 43
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.130    -2.472
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.724     2.921
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.467    -0.563
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.466     5.809
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.004     4.454
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.865     3.746
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.420    -0.059
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.579     4.851
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.067     8.289
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.037     5.308
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.171     1.395
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.714     4.896
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.905     0.294
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.216     0.298
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.825     5.544
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.964     3.599
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.658     4.679
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.054     3.045
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.424     1.353
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.328     8.238
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.334     4.362
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.128     4.176
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -6.026     5.425
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.455     2.078
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.280     4.685
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.431     6.158
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.477     5.925
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.250     4.889
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.358     4.181
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.000     8.089
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.168     4.806
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.008     7.002
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.584     3.003
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.121     7.078
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.329     4.711
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -5.058     3.285
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.341     6.309
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.266     4.418
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.241     4.605
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.592     5.630
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.208     4.736
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.179    16.862
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -1.179    16.862 85
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.984     4.910
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.964     1.977
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.073     4.154
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.582     5.540
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.557     5.302
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.155     4.447
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.451     7.463
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.049     3.659
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.913     4.664
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.342     1.781
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.321    10.917
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.717    10.253
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.929     0.676
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.786     4.355
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.275     5.112
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -7.034    11.706
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -7.034    11.706 101
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.035     3.617
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.712     2.276
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.179     0.972
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -6.906    12.989
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -6.906    12.989 105
EOF
gmt plot $proj $reg -SE -Gred << EOF
    -3.118     6.860     0.000     2.000     2.000
EOF
gmt text $proj $reg -F+f10p,Helvetica,black+a0.+jBL << EOF
   -28.000    28.000 Author: all ( 105 samples)
EOF
gmt text $proj $reg -F+f10p,Helvetica,black+a0.+jBL << EOF
   -28.000   -28.000 Average =   7.5 km at 335.6 deg; median =    6.1 km
EOF
gmt end
