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

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

# Mislocation vectors, with event numbers for the longer ones
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.791    -7.757
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
   -16.974    -3.608
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
   -16.974    -3.608 2
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -9.813    -3.572
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.116     7.342
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -5.271    -8.116
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    18.571     9.462
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    18.571     9.462 6
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.922    -8.776
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.444    11.763
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -6.712     9.586
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    16.720    13.063
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    16.720    13.063 10
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.617    10.495
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.623    -5.891
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.357    -6.883
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -7.020    -0.491
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -7.266    -2.789
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.093    -7.699
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.973    -0.632
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.208    -2.598
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    13.580    -6.046
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    13.580    -6.046 19
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     9.801     7.937
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     9.801     7.937 20
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    32.438     1.133
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    32.438     1.133 21
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.483    -0.482
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.574    -2.894
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    10.412     0.728
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -5.346    -2.160
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    17.768    10.258
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    17.768    10.258 26
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    12.513     0.875
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    12.513     0.875 27
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.798    -1.440
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     6.071    -3.648
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.211    -0.895
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -7.229     2.921
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     8.829    -2.869
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
   -21.262   -11.786
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
   -21.262   -11.786 33
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.980     1.349
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.944    -4.534
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    11.856    -2.090
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -6.222    -4.041
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.665    12.692
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -0.665    12.692 38
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     9.219     2.299
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.961     0.842
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     5.430     0.285
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.273     2.709
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.934    -3.219
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.546     1.424
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.833     9.524
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     6.439     1.969
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.435     1.746
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    11.335    -1.392
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     6.825    -6.591
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.433    -0.232
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     9.353    -0.983
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.549    -6.310
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.621     2.714
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.167     2.381
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.523     0.636
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -5.533     2.349
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.066     7.584
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.427     1.396
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.789     7.509
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.862    -0.821
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     8.697    -3.692
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     6.251     0.657
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.003     5.893
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -6.253     2.400
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.935     0.908
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.249    -7.844
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.096    -2.462
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     4.493     3.770
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.370    -3.010
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.056    -1.112
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.286    -0.120
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.151    12.775
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -4.151    12.775 72
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.117     1.430
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.315     1.809
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     5.335     2.837
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     7.777     5.650
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.078     1.156
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.941    -1.565
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    22.656     6.071
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    22.656     6.071 79
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.364    -0.167
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.314     0.133
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     5.077    -0.089
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.860    -0.790
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.405     4.323
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.831     5.248
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.428    -6.185
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.500     0.158
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -6.806    -7.299
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     5.111     0.993
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.854    -2.968
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.931    -1.827
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.317    -5.281
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -7.296    12.637
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -7.296    12.637 93
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.057   -12.084
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -5.775     2.942
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     7.664    -3.738
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.946     1.939
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    12.940    -6.311
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    12.940    -6.311 98
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.501    -9.334
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     6.356     4.287
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.041    -1.182
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.859    -1.254
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.579     8.436
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.687    -0.297
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    21.706     5.816
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    21.706     5.816 105
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.000    -2.098
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     6.154     3.698
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     4.214    -2.633
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.277     0.569
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.626     0.057
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.125    -2.208
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.329     0.894
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     6.490    -3.165
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -5.774     1.988
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.743     4.067
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.456     6.523
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -6.484     1.497
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -8.340    -1.621
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.629     3.384
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.532    -1.140
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.993    -2.955
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.100     2.252
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.345     5.366
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -8.316    -1.617
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -5.215     7.448
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.520     2.653
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.120     3.336
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.740     0.499
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.521     0.834
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.306     0.295
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -5.252     6.259
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.481     5.500
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.570    -3.219
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.462    20.053
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -2.462    20.053 134
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.286     0.573
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.917     2.817
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -6.230     4.868
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.174     7.857
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.604    10.468
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.883    -1.261
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.193     1.833
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.631    -1.982
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.448     5.301
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     5.091     7.270
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.293    -0.913
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     6.349    -2.961
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.210    -0.742
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.788     1.360
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.321     6.213
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.801     2.612
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    12.049    -5.365
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    12.049    -5.365 151
EOF
gmt plot $proj $reg -SE -Gred << EOF
     1.007     0.873     0.000     2.000     2.000
EOF
gmt text $proj $reg -F+f10p,Helvetica,black+a0.+jBL << EOF
   -31.000    31.000 Author: all ( 151 samples)
EOF
gmt text $proj $reg -F+f10p,Helvetica,black+a0.+jBL << EOF
   -31.000   -31.000 Average =   1.3 km at  49.1 deg; median =    6.1 km
EOF
gmt end
