#!/bin/bash
gmt begin samar2.55_xsec_1 pdf,ps
proj=-JX18.0/-12.3
reg=-R-123.7/112.3/0/161.0
gmt basemap $proj $reg -Bxa10f10+l'Distance (km)' -Bya10f10+l'Depth (km)' -BWesN+t'Cross-section samar2.55'
# Event    1 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    53.170    46.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    53.170    43.000
    53.170    49.000
EOF
# Event    2 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -65.632    18.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -65.632    14.000
   -65.632    22.000
EOF
# Event    3 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -35.003    19.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -35.003    16.000
   -35.003    22.000
EOF
# Event    4 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    31.144    51.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    31.144    48.000
    31.144    54.000
EOF
# Event    5 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -31.523     9.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -31.523     5.000
   -31.523    13.000
EOF
# Event    6 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    -5.230   105.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    -5.230   102.000
    -5.230   108.000
EOF
# Event    7 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -66.567     6.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -66.567     2.000
   -66.567    10.000
EOF
# Event    8 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -57.745     3.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -57.745     0.000
   -57.745     7.000
EOF
# Event    9 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    24.405    60.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    24.405    57.000
    24.405    63.000
EOF
# Event   10 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -52.042    31.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -52.042    27.000
   -52.042    35.000
EOF
# Event   11 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -61.013     4.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -61.013     0.000
   -61.013     8.000
EOF
# Event   12 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    58.080    29.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    58.080    25.000
    58.080    33.000
EOF
# Event   13 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -77.263     5.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -77.263     1.000
   -77.263     9.000
EOF
# Event   14 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -52.172    10.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -52.172     6.000
   -52.172    14.000
EOF
# Event   15 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -54.950     6.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -54.950     2.000
   -54.950    10.000
EOF
# Event   16 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    90.620    25.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    90.620    22.000
    90.620    28.000
EOF
# Event   17 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    96.812     7.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    96.812     3.000
    96.812    11.000
EOF
# Event   18 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    -7.814   117.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    -7.814   114.000
    -7.814   120.000
EOF
# Event   19 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    51.704    44.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    51.704    41.000
    51.704    47.000
EOF
# Event   20 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    84.197    32.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    84.197    28.000
    84.197    36.000
EOF
# Event   21 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    84.387    31.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    84.387    27.000
    84.387    35.000
EOF
# Event   22 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -57.272     7.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -57.272     3.000
   -57.272    11.000
EOF
# Event   23 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    91.228    39.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    91.228    36.000
    91.228    42.000
EOF
# Event   24 fixed depth, contraint code w
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -37.059    12.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -37.059     7.000
   -37.059    17.000
EOF
# Event   25 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -40.891    16.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -40.891    13.000
   -40.891    19.000
EOF
# Event   26 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    30.443    61.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    30.443    58.000
    30.443    64.000
EOF
# Event   27 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
  -118.745    14.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
  -118.745    10.000
  -118.745    18.000
EOF
# Event   28 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -23.177   128.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -23.177   125.000
   -23.177   131.000
EOF
# Event   29 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    32.616    51.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    32.616    48.000
    32.616    54.000
EOF
# Event   30 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    50.241    57.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    50.241    54.000
    50.241    60.000
EOF
# Event   31 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -57.039    11.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -57.039     7.000
   -57.039    15.000
EOF
# Event   32 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    58.129    35.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    58.129    32.000
    58.129    38.000
EOF
# Event   33 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -24.673   156.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -24.673   153.000
   -24.673   159.000
EOF
# Event   34 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    60.613    34.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    60.613    31.000
    60.613    37.000
EOF
# Event   35 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    92.944    32.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    92.944    29.000
    92.944    35.000
EOF
# Event   36 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    93.750    29.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    93.750    26.000
    93.750    32.000
EOF
# Event   37 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    88.588    23.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    88.588    19.000
    88.588    27.000
EOF
# Event   38 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -81.875    18.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -81.875    15.000
   -81.875    21.000
EOF
# Event   39 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    38.012    43.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    38.012    40.000
    38.012    46.000
EOF
# Event   40 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    91.387    28.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    91.387    25.000
    91.387    31.000
EOF
# Event   41 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    85.930    32.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    85.930    28.000
    85.930    36.000
EOF
# Event   42 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -53.409    13.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -53.409     9.000
   -53.409    17.000
EOF
# Event   43 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    60.496    40.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    60.496    37.000
    60.496    43.000
EOF
# Event   44 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    40.904    44.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    40.904    41.000
    40.904    47.000
EOF
# Event   45 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    40.474    46.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    40.474    43.000
    40.474    49.000
EOF
# Event   46 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    73.481    29.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    73.481    26.000
    73.481    32.000
EOF
# Event   47 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -17.785    19.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -17.785    15.000
   -17.785    23.000
EOF
# Event   48 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    31.267    60.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    31.267    57.000
    31.267    63.000
EOF
# Event   49 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    93.485     9.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    93.485     5.000
    93.485    13.000
EOF
# Event   50 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    47.790    52.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    47.790    49.000
    47.790    55.000
EOF
# Event   51 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -46.379     3.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -46.379     0.000
   -46.379     7.000
EOF
# Event   52 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -58.669     6.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -58.669     3.000
   -58.669     9.000
EOF
# Event   53 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -69.508     5.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -69.508     1.000
   -69.508     9.000
EOF
# Event   54 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    27.011    59.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    27.011    56.000
    27.011    62.000
EOF
# Event   55 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    65.972    39.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    65.972    36.000
    65.972    42.000
EOF
# Event   56 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    64.148    37.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    64.148    34.000
    64.148    40.000
EOF
# Event   57 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   107.264    27.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   107.264    23.000
   107.264    31.000
EOF
# Event   58 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    45.585    41.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    45.585    38.000
    45.585    44.000
EOF
# Event   59 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   101.370    27.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   101.370    23.000
   101.370    31.000
EOF
# Event   60 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -67.005    22.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -67.005    18.000
   -67.005    26.000
EOF
# Event   61 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    40.510    42.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    40.510    39.000
    40.510    45.000
EOF
# Event   62 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    28.019    21.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    28.019    18.000
    28.019    24.000
EOF
# Event   63 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    96.459    23.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    96.459    19.000
    96.459    27.000
EOF
# Event   64 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -11.590    25.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -11.590    22.000
   -11.590    28.000
EOF
# Event   65 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    44.918    44.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    44.918    41.000
    44.918    47.000
EOF
# Event   66 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    98.096    16.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    98.096    12.000
    98.096    20.000
EOF
# Event   67 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -53.950     4.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -53.950     0.000
   -53.950     8.000
EOF
# Event   68 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    -3.939     9.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    -3.939     5.000
    -3.939    13.000
EOF
# Event   69 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -34.259    12.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -34.259     9.000
   -34.259    15.000
EOF
# Event   70 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -53.634    12.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -53.634     8.000
   -53.634    16.000
EOF
# Event   71 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -40.011     8.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -40.011     4.000
   -40.011    12.000
EOF
# Event   72 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
  -115.202    13.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
  -115.202     9.000
  -115.202    17.000
EOF
# Event   73 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
  -115.062    15.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
  -115.062    11.000
  -115.062    19.000
EOF
# Event   74 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    80.872    28.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    80.872    25.000
    80.872    31.000
EOF
# Event   75 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    87.138    27.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    87.138    23.000
    87.138    31.000
EOF
# Event   76 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -58.558    21.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -58.558    17.000
   -58.558    25.000
EOF
# Event   77 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -54.214    12.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -54.214     9.000
   -54.214    15.000
EOF
# Event   78 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    33.027    46.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    33.027    43.000
    33.027    49.000
EOF
# Event   79 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -44.221    13.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -44.221    10.000
   -44.221    16.000
EOF
# Event   80 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    14.426    90.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    14.426    87.000
    14.426    93.000
EOF
# Event   81 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -58.781     6.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -58.781     2.000
   -58.781    10.000
EOF
# Event   82 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -48.162    26.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -48.162    23.000
   -48.162    29.000
EOF
# Event   83 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    -6.989   106.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    -6.989   103.000
    -6.989   109.000
EOF
# Event   84 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    22.622    47.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    22.622    44.000
    22.622    50.000
EOF
# Event   85 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -56.344    10.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -56.344     6.000
   -56.344    14.000
EOF
# Event   86 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -51.359    17.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -51.359    13.000
   -51.359    21.000
EOF
# Event   87 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -56.328    32.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -56.328    28.000
   -56.328    36.000
EOF
# Event   88 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -51.282     7.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -51.282     4.000
   -51.282    10.000
EOF
# Event   89 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -62.335     8.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -62.335     4.000
   -62.335    12.000
EOF
# Event   90 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -55.890     8.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -55.890     4.000
   -55.890    12.000
EOF
# Event   91 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -69.244    18.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -69.244    14.000
   -69.244    22.000
EOF
# Event   92 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -71.291    13.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -71.291    10.000
   -71.291    16.000
EOF
# Event   93 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -52.320     4.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -52.320     1.000
   -52.320     7.000
EOF
# Event   94 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -58.505    15.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -58.505    12.000
   -58.505    18.000
EOF
# Event   95 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -55.817    17.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -55.817    14.000
   -55.817    20.000
EOF
# Event   96 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -79.721    13.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -79.721    10.000
   -79.721    16.000
EOF
# Event   97 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -54.266    10.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -54.266     6.000
   -54.266    14.000
EOF
# Event   98 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    30.211    57.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    30.211    54.000
    30.211    60.000
EOF
# Event   99 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    88.957    30.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    88.957    27.000
    88.957    33.000
EOF
# Event  100 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -53.690    13.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -53.690    10.000
   -53.690    16.000
EOF
# Event  101 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
     6.903    96.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
     6.903    93.000
     6.903    99.000
EOF
# Event  102 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    62.974    27.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    62.974    24.000
    62.974    30.000
EOF
# Event  103 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -70.051     0.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -70.051     0.000
   -70.051     4.000
EOF
# Event  104 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   104.723    31.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   104.723    27.000
   104.723    35.000
EOF
# Event  105 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    33.109    54.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    33.109    51.000
    33.109    57.000
EOF
# Event  106 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -69.903    20.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -69.903    16.000
   -69.903    24.000
EOF
# Event  107 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -68.791     3.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -68.791     0.000
   -68.791     7.000
EOF
# Event  108 fixed depth, contraint code l
gmt plot $proj $reg  -Sc -Wthick,black << EOF
    81.435    25.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
    81.435    21.000
    81.435    29.000
EOF
# Event  109 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
   -91.097    16.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
   -91.097    13.000
   -91.097    19.000
EOF
# Event  110 fixed depth, contraint code n
gmt plot $proj $reg  -Sc -Wthick,black << EOF
     5.261    90.000     0.200
EOF
gmt plot $proj $reg -Wthin,black << EOF
     5.261    87.000
     5.261    93.000
EOF

# Event numbers
gmt text $proj $reg -F+f8p,Helvetica,black+a0.+jBL << EOF
    53.370    46.000 1
   -65.432    18.000 2
   -34.803    19.000 3
    31.344    51.000 4
   -31.323     9.000 5
    -5.030   105.000 6
   -66.367     6.000 7
   -57.545     3.000 8
    24.605    60.000 9
   -51.842    31.000 10
   -60.813     4.000 11
    58.280    29.000 12
   -77.063     5.000 13
   -51.972    10.000 14
   -54.750     6.000 15
    90.820    25.000 16
    97.012     7.000 17
    -7.614   117.000 18
    51.904    44.000 19
    84.397    32.000 20
    84.587    31.000 21
   -57.072     7.000 22
    91.428    39.000 23
   -36.859    12.000 24
   -40.691    16.000 25
    30.643    61.000 26
  -118.545    14.000 27
   -22.977   128.000 28
    32.816    51.000 29
    50.441    57.000 30
   -56.839    11.000 31
    58.329    35.000 32
   -24.473   156.000 33
    60.813    34.000 34
    93.144    32.000 35
    93.950    29.000 36
    88.788    23.000 37
   -81.675    18.000 38
    38.212    43.000 39
    91.587    28.000 40
    86.130    32.000 41
   -53.209    13.000 42
    60.696    40.000 43
    41.104    44.000 44
    40.674    46.000 45
    73.681    29.000 46
   -17.585    19.000 47
    31.467    60.000 48
    93.685     9.000 49
    47.990    52.000 50
   -46.179     3.000 51
   -58.469     6.000 52
   -69.308     5.000 53
    27.211    59.000 54
    66.172    39.000 55
    64.348    37.000 56
   107.464    27.000 57
    45.785    41.000 58
   101.570    27.000 59
   -66.805    22.000 60
    40.710    42.000 61
    28.219    21.000 62
    96.659    23.000 63
   -11.390    25.000 64
    45.118    44.000 65
    98.296    16.000 66
   -53.750     4.000 67
    -3.739     9.000 68
   -34.059    12.000 69
   -53.434    12.000 70
   -39.811     8.000 71
  -115.002    13.000 72
  -114.862    15.000 73
    81.072    28.000 74
    87.338    27.000 75
   -58.358    21.000 76
   -54.014    12.000 77
    33.227    46.000 78
   -44.021    13.000 79
    14.626    90.000 80
   -58.581     6.000 81
   -47.962    26.000 82
    -6.789   106.000 83
    22.822    47.000 84
   -56.144    10.000 85
   -51.159    17.000 86
   -56.128    32.000 87
   -51.082     7.000 88
   -62.135     8.000 89
   -55.690     8.000 90
   -69.044    18.000 91
   -71.091    13.000 92
   -52.120     4.000 93
   -58.305    15.000 94
   -55.617    17.000 95
   -79.521    13.000 96
   -54.066    10.000 97
    30.411    57.000 98
    89.157    30.000 99
   -53.490    13.000 100
     7.103    96.000 101
    63.174    27.000 102
   -69.851     0.000 103
   104.923    31.000 104
    33.309    54.000 105
   -69.703    20.000 106
   -68.591     3.000 107
    81.635    25.000 108
   -90.897    16.000 109
     5.461    90.000 110
EOF

# User-defined stars
gmt plot $proj $reg -Sa -Wthicker,red -Gred << EOF
    53.170    46.000     0.200
   -57.745     3.000     0.200
   -51.359    17.000     0.200
EOF

# Azimuth of the projection plane
gmt text $proj $reg -F+f10p,Helvetica,black+a0.+jTR << EOF
   112.164     0.100 Azimuth =   60
EOF
gmt end
