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

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

# Mislocation vectors, with event numbers for the longer ones
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.525     4.993
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     4.375    -1.091
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -5.062   -11.924
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -5.062   -11.924 3
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     7.516    -1.874
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     7.516    -1.874 4
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.882    -3.540
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.754    -4.132
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.856     0.381
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.632    -0.511
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.623    -3.535
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -5.356     4.184
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -5.356     4.184 10
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.431    -2.443
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.902     4.571
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -4.902     4.571 12
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.903    -5.119
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.282    -4.678
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.751    -3.591
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.225    -1.600
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.193    -3.672
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.202     3.132
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.496    -2.813
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.648     0.029
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.310    -0.635
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.134     3.831
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.376     4.321
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.679    -0.234
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.072    -0.168
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.386     0.121
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.000     2.128
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.265     2.040
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.548     1.134
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.639    -1.022
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.373    -0.998
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.219    -2.219
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.085     0.250
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.150    -0.496
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.910    -0.476
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.118     1.123
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.965    12.404
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     1.965    12.404 38
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.331     4.643
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.561     0.801
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.527    -0.292
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.675    -0.643
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.153    -0.420
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.823    -1.374
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.329    -1.399
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.391    -2.466
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.763    -1.322
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.215     0.730
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.769    -1.727
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.509     0.132
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.156     1.749
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.179     2.256
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.459     0.399
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.304     0.389
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.068     3.910
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.749    -0.535
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.148     3.622
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.481    -1.375
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.649    -3.052
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.040    -0.227
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.857    -1.252
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.780     0.397
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.576    -2.239
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.082    -0.799
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.545    -2.362
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.375     0.263
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.044    -1.336
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.834    -0.088
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.575    -0.984
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.549    -2.549
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.265    -3.488
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.194    -0.434
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.651    -0.749
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.871     0.398
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.975    -1.489
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.818    -2.818
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.352     0.172
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.165    -4.072
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.384    -0.507
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.199    -1.995
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.847    -1.467
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.366    -1.657
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.002    -2.147
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.548    -0.000
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.128    -4.362
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.963    -0.963
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.634    -5.345
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.034     0.389
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.403     0.108
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.175    -3.334
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.821    -1.622
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.153    -0.537
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.402    -3.307
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.015    -1.758
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.041     1.169
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.321    -1.510
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.585    -3.431
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.343    -6.031
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     3.343    -6.031 98
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.117    -0.025
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.774    -1.239
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.267    -2.348
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.457    -1.835
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.439    -4.246
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.043    -2.645
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.599    -2.181
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.378    -4.323
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.398     1.018
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.970    -1.524
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.874    -2.874
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.056     1.614
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.759     1.149
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.079    -2.759
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.736    -1.218
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.256    -0.516
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.482    -3.247
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.726    -1.348
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.111    -2.801
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.567     0.567
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.985     0.532
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.352    -1.974
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.023    -0.425
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.203    -3.283
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.442    -4.209
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.841    -5.058
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.338    -2.507
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.613    -1.504
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     6.330    -4.770
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     6.330    -4.770 127
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.128    -1.462
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.000    -2.115
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.168    -6.008
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.324    -1.523
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.894    -1.138
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.809    -3.880
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.266     5.077
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     4.569    -1.310
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -3.375     1.948
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.185    -2.138
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.208    -2.971
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.599    -4.003
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.240    -2.744
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.817    -5.278
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.236    -1.493
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.248    -2.829
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.471    -4.457
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     6.378    -0.558
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     6.378    -0.558 145
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.144     2.058
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -5.799    -5.799
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -5.799    -5.799 147
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.328    -1.858
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.807    -1.893
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     4.416    -3.839
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.262    -3.665
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -4.202    -1.285
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.062    -3.539
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     4.910    -2.950
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.445    -6.719
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -2.445    -6.719 155
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.923    -2.357
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.390    -3.274
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.584    -9.012
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -2.584    -9.012 158
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.270    -0.483
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.710   -13.555
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     0.710   -13.555 160
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.108    -0.747
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -1.461    -2.749
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -2.119    -6.932
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    -2.119    -6.932 163
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.111    -6.132
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     2.111    -6.132 164
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     3.438    -4.097
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     0.673    -7.689
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
     0.673    -7.689 166
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    -0.333    -4.768
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     1.202    -3.698
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
     2.603    -2.603
EOF
gmt plot $proj $reg -Wthin,black << EOF
     0.000     0.000
    17.053     5.541
EOF
gmt text $proj $reg -F+f8p+jCM << EOF
    17.053     5.541 170
EOF
gmt plot $proj $reg -SE -Gred << EOF
     0.931    -1.490     0.000     1.500     1.500
EOF
gmt text $proj $reg -F+f10p,Helvetica,black+a0.+jBL << EOF
   -17.500    17.500 Author: all ( 170 samples)
EOF
gmt text $proj $reg -F+f10p,Helvetica,black+a0.+jBL << EOF
   -17.500   -17.500 Average =   1.8 km at 148.0 deg; median =    3.1 km
EOF
gmt end
