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

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

# Mislocation vectors, with event numbers for the longer ones
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    13.668    18.138
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     8.435     0.590
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     6.660     3.541
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    11.834    -5.269
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     4.876    -7.509
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     5.908   -13.917
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.000   -11.400
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     8.422   -15.193
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.000   -11.195
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     6.568   -26.344
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    14.349   -11.620
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     6.315     0.000
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     9.002   -11.946
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    12.877    -4.687
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     8.922    -2.899
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.380   -21.755
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    19.128     1.002
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    10.832    -1.910
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    19.409    -4.481
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     5.464   -14.235
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.768   -13.140
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    10.296    -2.001
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    11.635    -5.180
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    26.300    -5.590
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     8.729    -8.429
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     5.655   -17.404
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.734   -24.794
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    14.762     2.075
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     6.010   -14.159
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.770   -13.033
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.066    -3.804
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.875    -6.541
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.745   -21.345
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    14.619   -13.163
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     5.648    -7.495
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    13.324    -9.680
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.514   -12.332
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.000   -18.868
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.065   -30.507
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     1.065   -30.507 39
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     6.418    -6.418
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.811    -3.028
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     5.814   -11.921
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     6.184   -13.890
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    18.036    -3.834
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.501   -16.469
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    20.494   -18.453
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    10.765   -24.178
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     4.458   -22.932
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    13.136    -8.861
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.903   -18.106
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     8.711    -8.123
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    12.322    -2.845
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     9.441   -15.109
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     6.367   -15.000
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.247   -11.869
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.450   -21.785
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    20.630    -8.335
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    23.151   -15.034
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    12.335    -5.492
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    11.233   -22.046
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.580   -12.485
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.147    -4.219
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.305     3.192
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -5.288     1.617
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     5.033    -7.751
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    13.562   -11.379
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     9.369   -21.044
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    13.923   -12.103
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     6.846   -14.036
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    12.457    -6.076
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     4.150    -6.153
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    15.047    -0.789
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    13.182    -2.802
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     7.572   -18.743
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.576   -16.499
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.769    -0.679
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     8.351    -9.952
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.701     6.349
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.999   -11.421
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.265    -6.507
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     4.896   -13.451
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.252   -15.301
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    14.127    -8.156
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     9.538   -16.520
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    10.717   -27.919
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    10.717   -27.919 85
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     6.500   -13.939
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    14.464   -16.639
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.906   -12.959
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    10.247   -14.635
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    12.527    -9.440
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.770   -14.700
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.822    -8.686
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    22.983   -35.391
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    22.983   -35.391 93
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     5.050    -8.747
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    12.691   -41.510
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    12.691   -41.510 95
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    19.710   -20.410
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     9.427   -13.975
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     8.154   -18.314
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    11.053   -15.214
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    27.578     6.367
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    20.740     3.285
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     8.674     7.540
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     9.066   -10.430
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    12.210   -11.791
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     9.870    -2.098
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     5.409    -0.189
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    11.705    -4.260
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     9.514    -4.848
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.650    -3.516
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     5.624     2.159
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    14.640   -11.032
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    17.333    -4.970
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    11.215     1.977
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    10.348    -5.502
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.779   -19.441
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     5.622    -1.195
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     6.342   -10.985
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    10.190    -1.797
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.683     0.207
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.802    -0.158
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.056    -3.561
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     6.106   -15.113
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     7.597   -13.159
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     8.925   -12.284
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     5.433    -0.860
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.160     7.533
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     5.905     4.781
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    10.646   -12.688
EOF
gmt plot $proj $reg -SE -Gred << EOF
     8.142    -9.638     0.000     2.000     2.000
EOF
gmt text $proj $reg -F+f10p,Helvetica,black+a0.+jBL << EOF
   -42.000    42.000 Author: all ( 128 samples)
EOF
gmt text $proj $reg -F+f10p,Helvetica,black+a0.+jBL << EOF
   -42.000   -42.000 Average =  12.6 km at 139.8 deg; median =   14.9 km
EOF
gmt end
