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

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

# Mislocation vectors, with event numbers for the longer ones
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.785    -3.147
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     8.284     4.405
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     8.284     4.405 2
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.084     4.788
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.590    -0.377
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     4.184     3.269
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.015    -5.952
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -4.015    -5.952 6
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.597     5.332
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -3.597     5.332 7
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.942     1.071
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.106     5.603
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     3.106     5.603 9
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.489     3.052
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.836     7.950
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     0.836     7.950 11
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.330    -0.122
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.155     5.941
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -1.155     5.941 13
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.193     1.097
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -8.464     0.741
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -8.464     0.741 15
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.342    -6.527
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     0.342    -6.527 16
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.986    -2.945
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.076    -3.455
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.744    -3.032
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.381    -0.298
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.459     1.459
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.129    -2.537
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.956    -2.016
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.682    -3.929
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -4.682    -3.929 24
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.283    -0.626
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -5.343    -2.268
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.255    -0.078
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.421    -7.030
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -2.421    -7.030 28
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.865    -2.253
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -9.378    -5.860
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -9.378    -5.860 30
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.421    -2.164
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.226     2.152
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.624    -1.339
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.925     0.556
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.005    -2.920
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.915    -5.187
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.301    -1.257
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.585     1.313
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.479     0.463
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.682     3.865
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.079     3.600
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.661    10.634
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -3.661    10.634 42
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -7.453     4.657
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -7.453     4.657 43
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.873    -1.641
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.691     1.751
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.933     2.198
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.127    -2.480
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.002    -3.083
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.366    -0.118
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.420    -0.187
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.654    -2.307
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.218    -1.504
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.972    -1.686
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.903    -3.369
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.374     0.291
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.047     0.280
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.304    -2.353
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.057    -2.904
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.399    -3.843
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.721     0.101
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.082    -2.156
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.199    -2.844
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.572    -2.693
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.862    -2.030
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.467     0.748
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.123     0.648
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.294    -1.199
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.957    -0.101
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.286    -2.758
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.000    -4.830
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.230    -1.574
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.785    -2.205
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.384    -2.044
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.552    -0.318
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.381     0.897
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.177    -0.283
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
   -11.845    -1.454
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
   -11.845    -1.454 77
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.544    -2.558
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.334    -3.177
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.033    -1.127
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.177    -1.507
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.638     0.878
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -5.954    -1.157
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -5.954    -1.157 83
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.032    -2.928
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.654     1.283
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.000    -1.165
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.116     0.261
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.529    -1.246
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.452    -1.218
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.355     0.424
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.582     0.969
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.862    -1.555
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.942    -2.486
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.266    -1.065
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.617    -0.789
EOF
gmt plot $proj $reg -SE -Gred << EOF
    -0.850    -0.596     0.000     1.500     1.500
EOF
gmt text $proj $reg -F+f10p,Helvetica,black+a0.+jBL << EOF
   -11.500    11.500 Author: all (  95 samples)
EOF
gmt text $proj $reg -F+f10p,Helvetica,black+a0.+jBL << EOF
   -11.500   -11.500 Average =   1.0 km at 235.0 deg; median =    3.0 km
EOF
gmt end
