aaspip / pyekfmm

A python package for 3D fast-marching-based traveltime calculation and its applications in seismology
GNU General Public License v3.0
47 stars 10 forks source link

Calculating incidence angles to receivers #7

Open abhargavashah opened 1 month ago

abhargavashah commented 1 month ago

Hello!

I had a question regarding the eikonal.c file in pyekfmm: I believe the calculation of the takeoff angles from the source occurs in line 914 (dip[i1+i2*n1+i3*n1*n2+is*n1*n2*n3]=atan2(sqrt(grad1 * grad1 + grad2 * grad2), -grad3) / cRPD;). I was also hoping to be able to calculate the incidence angles to the receivers. Would you by any chance know how best to modify or add to the current code in order to do this?

Thank you! I greatly appreciate any insight you might be able to provide!

-abhargavashah

chenyk1990 commented 1 month ago

Hi, abhargavashah,

Received. I'll help you solve the problem. I just returned from vacation, and work has piled up. I'll return with a solution ASAP.

Thanks, YC

On Thu, Jul 11, 2024 at 4:07 PM abhargavashah @.***> wrote:

Hello!

I had a question regarding the eikonal.c file in pyekfmm: I believe the calculation of the takeoff angles from the source occurs in line 914 (dip[i1+i2n1+i3n1n2+isn1n2n3]=atan2(sqrt(grad1

  • grad1 + grad2 * grad2), -grad3) / cRPD;). I was also hoping to be able to calculate the incidence angles to the receivers. Would you by any chance know how best to modify or add to the current code in order to do this?

Thank you! I greatly appreciate any insight you might be able to provide!

-abhargavashah

— Reply to this email directly, view it on GitHub https://github.com/aaspip/pyekfmm/issues/7, or unsubscribe https://github.com/notifications/unsubscribe-auth/ADHCWESQ3NINYG6PHQX2UULZL3XYZAVCNFSM6AAAAABKXXNZL6VHI2DSMVQWIX3LMV43ASLTON2WKOZSGQYDIMJSG43TAOI . You are receiving this because you are subscribed to this thread.Message ID: @.***>

abhargavashah commented 1 month ago

No worries at all, Yangkang! Thank you so much!

-abhargavashah