diff --git a/src/test/java/edu/sc/seis/TauP/MoonLowVelocityTest.java b/src/test/java/edu/sc/seis/TauP/MoonLowVelocityTest.java index f8b0aadc..8184e45e 100644 --- a/src/test/java/edu/sc/seis/TauP/MoonLowVelocityTest.java +++ b/src/test/java/edu/sc/seis/TauP/MoonLowVelocityTest.java @@ -2,6 +2,10 @@ import org.junit.jupiter.api.Test; +import java.io.IOException; +import java.util.List; + +import static org.junit.jupiter.api.Assertions.assertEquals; import static org.junit.jupiter.api.Assertions.assertNotNull; public class MoonLowVelocityTest { @@ -24,4 +28,17 @@ public void testTMod() throws Exception { TauModel tmoddepth = tmod.depthCorrect(933); assertNotNull(tmoddepth); } + + @Test + public void testMoonKhan_JGR_2014() throws SlownessModelException, VelocityModelException, IOException, TauModelException { + String lunarModelFile = "MoonKhan_JGR_2014_mod1.nd"; + VelocityModel vmod = VelocityModelTest.loadTestVelMod(lunarModelFile); + SlownessModel smod = new SphericalSModel(vmod); + TauModel tmod = new TauModel(smod); + SeismicPhase P_phase = SeismicPhaseFactory.createSinglePhase("P", tmod, 800); + DistanceRay dr = DistanceRay.ofDegrees(50); + List arrivalList = dr.calculate(P_phase); + assertEquals(1, arrivalList.size()); + assertEquals(521.75, arrivalList.get(0).getTime(), 0.01); + } } diff --git a/src/test/resources/edu/sc/seis/TauP/MoonKhan_JGR_2014_mod1.nd b/src/test/resources/edu/sc/seis/TauP/MoonKhan_JGR_2014_mod1.nd new file mode 100755 index 00000000..b83720c5 --- /dev/null +++ b/src/test/resources/edu/sc/seis/TauP/MoonKhan_JGR_2014_mod1.nd @@ -0,0 +1,76 @@ + 0.000 4.409 2.414 2.602 6750 6750 + 20.000 5.670 3.119 2.763 6750 6750 + 40.000 6.742 3.713 2.889 6750 6750 +mantle + 40.000 7.652 4.325 3.234 4000 9000 + 60.000 7.652 4.325 3.234 4000 9000 + 80.000 7.657 4.327 3.236 4000 9000 + 100.000 7.662 4.329 3.239 4000 9000 + 120.000 7.822 4.466 3.352 4000 9000 + 140.000 7.825 4.465 3.353 4000 9000 + 160.000 7.827 4.464 3.355 4000 9000 + 180.000 7.829 4.463 3.356 4000 9000 + 200.000 7.829 4.461 3.357 4000 9000 + 220.000 7.829 4.459 3.358 4000 9000 + 240.000 7.872 4.476 3.371 4000 9000 + 260.000 7.978 4.521 3.401 1500 3375 + 280.000 8.036 4.544 3.417 1500 3375 + 300.000 8.040 4.543 3.419 1500 3375 + 320.000 8.042 4.543 3.420 1500 3375 + 340.000 8.044 4.542 3.421 1500 3375 + 360.000 8.045 4.541 3.422 1500 3375 + 380.000 8.046 4.539 3.423 1500 3375 + 400.000 8.047 4.538 3.424 1500 3375 + 420.000 8.048 4.537 3.425 1500 3375 + 440.000 8.049 4.536 3.426 1500 3375 + 460.000 8.051 4.535 3.427 1500 3375 + 480.000 8.053 4.535 3.428 1500 3375 + 500.000 8.054 4.534 3.429 1500 3375 + 520.000 8.055 4.533 3.430 1500 3375 + 540.000 8.056 4.532 3.431 500 1125 + 560.000 8.057 4.531 3.432 500 1125 + 580.000 8.057 4.530 3.433 500 1125 + 600.000 8.057 4.529 3.433 500 1125 + 620.000 8.057 4.528 3.434 500 1125 + 640.000 8.057 4.526 3.435 500 1125 + 660.000 8.056 4.524 3.435 500 1125 + 680.000 8.055 4.523 3.436 500 1125 + 700.000 8.054 4.521 3.436 500 1125 + 720.000 8.052 4.519 3.436 500 1125 + 740.000 8.050 4.516 3.437 500 1125 + 760.000 8.048 4.514 3.437 300 675 + 780.000 8.047 4.512 3.437 300 675 + 800.000 8.045 4.510 3.437 300 675 + 820.000 8.043 4.508 3.438 300 675 + 840.000 8.042 4.506 3.438 300 675 + 860.000 8.040 4.504 3.438 300 675 + 880.000 8.038 4.502 3.438 300 675 + 900.000 8.037 4.500 3.439 300 675 + 920.000 8.036 4.499 3.439 300 675 + 940.000 8.035 4.498 3.439 300 675 + 960.000 8.034 4.496 3.440 300 675 + 980.000 8.031 4.494 3.440 300 675 + 1000.000 8.030 4.493 3.440 300 675 + 1020.000 8.028 4.491 3.440 300 675 + 1040.000 8.026 4.489 3.440 300 675 + 1060.000 8.024 4.487 3.440 300 675 + 1080.000 8.021 4.485 3.440 300 675 + 1100.000 8.016 4.482 3.439 300 675 + 1120.000 8.002 4.476 3.435 300 675 + 1140.000 7.982 4.469 3.428 300 675 + 1160.000 7.944 4.453 3.416 300 675 + 1180.000 7.883 4.437 3.395 300 675 + 1200.000 7.813 4.410 3.374 300 675 + 1220.000 7.664 4.302 3.353 300 675 + 1240.000 7.358 4.991 3.338 300 675 + 1240.000 7.709 5.335 3.317 150 500 + 1260.000 7.709 5.335 3.317 150 500 + 1280.000 8.107 5.747 3.344 150 500 + 1300.000 8.526 6.213 3.358 150 500 + 1407.415 8.526 6.213 3.368 150 500 +outer-core + 1407.415 2.622 0.000 4.495 0 10000 + 1636.000 2.622 0.000 4.495 0 10000 +inner-core + 1636.000 7.622 4.000 4.495 0 10000 + 1737.000 7.622 4.000 4.495 0 10000