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

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

# Mislocation vectors, with event numbers for the longer ones
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.916     5.197
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    18.475     5.298
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    18.475     5.298 2
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    19.719   -26.168
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    19.719   -26.168 3
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     8.502    -1.195
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     8.502    -1.195 4
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.480     0.854
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    22.364    -4.754
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    22.364    -4.754 6
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     7.115    -3.318
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     7.115    -3.318 7
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.286    -1.623
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.135    -0.612
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.906   -10.175
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     3.906   -10.175 10
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.253    -2.813
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.150    -1.002
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.144    -3.432
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.516     0.475
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.920    -5.811
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     4.339    -2.307
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.243    -1.989
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.414    -2.347
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     4.283    -5.684
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.977    -0.820
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.609    -4.151
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.581    -0.307
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.165     2.056
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.647     2.383
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     6.422    -4.666
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     6.422    -4.666 25
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     4.249    -1.981
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     6.045    -1.964
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.945    -4.366
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     5.608    -1.822
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     5.672    -3.972
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     7.249     1.409
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     7.249     1.409 31
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.605    -2.516
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.128    -1.460
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.479    -3.027
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.143    -6.168
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.784    -4.415
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     4.109    -2.004
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.882    -1.419
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.224    -2.224
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.983    -3.941
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.735    -1.731
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.884    -2.019
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.519    -2.082
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.085    -2.427
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.796     0.807
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.175    -2.509
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.813    -2.239
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.915     5.895
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.056     0.315
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.486     0.052
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     7.184    -2.903
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     7.184    -2.903 51
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.350     0.245
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.397    -2.251
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.796    -2.565
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.140    -1.335
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.267    -3.492
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     4.237    -2.648
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.783    -2.853
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.709    -4.171
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.233    -4.029
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.769    -3.300
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.534    -2.965
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.817    -1.471
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.125    -3.579
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.283     0.402
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.618    -2.908
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.614    -1.056
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.975    -2.774
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     5.660    -6.286
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     5.660    -6.286 69
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.622    -2.405
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.889    -2.315
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.201    -0.280
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     4.230    -2.542
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.937    -3.267
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.568    -2.716
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.883    -3.541
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.002    -0.861
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.905    -1.905
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.152     2.475
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.785    -3.920
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     4.415     5.079
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     4.979     1.911
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.499    -2.587
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.363    -0.265
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.126    -2.821
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.389    -1.047
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.726     1.309
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.209    -3.722
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.158    -1.654
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.525    -2.992
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.213    -1.364
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     4.858    -0.255
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.533    -1.827
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     4.089     0.071
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.207    -0.892
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.954    -1.420
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     4.800     0.420
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.185    -0.735
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.953    -1.220
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.242     0.000
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.142    -1.625
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.810    -0.839
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.222    -1.132
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     4.285    -0.989
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.947    -0.947
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.974    -5.891
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     5.962    -7.912
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     5.962    -7.912 107
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.995    -3.446
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.584    -0.139
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.709     0.760
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.871     1.957
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.733    -1.902
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     4.680    -3.039
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.533    -0.682
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.361    -3.432
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.145    -1.186
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.029     0.668
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     6.251    -1.675
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.435   -13.779
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -3.435   -13.779 119
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.218     3.167
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.030     1.586
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     6.234    -0.987
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.207    -1.049
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     5.425    -0.957
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.598    -5.989
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     9.047     3.473
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     9.047     3.473 126
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.070    -1.800
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.691    -1.800
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.370    -2.631
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.810    -0.546
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.100    -3.837
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     6.816    -0.837
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     4.537     2.115
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.467     3.122
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.899    -2.833
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     6.218    -0.109
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.940     1.809
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     7.982    -4.608
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     7.982    -4.608 138
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     4.085     1.650
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    10.869    -0.190
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    10.869    -0.190 140
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.890     0.989
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     4.334    -1.001
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.850     2.615
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.962     0.849
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.961     2.061
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     4.899    -0.688
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.257     1.018
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.877     0.339
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.482    -6.894
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     4.829     2.568
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     5.280     0.931
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.098    -1.861
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     4.782    -4.157
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     6.390     1.712
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.689     0.452
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     6.385    -2.580
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.428     2.701
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.495    -2.020
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     6.538     4.578
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     6.538     4.578 159
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.193     2.756
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.856     0.918
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.864    -1.275
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     8.964    -1.420
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     8.964    -1.420 163
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     4.555     0.399
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.252    -2.633
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.130     2.219
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.062     1.670
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.163     0.838
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.223    -3.193
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     5.831    -0.102
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.167     0.701
EOF
gmt plot $proj $reg -SE -Gred << EOF
     2.915    -1.473     0.000     2.000     2.000
EOF
gmt text $proj $reg -F+f10p,Helvetica,black+a0.+jBL << EOF
   -32.000    32.000 Author: all ( 171 samples)
EOF
gmt text $proj $reg -F+f10p,Helvetica,black+a0.+jBL << EOF
   -32.000   -32.000 Average =   3.3 km at 116.8 deg; median =    3.6 km
EOF
gmt end
