Welcome to SourceSpec Discussions!

[claudiodsf]

👋 Welcome!

We’re using Discussions as a place to connect with other members of our community. We hope that you:

• Ask questions you’re wondering about.
• Share ideas.
• Engage with other community members.
• Welcome others and are open-minded. Remember that this is a community we
  build together 💪.

[anniejerk]

Hi,
I’d like to try using a different velocity model instead of iasp91 to compute the takeoff angles when incorporating the radiation pattern and focal mechanism in the analysis. Is this currently possible?

[anniejerk]

@claudiodsf I've been working with NonLinLoc and have generated some grids, but since this is my first time using it, I'm having a bit of trouble getting your code to run smoothly with these grids.

When I attempt to run it, I encounter the following error messages:

input must be a scalar
NS.BER.00.HHE: Unable to get P arrival time: skipping trace
input must be a scalar
NS.BER.00.HHN: Unable to ...
The arrivals seemed to work previously, as they are stored in a separate .xml file, so I'm not quite sure what's going wrong now.

Would you happen to have any suggestions for what might be causing this? Also, is there a way to increase the verbosity level of the output so I can get a better idea of where the issue is occurring?

Thanks so much in advance for your help—I really appreciate your time and expertise!

[claudiodsf]

Hi @anniejerk. I don't have enough clues to help you. Could you try sharing somewhere a complete working example, with input files and your configuration?

[anniejerk]

sourcespec (3).zip

[claudiodsf]

Thanks,

I could run your test and found two problems:

1. A bug in SouceSpec preventing the grids to being correctly read. This is fixed in the latest development version, which you can install through pip install git+https://github.com/SeismicSource/sourcespec.git

2. Your grids are not correctly constructed: the epicenter is on the (0, 0) position, which is the lower-left corner of the grid. Since your model is 1D, you could construct a simpler 2D grid with something like VGGRID 2 250 10 0.0 0.0 -5.0 10.0 10.0 10.0 SLOW_LEN (note the first value of "2" for the number of x points). You should then use GTMODE GRID2D ANGLES_YES

[anniejerk]

Thank you so much!