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

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

# Mislocation vectors, with event numbers for the longer ones
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
   -16.848     2.069
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    11.422   -39.833
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    11.422   -39.833 2
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
   -14.244     9.250
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -7.720    -7.720
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
   -16.013     5.828
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.043     4.869
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
   -15.681    -3.333
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -6.450    -9.933
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
   -14.048     3.243
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.887    11.917
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -6.687    12.577
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    27.866     3.422
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    27.866     3.422 12
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
   -12.693    28.508
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
   -12.693    28.508 13
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
   -33.984     4.776
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
   -33.984     4.776 14
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.554    -3.145
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.678    -5.744
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     4.111    -3.574
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
   -17.185     2.415
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -8.214    10.144
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.962    -4.558
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.275     0.352
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.410    21.530
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     9.292     0.650
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -7.137    16.814
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
   -11.164     1.173
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     8.315    -7.228
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    15.751     6.046
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    21.189     0.740
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -7.998     2.911
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.071    21.849
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -6.091    21.243
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     8.909    -8.308
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     8.241    -7.958
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     5.740   -18.774
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
   -10.873    14.429
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     6.480     9.978
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    84.935   -11.937
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    84.935   -11.937 37
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
   -15.871    -4.852
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
   -10.844    15.487
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
   -15.634    11.359
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     9.950     7.498
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
   -20.405    21.881
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
   -20.405    21.881 42
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    31.947   -11.628
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    31.947   -11.628 43
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    16.661   -11.238
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
   -15.368    16.480
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.448     2.054
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.000     4.811
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -5.915     6.343
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.201    20.946
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.983    10.204
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     8.950    11.456
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -9.869    22.166
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -9.869    22.166 52
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.968    10.351
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    18.543     9.860
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.803    -0.857
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
   -13.182     3.043
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.334    19.135
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     5.956    -3.167
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.987    -2.251
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.160    -2.484
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.020    -3.799
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -6.512     4.731
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    46.903   -28.182
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    46.903   -28.182 63
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.205     5.876
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.750    -1.876
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.031     0.437
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.905    -2.352
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.698    10.721
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.307    -4.877
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
   -23.179     6.646
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.842    -5.916
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -9.042    -4.026
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     8.623     9.247
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     5.336     0.750
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.810    -1.088
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -6.668    11.097
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     4.527   -23.287
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.000     9.370
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    12.236     9.560
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     5.897    -4.129
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.877    -0.789
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.633    -2.539
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.842     3.340
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.031    -1.229
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -6.976    -5.854
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -7.241    17.058
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.188     1.532
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -8.803    -6.633
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
   -16.663    18.506
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
   -16.663    18.506 89
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    14.335   -20.473
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    14.335   -20.473 90
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    16.834    16.834
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
   -14.891    11.221
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
   -35.488    33.093
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
   -35.488    33.093 93
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     8.458     3.417
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     5.166    12.787
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.243     1.977
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.210    12.048
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    25.891     3.639
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    25.891     3.639 98
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -6.005    14.863
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
   -19.979     6.879
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -7.593     4.037
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
   -13.246     1.862
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -6.429     2.340
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.032   -12.161
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
   -13.521    11.753
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -7.563     0.529
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
   -11.689    13.930
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     7.306     5.506
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.864     5.121
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     4.854     1.767
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.649     3.683
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.985    10.211
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     4.764     7.336
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     4.153    12.782
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     7.764     8.931
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.295     8.430
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.089    10.774
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.341    13.359
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.078     2.667
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.444     5.303
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.801     3.000
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     4.996    11.770
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.552     5.254
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     5.120    10.048
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     7.517    11.145
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.850    27.392
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     3.850    27.392 126
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -5.950    10.305
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
   -13.141    11.026
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.037    16.591
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.309    17.727
EOF
gmt plot $proj $reg -SE -Gred << EOF
    -0.523     4.556     0.000     2.000     2.000
EOF
gmt text $proj $reg -F+f10p,Helvetica,black+a0.+jBL << EOF
   -85.000    85.000 Author: all ( 130 samples)
EOF
gmt text $proj $reg -F+f10p,Helvetica,black+a0.+jBL << EOF
   -85.000   -85.000 Average =   4.6 km at 353.4 deg; median =   12.1 km
EOF
gmt end
