#ifndef ONE_D_MODELS_H #define ONE_D_MODELS_H #include #include "config.h" // ak135 http://rses.anu.edu.au/seismology/ak135/ak135t.html // depth(km), P velocity(km/s) inline const std::vector> model_1d_ak135 \ {{ 0.000, 5.8000}, { 20.000, 5.8000}, { 20.000, 6.5000}, { 35.000, 6.5000}, { 35.000, 8.0400}, { 77.500, 8.0450}, { 120.000, 8.0500}, { 165.000, 8.1750}, { 210.000, 8.3000}, { 210.000, 8.3000}, { 260.000, 8.4825}, { 310.000, 8.6650}, { 360.000, 8.8475}, { 410.000, 9.0300}, { 410.000, 9.3600}, { 460.000, 9.5280}, { 510.000, 9.6960}, { 560.000, 9.8640}, { 610.000, 10.0320}, { 660.000, 10.2000}, { 660.000, 10.7900}, { 710.000, 10.9229}, { 760.000, 11.0558}, { 809.500, 11.1353}, { 859.000, 11.2221}, { 908.500, 11.3068}, { 958.000, 11.3896}, { 1007.500, 11.4705}, { 1057.000, 11.5495}, { 1106.500, 11.6269}, { 1156.000, 11.7026}, { 1205.500, 11.7766}, { 1255.000, 11.8491}, { 1304.500, 11.9200}, { 1354.000, 11.9895}, { 1403.500, 12.0577}, { 1453.000, 12.1245}, { 1502.500, 12.1912}, { 1552.000, 12.2550}, { 1601.500, 12.3185}, { 1651.000, 12.3819}, { 1700.500, 12.4426}, { 1750.000, 12.5031}, { 1799.500, 12.5631}, { 1849.000, 12.6221}, { 1898.500, 12.6804}, { 1948.000, 12.7382}, { 1997.500, 12.7956}, { 2047.000, 12.8526}, { 2096.500, 12.9096}, { 2146.000, 12.9668}, { 2195.500, 13.0222}, { 2245.000, 13.0783}, { 2294.500, 13.1336}, { 2344.000, 13.1894}, { 2393.500, 13.2465}, { 2443.000, 13.3018}, { 2492.500, 13.3585}, { 2542.000, 13.4156}, { 2591.500, 13.4741}, { 2640.000, 13.5312}, { 2690.000, 13.5900}, { 2740.000, 13.6494}, { 2740.000, 13.6494}, { 2789.670, 13.6530}, { 2839.330, 13.6566}, { 2891.500, 13.6602}, { 2891.500, 8.0000}, { 2939.330, 8.0382}, { 2989.660, 8.1283}, { 3039.990, 8.2213}, { 3090.320, 8.3122}, { 3140.660, 8.4001}, { 3190.990, 8.4861}, { 3241.320, 8.5692}, { 3291.650, 8.6496}, { 3341.980, 8.7283}, { 3392.310, 8.8036}, { 3442.640, 8.8761}, { 3492.970, 8.9461}, { 3543.300, 9.0138}, { 3593.640, 9.0792}, { 3643.970, 9.1426}, { 3694.300, 9.2042}, { 3744.630, 9.2634}, { 3794.960, 9.3205}, { 3845.290, 9.3760}, { 3895.620, 9.4297}, { 3945.950, 9.4814}, { 3996.280, 9.5306}, { 4046.620, 9.5777}, { 4096.950, 9.6232}, { 4147.280, 9.6673}, { 4197.610, 9.7100}, { 4247.940, 9.7513}, { 4298.270, 9.7914}, { 4348.600, 9.8304}, { 4398.930, 9.8682}, { 4449.260, 9.9051}, { 4499.600, 9.9410}, { 4549.930, 9.9761}, { 4600.260, 10.0103}, { 4650.590, 10.0439}, { 4700.920, 10.0768}, { 4801.580, 10.1415}, { 4851.910, 10.1739}, { 4902.240, 10.2049}, { 4952.580, 10.2329}, { 5002.910, 10.2565}, { 5053.240, 10.2745}, { 5103.570, 10.2854}, { 5153.500, 10.2890}, { 5153.500, 11.0427}, { 5204.610, 11.0585}, { 5255.320, 11.0718}, { 5306.040, 11.0850}, { 5356.750, 11.0983}, { 5407.460, 11.1166}, { 5458.170, 11.1316}, { 5508.890, 11.1457}, { 5559.600, 11.1590}, { 5610.310, 11.1715}, { 5661.020, 11.1832}, { 5711.740, 11.1941}, { 5813.160, 11.2134}, { 5863.870, 11.2219}, { 5914.590, 11.2295}, { 5965.300, 11.2364}, { 6016.010, 11.2424}, { 6066.720, 11.2477}, { 6117.440, 11.2521}, { 6168.150, 11.2557}, { 6218.860, 11.2586}, { 6269.570, 11.2606}, { 6320.290, 11.2618}, { 6371.000, 11.2622}}; // iasp91 model from http://ds.iris.edu/spud/earthmodel/9991809 inline const std::vector> model_1d_iasp91 \ {{0.00, 5.8000}, {1.00, 5.8000}, {2.00, 5.8000}, {3.00, 5.8000}, {4.00, 5.8000}, {5.00, 5.8000}, {6.00, 5.8000}, {7.00, 5.8000}, {8.00, 5.8000}, {9.00, 5.8000}, {10.00, 5.8000}, {11.00, 5.8000}, {12.00, 5.8000}, {13.00, 5.8000}, {14.00, 5.8000}, {15.00, 5.8000}, {16.00, 5.8000}, {17.00, 5.8000}, {18.00, 5.8000}, {19.00, 5.8000}, {20.00, 5.8000}, {20.00, 6.5000}, {21.00, 6.5000}, {22.00, 6.5000}, {23.00, 6.5000}, {24.00, 6.5000}, {25.00, 6.5000}, {26.00, 6.5000}, {27.00, 6.5000}, {28.00, 6.5000}, {29.00, 6.5000}, {30.00, 6.5000}, {31.00, 6.5000}, {32.00, 6.5000}, {33.00, 6.5000}, {34.00, 6.5000}, {35.00, 6.5000}, {35.00, 8.0400}, {40.00, 8.0406}, {45.00, 8.0412}, {50.00, 8.0418}, {60.00, 8.0429}, {70.00, 8.0441}, {80.00, 8.0453}, {90.00, 8.0465}, {100.00, 8.0476}, {110.00, 8.0488}, {120.00, 8.0500}, {120.00, 8.0500}, {130.00, 8.0778}, {140.00, 8.1056}, {150.00, 8.1333}, {160.00, 8.1611}, {170.00, 8.1889}, {180.00, 8.2167}, {190.00, 8.2444}, {200.00, 8.2722}, {210.00, 8.3000}, {210.00, 8.3000}, {220.00, 8.3365}, {230.00, 8.3730}, {240.00, 8.4095}, {250.00, 8.4460}, {260.00, 8.4825}, {270.00, 8.5190}, {280.00, 8.5555}, {290.00, 8.5920}, {300.00, 8.6285}, {310.00, 8.6650}, {320.00, 8.7015}, {330.00, 8.7380}, {340.00, 8.7745}, {350.00, 8.8110}, {360.00, 8.8475}, {370.00, 8.8840}, {380.00, 8.9205}, {390.00, 8.9570}, {400.00, 8.9935}, {410.00, 9.0300}, {410.00, 9.3600}, {420.00, 9.3936}, {430.00, 9.4272}, {440.00, 9.4608}, {450.00, 9.4944}, {460.00, 9.5280}, {470.00, 9.5616}, {480.00, 9.5952}, {490.00, 9.6288}, {500.00, 9.6624}, {510.00, 9.6960}, {520.00, 9.7296}, {530.00, 9.7632}, {540.00, 9.7968}, {550.00, 9.8304}, {560.00, 9.8640}, {570.00, 9.8976}, {580.00, 9.9312}, {590.00, 9.9648}, {600.00, 9.9984}, {610.00, 10.0320}, {620.00, 10.0656}, {630.00, 10.0992}, {640.00, 10.1328}, {650.00, 10.1664}, {660.00, 10.2000}, {660.00, 10.7900}, {670.00, 10.8166}, {680.00, 10.8432}, {690.00, 10.8697}, {700.00, 10.8963}, {710.00, 10.9229}, {720.00, 10.9495}, {730.00, 10.9761}, {740.00, 11.0026}, {750.00, 11.0292}, {760.00, 11.0558}, {760.00, 11.0558}, {770.00, 11.0738}, {780.00, 11.0917}, {790.00, 11.1095}, {800.00, 11.1272}, {900.00, 11.2997}, {1000.00, 11.4640}, {1100.00, 11.6208}, {1200.00, 11.7707}, {1300.00, 11.9142}, {1400.00, 12.0521}, {1500.00, 12.1849}, {2000.00, 12.7944}, {2500.00, 13.3697}, {2700.00, 13.6076}, {2740.00, 13.6564}, {2740.00, 13.6564}, {2750.00, 13.6587}, {2800.00, 13.6703}, {2850.00, 13.6818}, {2889.00, 13.6908}, {2889.00, 8.0088}, {2900.00, 8.0280}, {3000.00, 8.1995}, {3100.00, 8.3642}, {3200.00, 8.5222}, {3300.00, 8.6735}, {3400.00, 8.8180}, {3500.00, 8.9558}, {4000.00, 9.5437}, {4500.00, 9.9633}, {5153.90, 10.2578}, {5153.90, 11.0914}, {5500.00, 11.1644}, {6000.00, 11.2270}, {6371.00, 11.2409}}; // fortran test model inline const std::vector> model_1d_prem \ = {{ 0.00, 5.80000}, { 15.00, 5.80000}, { 15.00, 6.80000}, { 24.40, 6.80000}, { 24.40, 8.11061}, { 40.00, 8.10119}, { 60.00, 8.08907}, { 80.00, 8.07688}, { 115.00, 8.05540}, { 150.00, 8.03370}, { 185.00, 8.01180}, { 220.00, 7.98970}, { 220.00, 8.55896}, { 265.00, 8.64552}, { 310.00, 8.73209}, { 355.00, 8.81867}, { 400.00, 8.90522}, { 400.00, 9.13397}, { 450.00, 9.38990}, { 500.00, 9.64588}, { 550.00, 9.90185}, { 600.00, 10.15782}, { 635.00, 10.21203}, { 670.00, 10.26622}, { 670.00, 10.75131}, { 721.00, 10.91005}, { 771.00, 11.06557}, { 871.00, 11.24490}, { 971.00, 11.41560}, { 1071.00, 11.57828}, { 1171.00, 11.73357}, { 1271.00, 11.88209}, { 1371.00, 12.02445}, { 1471.00, 12.16126}, { 1571.00, 12.29316}, { 1671.00, 12.42075}, { 1771.00, 12.54466}, { 1871.00, 12.66550}, { 1971.00, 12.78389}, { 2071.00, 12.90045}, { 2171.00, 13.01579}, { 2271.00, 13.13055}, { 2371.00, 13.24532}, { 2471.00, 13.36074}, { 2571.00, 13.47742}, { 2671.00, 13.59597}, { 2741.00, 13.68041}, { 2771.00, 13.68753}, { 2871.00, 13.71168}, { 2891.00, 13.71660}, { 2891.00, 8.06482}, { 2971.00, 8.19939}, { 3071.00, 8.36019}, { 3171.00, 8.51298}, { 3271.00, 8.65805}, { 3371.00, 8.79573}, { 3471.00, 8.92632}, { 3571.00, 9.05015}, { 3671.00, 9.16752}, { 3771.00, 9.27867}, { 3871.00, 9.38418}, { 3971.00, 9.48409}, { 4071.00, 9.57881}, { 4171.00, 9.66865}, { 4271.00, 9.75393}, { 4371.00, 9.83496}, { 4471.00, 9.91206}, { 4571.00, 9.98554}, { 4671.00, 10.05572}, { 4771.00, 10.12291}, { 4871.00, 10.18743}, { 4971.00, 10.24959}, { 5071.00, 10.30971}, { 5149.50, 10.35568}, { 5149.50, 11.02827}, { 5171.00, 11.03643}, { 5271.00, 11.07249}, { 5371.00, 11.10542}, { 5471.00, 11.13521}, { 5571.00, 11.16186}, { 5671.00, 11.18538}, { 5771.00, 11.20576}, { 5871.00, 11.22301}, { 5971.00, 11.23712}, { 6071.00, 11.24809}, { 6171.00, 11.25593}, { 6271.00, 11.26064}, { 6371.00, 11.26220}}; #endif