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

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

# Mislocation vectors, with event numbers for the longer ones
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -6.633     1.777
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -6.930     1.473
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.100     1.164
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -7.900     3.853
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -7.900     3.853 4
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.697    -4.201
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     8.334    -4.247
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     8.334    -4.247 6
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.551     3.921
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -5.151    -1.977
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.664     1.446
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.623    -4.436
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
   -10.922    -4.636
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
   -10.922    -4.636 11
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -7.162    -2.190
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.889    -5.161
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.604    -0.091
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.211     3.693
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.131    -3.753
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
   -13.214     8.582
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
   -13.214     8.582 17
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -6.825    -3.039
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -5.286     2.693
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -5.318     0.279
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -8.368    -0.732
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -8.368    -0.732 21
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -6.605     1.046
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.538    -2.628
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.700    -1.420
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.658    -6.258
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -7.691    -9.844
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -7.691    -9.844 26
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -9.007    -1.915
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -9.007    -1.915 27
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.084    -0.286
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -5.570    -2.250
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -5.165    -0.634
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.308    -5.506
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -6.041    -3.212
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -5.389     6.199
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -5.389     6.199 33
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -7.850     0.687
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -7.142     2.320
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -6.129    -3.259
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -6.751     2.866
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -6.250     2.914
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -5.144     2.399
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.857    -0.150
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.141     2.689
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.042     5.173
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.757     2.003
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.163    -2.308
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.796    -0.983
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.589     4.132
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.883     3.542
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -7.673     2.346
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -7.673     2.346 48
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -6.781     1.817
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.301     4.544
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -7.108    -1.511
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.818     1.473
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -7.691     2.351
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -7.691     2.351 53
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -6.350     2.311
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -7.315     2.808
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -6.214     7.954
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -6.214     7.954 56
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.890    -0.547
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.901     0.066
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.403    -3.308
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.491     2.267
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.121    -6.925
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.071     1.245
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.947    -4.104
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -6.975     0.244
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -5.963    -0.945
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -6.363    -2.191
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -9.884    -1.921
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -9.884    -1.921 67
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.864    -7.034
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -7.540    -1.603
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.664     2.876
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.691    -0.585
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.017    -2.163
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -6.456     2.350
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.553    -1.889
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.553     2.644
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.098    -5.642
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.982    -1.429
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.307    -2.403
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.395    -1.178
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.934    -5.833
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.036    -1.021
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -5.258     2.915
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.861    -2.079
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.228    -2.096
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.895     0.472
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.392    -1.014
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.097     1.211
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -8.987    -2.408
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -8.987    -2.408 88
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.263     0.536
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.657    -2.551
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.553    -2.302
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -7.528     2.159
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.262     1.736
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.488    -1.596
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -7.273     5.683
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -7.273     5.683 95
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.557     1.173
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.550     0.534
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.236     2.664
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.920     0.624
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.293     1.606
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.567     0.060
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.302     2.063
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.593    -0.896
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -7.696     9.172
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -7.696     9.172 104
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.440     2.121
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.147     2.286
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -5.522     0.974
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.308     0.763
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.122     0.383
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -5.390     4.365
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.660    -0.909
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.272    -1.145
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.444     2.811
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.792     1.131
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.914     1.584
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.286     1.532
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.015     0.319
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.066     2.093
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -5.441     6.719
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -5.441     6.719 119
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.284     1.896
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.370     4.230
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -5.796     0.815
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -6.074     2.211
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.183     1.484
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.437     0.727
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.481     0.903
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.019     0.109
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.382    -1.050
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.084    -0.211
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.117    -1.974
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -5.002     0.526
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.449     1.802
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.966     2.578
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.129    -1.985
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.208    -0.698
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.482    -0.569
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.835    -0.607
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.896    -0.518
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.815    -1.533
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.363    -0.000
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.495     1.728
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.681     0.894
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.314    -0.504
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.381    -2.739
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.414    -0.850
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.076     0.038
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.789    -1.348
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.276    -2.294
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.053     3.065
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.207     0.975
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.680    -1.974
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.794    -0.337
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.823    -0.811
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.240    -0.558
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.248     0.647
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.451    -0.389
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.107    -2.511
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.404    -0.049
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.607     3.460
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.344    -1.992
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.245    -0.325
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.529    -1.308
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.282    -4.027
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.518     1.513
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.263     0.286
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.214     0.295
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.669     2.866
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.796     2.604
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.772    -0.665
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.478    -3.410
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.970    -1.646
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -6.221     0.874
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.940     2.007
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.140    -0.130
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.454     1.124
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.242     0.437
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.446     2.839
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.110     0.999
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.590     5.746
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     5.953     4.995
EOF
gmt plot $proj $reg -SE -Gred << EOF
    -2.743     0.195     0.000     1.500     1.500
EOF
gmt text $proj $reg -F+f10p,Helvetica,black+a0.+jBL << EOF
   -15.500    15.500 Author: all ( 180 samples)
EOF
gmt text $proj $reg -F+f10p,Helvetica,black+a0.+jBL << EOF
   -15.500   -15.500 Average =   2.7 km at 274.1 deg; median =    3.9 km
EOF
gmt end
