"""
# `EK137`

1-D layered model of the Earth, containing density.

## Reference

Kennett, B.L.N., 2020. Radial earth models revisited.
Geophysical Journal International 222, 2189–2204.
https://doi.org/10.1093/gji/ggaa298
"""
const EK137 = LinearLayeredModel(
    r = 6371. .- reverse([0.000, 20.000, 20.000, 35.000, 35.000, 77.500, 120.000, 165.000,
        210.000, 210.000, 260.000, 310.000, 360.000, 410.000, 410.000, 460.000,
        510.000, 560.000, 610.000, 660.000, 660.000, 710.000, 760.000, 809.500,
        859.000, 908.500, 958.000, 1007.500, 1057.000, 1106.500, 1156.000, 1205.500,
        1255.000, 1304.500, 1354.000, 1403.500, 1453.000, 1502.500, 1552.000, 1601.500,
        1651.000, 1700.500, 1750.000, 1799.500, 1849.000, 1898.500, 1948.000, 1997.500,
        2047.000, 2096.500, 2146.000, 2195.500, 2245.000, 2294.500, 2344.000, 2393.500,
        2443.000, 2492.500, 2542.000, 2591.500, 2640.000, 2690.000, 2740.000, 2740.000,
        2789.670, 2839.330, 2891.500, 2891.500, 2939.330, 2989.660, 3039.990, 3090.320,
        3140.660, 3190.990, 3241.320, 3291.650, 3341.980, 3392.310, 3442.640, 3492.970,
        3543.300, 3593.640, 3643.970, 3694.300, 3744.630, 3794.960, 3845.290, 3895.620,
        3945.950, 3996.280, 4046.620, 4096.950, 4147.280, 4197.610, 4247.940, 4298.270,
        4348.600, 4398.930, 4449.260, 4499.600, 4549.930, 4600.260, 4650.590, 4700.920,
        4751.250, 4801.580, 4851.910, 4902.240, 4952.580, 5002.910, 5053.240, 5103.570,
        5153.500, 5153.500, 5204.610, 5255.320, 5306.040, 5356.750, 5407.460, 5458.170,
        5508.890, 5559.600, 5610.310, 5661.020, 5711.740, 5762.450, 5813.160, 5863.870,
        5914.590, 5965.300, 6016.010, 6066.720, 6117.440, 6168.150, 6218.860, 6269.570,
        6320.290, 6371.000]),

    vp = reverse([5.8000, 5.8000, 6.5000, 6.5000, 8.0400, 8.0450, 8.0500, 8.1750,
        8.3000, 8.3000, 8.4825, 8.6640, 8.8445, 9.0200, 9.3800, 9.5320,
        9.6960, 9.8640, 10.0320, 10.2000, 10.7900, 10.9229, 11.0558, 11.1353,
        11.2221, 11.3068, 11.3896, 11.4705, 11.5499, 11.6277, 11.7038, 11.7782,
        11.8511, 11.9224, 11.9923, 12.0609, 12.1281, 12.1952, 12.2586, 12.3217,
        12.3847, 12.4450, 12.5051, 12.5647, 12.6233, 12.6812, 12.7386, 12.7956,
        12.8526, 12.9086, 12.9648, 13.0192, 13.0743, 13.1286, 13.1834, 13.2395,
        13.2928, 13.34851, 13.4056, 13.4641, 13.5212, 13.5800, 13.6394, 13.6394,
        13.6410, 13.6426, 13.64421, 8.0090, 8.0941, 8.1816, 8.2666, 8.3496,
        8.4305, 8.5084, 8.5848, 8.6591, 8.7320, 8.8024, 8.8714, 8.9382,
        9.0034, 9.0664, 9.1286, 9.1889, 9.2470, 9.3047, 9.3604, 9.4145,
        9.4669, 9.5180, 9.5677, 9.6159, 9.6623, 9.7080, 9.7513, 9.7914,
        9.8304, 9.8682, 9.9051, 9.9410, 9.9761, 10.0103, 10.0439, 10.0768,
        10.1092, 10.1415, 10.1739, 10.2049, 10.2329, 10.2565, 10.2745, 10.2854,
        10.2890, 11.0427, 11.0585, 11.0718, 11.0850, 11.0983, 11.1166, 11.1316,
        11.1457, 11.1590, 11.1715, 11.1832, 11.1941, 11.2041, 11.2134, 11.2219,
        11.2295, 11.2364, 11.2424, 11.2477, 11.2521, 11.2557, 11.2586, 11.2606,
        11.2618, 11.2622]),

    vs = reverse([3.4200, 3.4200, 3.8100, 3.8100, 4.4800, 4.4900, 4.5000, 4.5090,
        4.5180, 4.5230, 4.6090, 4.6960, 4.7880, 4.8600, 5.0900, 5.1910,
        5.2920, 5.3980, 5.5040, 5.6100, 6.0100, 6.1097, 6.2095, 6.2426,
        6.2798, 6.3160, 6.3512, 6.3854, 6.4187, 6.4510, 6.4828, 6.5138,
        6.5439, 6.5727, 6.6008, 6.6285, 6.6555, 6.6815, 6.7073, 6.7326,
        6.7573, 6.7815, 6.8052, 6.8286, 6.8515, 6.8742, 6.8972, 6.9194,
        6.9418, 6.9627, 6.9855, 7.0063, 7.0281, 7.0490, 7.0700, 7.0901,
        7.1104, 7.1319, 7.1526, 7.1737, 7.1951, 7.2168, 7.2390, 7.2390,
        7.2497, 7.2594, 7.2711, 0.0000, 0.0000, 0.0000, 0.0000, 0.0000,
        0.0000, 0.0000, 0.0000, 0.0000, 0.0000, 0.0000, 0.0000, 0.0000,
        0.0000, 0.0000, 0.0000, 0.0000, 0.0000, 0.0000, 0.0000, 0.0000,
        0.0000, 0.0000, 0.0000, 0.0000, 0.0000, 0.0000, 0.0000, 0.0000,
        0.0000, 0.0000, 0.0000, 0.0000, 0.0000, 0.0000, 0.0000, 0.0000,
        0.0000, 0.0000, 0.0000, 0.0000, 0.0000, 0.0000, 0.0000, 0.0000,
        0.0000, 3.4243, 3.4387, 3.4514, 3.4635, 3.4751, 3.4861, 3.4965,
        3.5064, 3.5157, 3.5244, 3.5326, 3.5402, 3.5472, 3.5537, 3.5596,
        3.5650, 3.5698, 3.5740, 3.5777, 3.5808, 3.5833, 3.5867, 3.5875,
        3.5878, 3.5878]),

    density = reverse([2.6510, 2.6510, 2.8630, 2.8630, 3.2900, 3.2940, 3.2960, 3.3425,
        3.3889, 3.3889, 3.4512, 3.5203, 3.5894, 3.6586, 3.7836, 3.8482,
        3.9108, 3.9735, 4.0371, 4.1008, 4.3253, 4.3756, 4.4259, 4.4561,
        4.4899, 4.5220, 4.5654, 4.5926, 4.6198, 4.6467, 4.6736, 4.7001,
        4.7266, 4.7528, 4.7790, 4.8050, 4.8307, 4.8562, 4.8817, 4.9069,
        4.9321, 4.9570, 4.9817, 5.0062, 5.0306, 5.0548, 5.0789, 5.1027,
        5.1264, 5.1499, 5.1732, 5.1963, 5.2192, 5.2420, 5.2647, 5.2870,
        5.3092, 5.3313, 5.3532, 5.3748, 5.4012, 5.4276, 5.4537, 5.6114,
        5.6476, 5.6828, 5.7591, 9.9145, 9.9942, 10.0722, 10.1485, 10.2233,
        10.2964, 10.3679, 10.4378, 10.5062, 10.5731, 10.6385, 10.7023, 10.7647,
        10.8257, 10.8852, 10.9434, 11.0001, 11.0555, 11.1095, 11.1623, 11.2137,
        11.2639, 11.3127, 11.3604, 11.4069, 11.4521, 11.4962, 11.5391, 11.5809,
        11.6216, 11.6612, 11.6998, 11.7373, 11.7737, 11.8092, 11.8437, 11.8772,
        11.9098, 11.9414, 11.9722, 12.0001, 12.0311, 12.0593, 12.0867, 12.1133,
        12.1391, 12.7037, 12.7289, 12.7530, 12.7760, 12.7980, 12.8188, 12.8387,
        12.8574, 12.8751, 12.8917, 12.9072, 12.9217, 12.9351, 12.9474, 12.9586,
        12.9688, 12.9779, 12.9859, 12.9929, 12.9988, 13.0036, 13.0074, 13.0100,
        13.0117, 13.0122])
    )
