-
Notifications
You must be signed in to change notification settings - Fork 22
-
Notifications
You must be signed in to change notification settings - Fork 22
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Handle multiple sources and receivers with varying orientations in 'bipole' #74
Comments
Thanks for reporting @AndyMcAliley , this seems to be a bug. I'll look into it! |
I think # Remove if source is x- or y-directed
check = np.atleast_1d(srcazm)[0]
if np.allclose(srcazm % (np.pi/2), 0): # if all angles are multiples of 90
if np.isclose(check // (np.pi/2) % 2, 0): # Multiples of pi (180)
bab[:, 1] *= False # x-directed source, remove y
else: # Multiples of pi/2 (90)
bab[:, 0] *= False # y-directed source, remove x |
Yes, you are right. You touched a scenario that I never used so far...! I fixed it and will soon push an update. |
And it obviously got not caught by the tests, so I'll have to add a test too... |
OK, a fix is in the branch azimuth-dip. Could you test it ( If it works I will push later this week a new release. |
Looks good to me, for both sources and receivers! Thanks so much! |
In
empymod.bipole
, whensrc
orrec
is formatted as[x0, x1, y0, y1, z0, z1]
with multiple bipoles with different azimuths, the fields are computed as if all bipoles have the same orientation as the first in the list. For example, for two orthogonal sources:All entries of
bex[:,1]
are0+0j
. Similarly, for two orthogonal receivers,Output:
0j
.The text was updated successfully, but these errors were encountered: