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

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

# Mislocation vectors, with event numbers for the longer ones
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    10.567     8.557
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    10.567     8.557 1
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     9.828     2.633
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     9.828     2.633 2
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    26.884   -13.112
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    26.884   -13.112 3
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    13.444   -16.602
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    13.444   -16.602 4
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.751     7.167
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     2.751     7.167 5
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     5.130     6.114
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     5.130     6.114 6
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.204    -7.686
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -2.204    -7.686 7
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.933    -0.148
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.682     3.086
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.757     4.926
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -4.757     4.926 10
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     4.302     5.709
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     4.302     5.709 11
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.819     9.367
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -0.819     9.367 12
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.093     3.576
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.135    -2.456
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.049     3.430
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.771     0.062
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.744    -6.996
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -1.744    -6.996 17
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.284    -0.558
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.767     0.121
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.820     3.882
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.038    -1.482
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.257    -1.729
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.305    -2.560
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.895    -3.588
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.190     0.465
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.972    -0.759
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.922    -1.556
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.327    -0.093
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.902    -0.841
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.029    -0.721
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.052    -2.966
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -6.608    -0.115
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -6.608    -0.115 32
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.715    -4.427
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.172     2.403
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.183    -2.611
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.622    -0.898
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.266    -0.694
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.130    -3.711
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.011    -1.342
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.321     0.458
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.046     2.608
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.670    -1.043
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.975     4.887
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.605    -3.430
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.641     4.558
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.846    -0.275
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.319    -2.271
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.366     2.206
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.000    -2.215
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.046    -2.623
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.476     0.752
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.678     0.678
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.340    -2.417
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.448    -6.377
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -2.448    -6.377 54
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.947    -1.816
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.301    -0.448
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.221    -2.504
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.137    -2.615
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.449     1.381
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.222    -0.863
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.762    -2.932
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.422    -0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.532    -1.643
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.069    -2.575
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.108    -2.058
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.073     4.171
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.611    -0.656
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.250    -1.544
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.135    -1.551
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.895     2.180
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.481    -2.727
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -5.608    -2.266
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -5.608    -2.266 72
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.842     2.194
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.893    -3.572
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.533    -2.742
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -5.496    -0.969
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.057    -5.439
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.387    -0.251
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.643    -2.401
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.734     2.640
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.341    -0.248
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.306    -3.838
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.663     1.299
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.835    -0.492
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.171     2.440
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.831    -2.789
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.395    -3.219
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.463     2.624
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.360    -0.416
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.650    -2.023
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.075    -6.035
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -3.075    -6.035 91
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.360    -1.224
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.958     0.833
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.139    -1.926
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.171     0.761
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.038    -2.199
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.748     1.244
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.932     2.307
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.252    -1.139
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     5.443    -7.773
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     5.443    -7.773 100
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.922    -0.625
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.612    -4.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.087    -4.095
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.547    -1.253
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.161     1.535
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.555     0.269
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.035     0.239
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     6.927    -0.728
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     6.927    -0.728 108
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.301    -0.334
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     5.172     0.819
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.392    -0.738
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.783     3.693
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.667     3.391
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.367    -0.721
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.445    -4.237
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.545    -2.564
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.845    -3.389
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.957     0.069
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.295    -2.811
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.366     0.112
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.918     1.882
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.869    -2.670
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     4.708    -5.610
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     4.708    -5.610 123
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -6.590    -3.653
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -6.590    -3.653 124
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.148    -2.372
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.549     5.068
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.886    -3.095
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.064     0.876
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.633    -5.025
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -6.694    -4.863
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -6.694    -4.863 130
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     8.896     4.730
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     8.896     4.730 131
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.393    -0.883
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.702    -2.592
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.197     3.683
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.526    -1.835
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.123    -2.355
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -5.498    11.791
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -5.498    11.791 137
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.549    -1.394
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.109    -1.248
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.235    -1.324
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.815     7.758
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     0.815     7.758 141
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     4.438     3.594
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.424    -1.581
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.599    -0.489
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.078    -4.452
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.500     9.331
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     2.500     9.331 146
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.895     3.885
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.524     2.271
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.155    -1.650
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.941    -1.162
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.250     3.250
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.894    -0.924
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.830    -2.733
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.920    -2.830
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.149     1.641
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     4.303    -2.192
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -5.001    -2.772
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.712     1.763
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.306     3.690
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.302    -2.168
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.825    -0.674
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.335    -0.076
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     4.672     1.165
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.858    -0.419
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.315    -2.473
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.248     0.429
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.453     2.861
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.300    -0.703
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     6.643     0.232
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     6.643     0.232 169
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.905     5.701
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     2.905     5.701 170
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.292     9.193
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     2.292     9.193 171
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.009    -1.880
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.864     1.076
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.378     0.474
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -7.300    -2.949
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -7.300    -2.949 175
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.039    -3.624
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.858    -2.760
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     9.584     7.488
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     9.584     7.488 178
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     4.202     5.378
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     4.202     5.378 179
EOF
gmt plot $proj $reg -SE -Gred << EOF
     0.090    -0.366     0.000     2.000     2.000
EOF
gmt text $proj $reg -F+f10p,Helvetica,black+a0.+jBL << EOF
   -29.000    29.000 Author: all ( 179 samples)
EOF
gmt text $proj $reg -F+f10p,Helvetica,black+a0.+jBL << EOF
   -29.000   -29.000 Average =   0.4 km at 166.2 deg; median =    3.0 km
EOF
gmt end
