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

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

# Mislocation vectors, with event numbers for the longer ones
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.787     0.156
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -9.807     3.570
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -9.807     3.570 2
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.923    -0.599
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.415     0.249
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.751     0.864
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -5.518    -2.934
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -5.518    -2.934 6
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.449    -0.239
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -7.656    -1.076
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -7.656    -1.076 8
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.725    -0.060
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.072     1.531
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
   -13.088    -4.764
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
   -13.088    -4.764 11
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.140     3.367
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -3.140     3.367 12
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.688    -0.516
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.084    -0.270
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.812    -0.071
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.721     0.262
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.906    -0.160
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.775     0.014
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     4.708     5.611
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     4.708     5.611 19
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.708    -0.112
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.554     0.235
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.107     0.874
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.060    -0.056
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.387     2.309
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -1.387     2.309 24
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.341     0.525
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.954     0.151
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.759     2.251
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     1.759     2.251 27
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.098     0.933
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.248     1.767
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.403     0.326
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.616     0.370
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.622     0.359
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.318     0.874
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.403     0.123
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.653    -0.602
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.216    -0.000
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.118     0.020
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.023     0.036
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.413     1.144
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.590     2.364
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     0.590     2.364 40
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.884     3.083
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     0.884     3.083 41
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.787     1.212
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.064     1.392
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -2.064     1.392 43
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.624     1.625
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.777     0.544
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.024     1.375
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.048     1.377
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.120     1.940
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.918     4.724
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     0.918     4.724 49
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.476     1.558
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.042    -0.236
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.152     0.868
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.328     0.616
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.224     0.527
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.383    -0.527
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.367     1.886
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.014    -0.389
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.505    -0.840
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.646     0.198
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.394    -1.374
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.512    -0.239
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.583     0.695
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.535     0.764
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.560    -0.204
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.068     0.521
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.938     0.066
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.458    -0.427
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.712    -0.164
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.532    -2.451
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -1.532    -2.451 69
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.829     1.863
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.924    -1.503
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -1.924    -1.503 71
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.386    -0.291
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.284    -1.335
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.352     3.523
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     1.352     3.523 74
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.071     0.806
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.295     2.101
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.211    -1.501
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.379     0.193
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.006     0.317
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.156     0.334
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.501     3.220
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     1.501     3.220 81
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.303    -0.374
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.603    -0.817
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.567     0.253
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.483     1.938
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.593     0.659
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.797     0.112
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.761     3.580
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     0.761     3.580 88
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.024    -1.355
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.259     1.334
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.438     0.523
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.495    -1.603
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.508     0.108
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.724     0.560
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.951     0.572
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.761     0.907
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.707    -0.012
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.781     0.547
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.426     0.382
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.643     0.990
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.201     1.914
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.138    -0.348
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.678     0.407
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.416    -0.426
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -2.416    -0.426 104
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.416    -1.029
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.716     0.483
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.913    -1.088
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.092    -0.011
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.018     0.273
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.311    -0.115
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.695     1.756
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -1.695     1.756 111
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.271    -1.275
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.189    -0.115
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.238     1.503
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.659    -0.321
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.449    -0.104
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.890     0.955
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.353    -0.318
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.968    -1.075
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.491    -0.398
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.489     0.104
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.081     2.323
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.017    -0.712
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.309     0.023
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.149     0.141
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.679     0.450
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.830     0.207
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.745    -1.344
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.764    -0.879
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.470     0.126
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.093    -0.154
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.907     0.111
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.184     0.029
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.396    -0.099
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.422    -0.015
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.545     0.278
EOF
gmt plot $proj $reg -SE -Gred << EOF
    -0.762     0.399     0.000     1.500     1.500
EOF
gmt text $proj $reg -F+f10p,Helvetica,black+a0.+jBL << EOF
   -13.500    13.500 Author: all ( 136 samples)
EOF
gmt text $proj $reg -F+f10p,Helvetica,black+a0.+jBL << EOF
   -13.500   -13.500 Average =   0.9 km at 297.6 deg; median =    1.2 km
EOF
gmt end
