Skip to content

Commit

Permalink
test for moon modle with LVZ in mid mantle
Browse files Browse the repository at this point in the history
  • Loading branch information
@crotwell
crotwell committed Oct 28, 2024
1 parent 77a52cd commit 5711a8d
Show file tree
Hide file tree
Showing 2 changed files with 93 additions and 0 deletions.
17 changes: 17 additions & 0 deletions src/test/java/edu/sc/seis/TauP/MoonLowVelocityTest.java
View file
Original file line number Diff line number Diff line change
Expand Up @@ -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 {
Expand All @@ -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<Arrival> arrivalList = dr.calculate(P_phase);
assertEquals(1, arrivalList.size());
assertEquals(521.75, arrivalList.get(0).getTime(), 0.01);
}
}
76 changes: 76 additions & 0 deletions src/test/resources/edu/sc/seis/TauP/MoonKhan_JGR_2014_mod1.nd
View file
Original file line number Diff line number Diff line change
@@ -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

0 comments on commit 5711a8d

Please sign in to comment.