#!/usr/bin/env bash
#
# Generated by:
# taup curve --mod prem -h 500 -p s,S,ScS,Sdiff --gmt -o stdout 
# The TauP Toolkit
# version 3.1.0
# https://github.com/crotwell/TauP
#
# for use with GMT6, modern syntax
# https://docs.generic-mapping-tools.org/latest/reference/introduction.html#modern-and-classic-mode
#
#
gmt begin taup_curve

gmt basemap -JX6.0i/6.0i -R0.0/160.87334025390402/108.95337178722247/1929.835529979536 -Bxa20.0+l'Degrees' -Bya200.0+l'Time (s)' -BWSne+t'prem (h=500.0 km)' 
# s
gmt plot  -l's' -W,skyblue <<END
> 1/1 s degree180/time s for a source depth of 500.0 kilometers in the prem model
10.6403  264.191
8.86877  229.559
8.31289  218.811
8.26321  217.855
8.10656  214.847
8.09936  214.710
8.09220  214.572
8.07795  214.299
8.04976  213.759
7.99463  212.704
7.98169  212.457
7.96883  212.211
7.94337  211.724
7.89343  210.771
7.54849  204.216
7.33417  200.173
7.02813  194.447
6.91584  192.361
6.75804  189.445
6.68233  188.052
6.66925  187.812
6.54199  185.481
6.51542  184.997
6.46946  184.159
6.41987  183.257
6.30567  181.189
6.29465  180.990
6.28406  180.799
6.27351  180.609
6.18100  178.945
6.09174  177.347
6.07808  177.103
6.06448  176.860
6.03751  176.379
5.98439  175.435
5.93235  174.512
5.88133  173.610
5.78221  171.865
5.68672  170.195
5.59457  168.592
5.54967  167.815
5.50553  167.054
5.41937  165.574
5.33590  164.149
5.22828  162.325
5.20127  161.869
5.12479  160.585
5.08734  159.959
5.05041  159.344
4.97799  158.144
4.86238  156.243
4.83598  155.811
4.83409  155.781
4.83219  155.750
4.82841  155.688
4.82086  155.565
4.82039  155.557
4.81603  155.486
4.81168  155.415
4.80299  155.274
4.78569  154.993
4.75142  154.438
4.61755  152.286
4.52884  150.876
4.44284  149.523
4.43146  149.345
4.42013  149.168
4.39759  148.816
4.35304  148.124
4.31387  147.518
4.29449  147.219
4.27524  146.923
4.23711  146.339
4.21824  146.051
4.19949  145.766
4.15968  145.162
4.13997  144.864
4.12040  144.569
4.08557  144.046
4.06832  143.787
4.05116  143.531
4.01716  143.025
4.00030  142.775
3.98354  142.527
3.94805  142.004
3.91299  141.490
3.85064  140.582
3.84742  140.535
3.78958  139.701
3.73285  138.890
3.67718  138.102
3.61358  137.209
3.58228  136.774
3.55130  136.345
3.49027  135.507
3.48937  135.495
3.48847  135.482
3.48667  135.458
3.48487  135.433
3.48307  135.409
3.48077  135.377
3.47847  135.346
3.47617  135.315
3.47388  135.283
3.46930  135.221
3.46473  135.159
3.45859  135.075
3.45247  134.992
3.45212  134.988
3.44882  134.943
3.44487  134.889
3.44093  134.836
3.43305  134.730
3.42523  134.624
3.41742  134.518
3.40966  134.414
3.40191  134.310
3.40107  134.299
3.40024  134.287
3.39856  134.265
3.39689  134.242
3.39522  134.220
3.39187  134.175
3.38853  134.130
3.38187  134.041
3.37522  133.952
3.36190  133.774
3.34865  133.598
3.32618  133.300
3.31501  133.152
3.30388  133.005
3.28173  132.714
3.25974  132.425
3.21621  131.859
3.19096  131.533
3.16590  131.211
3.11638  130.580
3.06761  129.965
3.06377  129.917
3.05994  129.868
3.05228  129.773
3.03703  129.582
3.03284  129.530
3.02867  129.478
3.01011  129.248
3.00088  129.134
2.99167  129.020
2.97076  128.763
2.94999  128.509
2.91890  128.131
2.91746  128.114
2.91601  128.096
2.88238  127.691
2.84908  127.292
2.82710  127.031
2.80183  126.733
2.78203  126.501
2.77923  126.468
2.77644  126.435
2.77086  126.370
2.75972  126.241
2.73757  125.984
2.73573  125.962
2.73390  125.941
2.71370  125.708
2.69362  125.478
2.69076  125.446
2.68791  125.413
2.66514  125.154
2.64252  124.899
2.59772  124.397
2.59068  124.319
2.57644  124.161
2.56226  124.005
2.52715  123.620
2.49238  123.244
2.44724  122.761
2.42694  122.546
2.41255  122.394
2.39821  122.244
2.36971  121.947
2.34142  121.656
2.34019  121.643
2.32393  121.477
2.30773  121.312
2.27555  120.988
2.24364  120.669
2.23785  120.612
2.23208  120.555
2.22056  120.441
2.21640  120.401
2.21225  120.360
2.19681  120.209
2.18144  120.060
2.17951  120.041
2.17758  120.022
2.17373  119.985
2.16435  119.894
2.15499  119.804
2.14566  119.715
2.13635  119.626
2.11779  119.450
2.09932  119.276
2.06263  118.934
2.05873  118.898
2.05483  118.862
2.04235  118.747
2.02991  118.633
2.00513  118.408
1.97575  118.144
1.95113  117.925
1.93640  117.796
1.92173  117.667
1.91244  117.587
1.90781  117.546
1.90318  117.506
1.88889  117.383
1.88176  117.322
1.87464  117.261
1.85828  117.122
1.84197  116.984
1.83712  116.944
1.82983  116.882
1.82900  116.876
1.82817  116.869
1.82652  116.855
1.82321  116.827
1.81990  116.800
1.81659  116.772
1.81329  116.745
1.80999  116.717
1.80339  116.663
1.79461  116.590
1.78585  116.519
1.77711  116.447
1.76838  116.376
1.76043  116.312
1.75250  116.248
1.73668  116.121
1.71777  115.970
1.69893  115.822
1.68640  115.724
1.67391  115.627
1.65842  115.508
1.64299  115.390
1.62427  115.248
1.62306  115.239
1.62186  115.230
1.61944  115.212
1.61461  115.176
1.60545  115.107
1.59631  115.039
1.57808  114.905
1.55834  114.761
1.55352  114.726
1.52907  114.550
1.50474  114.378
1.48053  114.209
1.47825  114.193
1.47100  114.143
1.46225  114.083
1.45351  114.023
1.42888  113.857
1.40437  113.694
1.39423  113.627
1.37996  113.534
1.36851  113.459
1.35708  113.386
1.33973  113.276
1.32244  113.167
1.31144  113.099
1.28467  112.934
1.26728  112.829
1.25816  112.774
1.25498  112.755
1.23968  112.665
1.22442  112.575
1.20785  112.480
1.19132  112.385
1.18702  112.361
1.18273  112.337
1.16238  112.223
1.14209  112.112
1.13785  112.089
1.13442  112.070
1.13099  112.051
1.11668  111.975
1.10954  111.937
1.10240  111.899
1.08800  111.823
1.07362  111.749
1.06859  111.723
1.06356  111.698
1.04451  111.601
1.02550  111.507
1.02201  111.490
1.02043  111.482
1.01886  111.474
1.00348  111.400
0.988139  111.326
0.966397  111.224
0.962373  111.205
0.957966  111.185
0.953561  111.165
0.942465  111.114
0.931385  111.064
0.924295  111.032
0.918461  111.007
0.917848  111.004
0.917234  111.001
0.916008  110.996
0.914511  110.989
0.913014  110.983
0.910021  110.969
0.907642  110.959
0.905263  110.949
0.900509  110.928
0.895757  110.907
0.891008  110.887
0.885115  110.862
0.879226  110.836
0.867460  110.787
0.856854  110.743
0.851556  110.721
0.846261  110.699
0.833400  110.647
0.820559  110.595
0.801334  110.520
0.797261  110.504
0.797201  110.504
0.775132  110.420
0.753055  110.338
0.731029  110.259
0.709098  110.182
0.709051  110.182
0.695455  110.135
0.688115  110.111
0.666766  110.040
0.646272  109.975
0.625410  109.910
0.605395  109.850
0.591949  109.811
0.585009  109.791
0.565390  109.736
0.545769  109.683
0.526176  109.632
0.508281  109.586
0.490409  109.543
0.487678  109.536
0.468744  109.492
0.449833  109.450
0.420545  109.387
0.420141  109.386
0.419737  109.385
0.419332  109.385
0.418928  109.384
0.418120  109.382
0.417312  109.381
0.416504  109.379
0.415696  109.377
0.414081  109.374
0.412465  109.371
0.410850  109.367
0.409234  109.364
0.407619  109.361
0.406004  109.358
0.404389  109.355
0.402775  109.351
0.401160  109.348
0.399546  109.345
0.397931  109.342
0.394703  109.336
0.394174  109.335
0.393644  109.334
0.393115  109.333
0.392586  109.331
0.391527  109.329
0.390468  109.327
0.388943  109.325
0.387418  109.322
0.385313  109.318
0.383209  109.314
0.379001  109.306
0.374793  109.298
0.370587  109.290
0.366382  109.283
0.362177  109.275
0.357974  109.268
0.353771  109.261
0.345446  109.246
0.337123  109.232
0.328804  109.219
0.320489  109.206
0.313308  109.194
0.306130  109.183
0.298953  109.173
0.291779  109.162
0.287534  109.156
0.279364  109.145
0.271197  109.134
0.263032  109.123
0.254870  109.113
0.238662  109.093
0.222462  109.075
0.208235  109.060
0.194013  109.046
0.190276  109.042
0.174274  109.028
0.158278  109.015
0.142355  109.003
0.126437  108.993
0.111644  108.984
0.0968558  108.976
0.0947195  108.975
0.0789058  108.969
0.0630949  108.963
0.0473127  108.959
0.0315321  108.956
0.0157658  108.954
0.00000  108.953
END
# end of s
# S
gmt plot  -l'S' -W,olivedrab <<END
> 1/1 S degree180/time S for a source depth of 500.0 kilometers in the prem model
10.6403  264.191
11.6826  284.616
12.7250  304.936
13.5409  320.711
13.6207  322.246
13.8707  327.046
13.8828  327.278
13.8945  327.503
13.9179  327.950
13.9641  328.835
14.0551  330.577
14.0768  330.992
14.0983  331.403
14.1409  332.217
14.2252  333.826
14.8302  345.324
15.2341  352.943
15.8439  364.352
16.0833  368.800
16.4266  375.145
16.5979  378.297
16.6282  378.854
16.9214  384.221
18.6854  416.420
19.9185  438.889
20.8992  456.718
22.6867  489.095
22.8423  491.903
22.9884  494.541
23.1319  497.129
24.3216  518.532
25.3907  537.675
24.2417  517.152
23.0928  496.637
22.2306  481.254
21.0958  461.026
19.6465  435.248
18.6446  417.483
17.8657  403.711
16.6763  382.776
15.7718  366.951
15.0383  354.196
14.7171  348.636
14.4198  343.505
13.8841  334.302
13.4110  326.227
12.8540  316.786
12.7221  314.562
12.3635  308.540
13.4726  327.098
14.5817  345.628
15.4149  359.509
16.5257  377.930
17.7874  398.682
18.0369  402.757
18.0546  403.045
18.0719  403.327
18.1062  403.886
18.1739  404.990
18.1781  405.059
18.2168  405.689
18.2553  406.315
18.3311  407.548
18.4792  409.955
18.7631  414.557
19.7796  430.898
20.4140  440.980
21.0061  450.299
22.0282  466.305
23.0503  482.301
23.9417  496.236
25.2462  516.589
27.1767  546.594
28.5706  568.150
29.2014  577.871
29.7984  587.048
30.9128  604.122
31.4402  612.172
31.9491  619.921
33.1166  637.637
33.6453  645.628
34.1558  653.323
35.0368  666.557
35.4646  672.960
35.8831  679.210
36.6963  691.313
37.0945  697.217
37.4856  703.004
38.4227  716.818
39.2667  729.196
40.7141  750.272
40.7904  751.378
42.0957  770.204
43.5089  790.409
44.7923  808.593
46.2264  828.714
47.0400  840.040
47.7744  850.207
49.1832  869.553
49.2058  869.863
49.2271  870.153
49.2691  870.726
49.3108  871.295
49.3523  871.861
49.4053  872.584
49.4581  873.304
49.5109  874.023
49.5636  874.740
49.6683  876.166
49.7728  877.587
49.9129  879.490
50.0525  881.386
50.0719  881.649
50.1729  883.018
50.2793  884.460
50.3815  885.844
50.5799  888.527
50.7736  891.142
50.9644  893.715
51.1525  896.249
51.3390  898.756
51.3593  899.029
51.3794  899.300
51.4196  899.839
51.4597  900.378
51.4997  900.915
51.5796  901.986
51.6592  903.054
51.8178  905.177
51.9757  907.290
52.2901  911.487
52.6023  915.646
53.1287  922.632
53.3909  926.101
53.6510  929.534
54.2532  937.457
54.8006  944.626
55.8587  958.394
56.4701  966.291
57.0722  974.029
58.3777  990.674
59.5935  1006.01
59.6920  1007.24
59.7877  1008.44
59.9775  1010.82
60.3543  1015.52
60.4915  1017.23
60.6085  1018.68
61.1022  1024.81
61.3440  1027.80
61.5827  1030.74
62.1191  1037.34
62.6499  1043.83
63.4395  1053.44
63.4777  1053.90
63.5148  1054.35
64.4510  1065.64
65.3300  1076.15
65.9074  1083.00
66.5670  1090.79
67.1416  1097.53
67.2195  1098.44
67.2962  1099.34
67.4485  1101.12
67.7502  1104.63
68.3460  1111.54
68.3965  1112.13
68.4461  1112.70
68.9863  1118.92
69.5235  1125.08
69.6201  1126.18
69.7046  1127.15
70.3442  1134.42
70.9690  1141.48
72.1932  1155.19
72.3905  1157.38
72.8150  1162.08
73.2175  1166.52
74.1997  1177.27
75.1712  1187.80
76.4780  1201.78
77.0603  1207.95
77.4710  1212.27
77.8792  1216.55
78.7136  1225.24
79.5309  1233.66
79.5684  1234.05
80.0382  1238.85
80.5065  1243.62
81.4496  1253.13
82.3823  1262.43
82.5550  1264.14
82.7247  1265.82
83.0623  1269.15
83.1850  1270.35
83.3071  1271.55
83.7593  1275.97
84.2119  1280.37
84.2700  1280.93
84.3271  1281.49
84.4409  1282.59
84.7174  1285.25
84.9941  1287.91
85.2705  1290.56
85.5467  1293.20
86.0870  1298.33
86.6337  1303.48
87.7229  1313.64
87.8441  1314.76
87.9617  1315.84
88.3195  1319.13
88.6875  1322.50
89.4242  1329.20
90.3076  1337.13
91.0267  1343.52
91.4703  1347.42
91.9135  1351.30
93.0378  1361.10
94.1621  1370.88
94.7424  1375.91
95.2527  1380.33
96.6088  1392.02
97.2052  1397.14
97.7675  1401.95
98.9707  1412.18
100.081  1421.55
100.403  1424.25
100.873  1428.20
END
# end of S
# ScS
gmt plot  -l'ScS' -W,goldenrod <<END
> 1/1 ScS degree180/time ScS for a source depth of 500.0 kilometers in the prem model
100.873  1428.20
99.2896  1414.95
97.7059  1401.72
96.4245  1391.01
94.6426  1376.13
93.4121  1365.86
92.1816  1355.60
90.3389  1340.25
88.8151  1327.58
87.4945  1316.62
86.3178  1306.87
85.2924  1298.38
84.2671  1289.91
83.1125  1280.39
81.9579  1270.88
79.9729  1254.60
78.2160  1240.24
76.6302  1227.35
75.3018  1216.58
74.0654  1206.61
72.9399  1197.56
71.8144  1188.55
70.6149  1178.99
69.4155  1169.47
68.3405  1160.99
67.2655  1152.54
65.9414  1142.20
64.6910  1132.50
63.2251  1121.20
61.8445  1110.66
60.2646  1098.69
60.1656  1097.95
60.0671  1097.21
59.8711  1095.74
59.4837  1092.83
58.7642  1087.45
58.0642  1082.25
56.7181  1072.32
55.3276  1062.17
54.9979  1059.78
53.3792  1048.15
51.8483  1037.30
50.3943  1027.15
50.2608  1026.23
49.8396  1023.32
49.3379  1019.87
48.8445  1016.50
47.4911  1007.35
46.1942  998.712
45.6712  995.270
44.9480  990.547
44.3775  986.854
43.8168  983.251
42.9811  977.935
42.1658  972.808
41.6563  969.637
40.4421  962.178
39.6722  957.524
39.2742  955.141
39.1365  954.320
38.4793  950.430
37.8337  946.653
37.1435  942.663
36.4657  938.794
36.2913  937.807
36.1176  936.827
35.3031  932.277
34.5051  927.892
34.3398  926.992
34.2068  926.271
34.0742  925.554
33.5243  922.602
33.2522  921.155
32.9818  919.725
32.4404  916.889
31.9059  914.124
31.7203  913.172
31.5355  912.229
30.8406  908.719
30.1566  905.323
30.0320  904.711
29.9758  904.435
29.9197  904.161
29.3752  901.517
28.8371  898.942
28.0834  895.400
27.9450  894.758
27.7938  894.059
27.6431  893.366
27.2651  891.640
26.8902  889.948
26.6515  888.880
26.4558  888.011
26.4352  887.920
26.4147  887.830
26.3737  887.648
26.3236  887.427
26.2736  887.207
26.1737  886.768
26.0945  886.421
26.0154  886.075
25.8575  885.387
25.7001  884.705
25.5432  884.028
25.3490  883.196
25.1556  882.372
24.7710  880.750
24.4262  879.313
24.2546  878.605
24.0836  877.903
23.6702  876.224
23.2598  874.582
22.6501  872.188
22.5216  871.690
22.5197  871.683
21.8275  869.046
21.1417  866.505
20.4636  864.063
19.7945  861.722
19.7930  861.717
19.3810  860.311
19.1594  859.566
18.5184  857.454
17.9075  855.502
17.2900  853.590
16.7016  851.826
16.3083  850.678
16.1060  850.097
15.5363  848.499
14.9697  846.963
14.4071  845.490
13.8957  844.198
13.3874  842.958
13.3099  842.772
12.7741  841.520
12.2413  840.323
11.4205  838.574
11.4092  838.551
11.3979  838.528
11.3867  838.505
11.3754  838.482
11.3528  838.435
11.3302  838.389
11.3077  838.343
11.2851  838.297
11.2401  838.206
11.1950  838.115
11.1499  838.024
11.1049  837.933
11.0599  837.843
11.0149  837.754
10.9699  837.665
10.9250  837.576
10.8800  837.487
10.8351  837.399
10.7902  837.311
10.7005  837.137
10.6857  837.109
10.6710  837.080
10.6563  837.052
10.6416  837.024
10.6122  836.967
10.5828  836.911
10.5405  836.830
10.4981  836.750
10.4397  836.639
10.3813  836.529
10.2647  836.311
10.1481  836.096
10.0317  835.884
9.91541  835.674
9.79920  835.466
9.68311  835.261
9.56713  835.059
9.33762  834.666
9.10852  834.283
8.87982  833.909
8.65150  833.546
8.45459  833.240
8.25795  832.941
8.06158  832.650
7.86547  832.365
7.74951  832.201
7.52655  831.891
7.30390  831.590
7.08156  831.298
6.85951  831.016
6.41921  830.483
5.97996  829.986
5.59482  829.578
5.21037  829.199
5.10942  829.103
4.67758  828.717
4.24647  828.365
3.81783  828.048
3.38978  827.765
2.99237  827.533
2.59537  827.330
2.53804  827.303
2.11386  827.122
1.68999  826.974
1.26709  826.859
0.844387  826.777
0.422162  826.728
0.00000  826.712
END
# end of ScS
# Sdiff
gmt plot  -l'Sdiff' -W,firebrick <<END
> 1/1 Sdiff degree180/time Sdiff for a source depth of 500.0 kilometers in the prem model
100.873  1428.20
102.873  1444.92
104.873  1461.64
106.873  1478.36
108.873  1495.08
110.873  1511.80
112.873  1528.52
114.873  1545.24
116.873  1561.97
118.873  1578.69
120.873  1595.41
122.873  1612.13
124.873  1628.85
126.873  1645.57
128.873  1662.29
130.873  1679.02
132.873  1695.74
134.873  1712.46
136.873  1729.18
138.873  1745.90
140.873  1762.62
142.873  1779.34
144.873  1796.06
146.873  1812.79
148.873  1829.51
150.873  1846.23
152.873  1862.95
154.873  1879.67
156.873  1896.39
158.873  1913.11
160.873  1929.84
END
# end of Sdiff
gmt end 
