vel_inc problem

[emsutti]

Hello Mark,

I'm running HASHpy with python3. I modified some the I/O functions so I matches with the new release. Neverthless I get this error:

*** event outside of velocity table depth range ***

I think that probaly the vel_inc script is not reading the velocity model properly.
Please, give me a hint where can I modify the script.

PD: I've attached my velocity model.
model.txt

[markcwill]

That error message is from the Fortran code which I did not write, it's the original Hardebeck & Shearer code, which has been re-syntaxed to F95 from F77. You will have to check the Fortran source for more info. We use the exact file formats of the original Fortran code, so you may want to compile the original HASH driver in Fortran and try your model, and see if you get an error. If you don't, that would indicate the problem is with this wrapper code, and not the Fortran libs.

But at a quick glance, you have a lot of layers, and your depths are outside the range of the array limits and sizes. You can look at the default values in the Fortran include files, e.g. "vel.inc". If you want to change them, you would need to edit all of those values and re-compile the library & wrappers.

[emsutti]

Hi Mark!

I figured it out finally! Thanks a lot for your help!

Best!

Emilio.