I'm writing this post anxiously awaiting my five hour code to finish running. As I've mentioned previously, the tomography code is not yet working for my project. That being said, I've made what is hopefully a huge step forward this morning. I guess moving to Purdue provided a lucky spark!
The essence of the problem was that our geometry was not working out correctly. Because we are dealing with P-wave tomography, one of the complications is that we need all of our waves to enter through the bottom of our model, not the sides. It's not a straightforward calculation to compute this, and you can't just overestimate because the larger your grid becomes, the slower your program runs. Without this geometry behaving correctly, the travel times for the stations are computed incorrectly.
To visually display this information, Dr. Pavlis wrote some code that turns the P-Wave arrival times files into a format that is viewable by Paraview, a 3D visualization software. When the calculations are correct, the arrivals are supposed to be sphere/ellipsoid in shape. When they are not correct, well I'll just let you take a look...
Needless to say, not spherical in shape. However, after making a few adjustments to the buffer zone on the depth grid geometry, I have images looking like this:
Spheres! Not stair step messy jumbles! Now, this is just the travel times for one specific station. I decided to test just one station this time before running the four+ hour code that computes every station, which I made the mistake of doing last time. Now that I have this image however, the code for all the stations is currently running. After this is finished, I'll try the inversion code one more time. If it works this time, I may actually do a little victory dance. This is potentially a huge huge step towards completing tomography model version 1.0. Needless to say, I'm a little excited. Because apparently clean geometric shapes do that to you when you are a scientist.
Hopefully I'll be reporting back later this week with even more successes. C'mon Purdue, I know you have more luck in you!
You must be logged into the CMS to post a comment.