#!/bin/bash
gmt begin samar2.55_epicenter_mislocation pdf,ps
proj=-JX12
reg=-R-23/23/-23/23
gmt basemap $proj $reg -Bxa10f5+l'km' -Bya10f5+l'km' -BWeSn+t'Epicentral Mislocation samar2.55'

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

# Mislocation vectors, with event numbers for the longer ones
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.525     2.699
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.414    -2.552
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.852     3.349
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.945    -4.872
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.850     5.364
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     6.777    -0.474
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.349    -0.689
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.525     3.972
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.124    -3.565
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
   -13.385    -7.728
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
   -13.385    -7.728 10
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -5.136    -1.669
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.435    -6.223
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.553     0.970
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.665    -0.664
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.506     5.138
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     5.792    -0.814
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.546    -2.190
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.734     9.831
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    10.206     8.264
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    10.206     8.264 19
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     6.126     1.991
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.694     4.562
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.939    -5.769
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.320     0.754
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -5.063    -3.164
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.697     5.530
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.358     0.071
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.186     1.531
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.245     4.672
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.170     4.855
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     7.984     1.122
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -7.988    -0.981
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.208     3.440
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.031     1.961
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     5.691     7.028
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    15.657     9.040
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    15.657     9.040 35
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.006    -1.079
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.867     0.506
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.249     1.379
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.050     7.474
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    16.339    14.203
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    16.339    14.203 40
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     9.279    -0.324
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.044    -2.116
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.802     0.656
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.722     5.139
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.408     4.867
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    12.955    11.665
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    12.955    11.665 46
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     5.724    -1.113
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.945     0.670
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     5.087     4.120
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.474     4.320
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.898    -2.937
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
   -10.731    -5.004
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
   -10.731    -5.004 52
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.799     1.977
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.018    -0.177
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    13.461     3.108
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    13.461     3.108 55
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     8.083    12.936
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     8.083    12.936 56
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.904     0.617
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -8.164     4.160
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     4.943     1.997
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.953     1.845
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.124     3.919
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
   -16.218    -3.153
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
   -16.218    -3.153 62
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     7.919     6.884
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.924     3.197
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.229    -6.323
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     6.599     9.783
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     6.599     9.783 66
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.525    -3.525
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.942   -10.259
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.256    -5.256
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.913     1.060
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.723    -6.877
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.604    -0.897
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -5.025    -0.528
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     4.851     2.915
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     7.760     3.294
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -8.517     1.810
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -6.412    -3.410
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.550     7.614
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -7.148    -3.333
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     6.412     4.164
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.776    -2.329
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     6.283     0.660
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     7.227     2.488
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.320     2.239
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.519     1.229
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -7.818    -1.099
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.382     4.177
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.179     3.072
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.509     1.065
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.608     1.096
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     6.973     1.355
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     6.199     9.546
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.173     3.173
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    10.653     0.745
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.748     2.303
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -7.530     4.004
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.381    -5.155
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.738     5.981
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    14.606    11.007
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    14.606    11.007 99
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.088     3.198
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     5.929     0.728
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.428     4.979
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.177     5.071
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    13.995     8.745
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    13.995     8.745 104
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.645     0.358
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.824     1.230
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     7.954     9.150
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     7.954     9.150 107
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     4.302     5.313
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.656     1.930
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.194    -1.216
EOF
gmt plot $proj $reg -SE -Gred << EOF
     1.009     1.816     0.000     1.500     1.500
EOF
gmt text $proj $reg -F+f10p,Helvetica,black+a0.+jBL << EOF
   -21.500    21.500 Author: all ( 110 samples)
EOF
gmt text $proj $reg -F+f10p,Helvetica,black+a0.+jBL << EOF
   -21.500   -21.500 Average =   2.1 km at  29.1 deg; median =    5.8 km
EOF
gmt end
