search

trichter / rf

Receiver function calculation in seismology
MIT License
106 stars 62 forks source link

moveout for Sp #2

Closed tbelinic closed 8 years ago

tbelinic commented 8 years ago

Great job, but moveout for Sp doesn't work properly:

In [177]: new_t = np.interp(t, t0, t1, left=0, right=0) In [178]: new_t Out[178]: array([ 0., 0., 0., ..., 0., 0., 0.])

For Sp you have to change right to None: new_t = np.interp(t, t0, t1, left=0, right=None)

trichter commented 8 years ago

Thanks for the bug report and the suggested solution. I'll have a look into it.

@qwertz07 I cannot reproduce the error. Could you please provide a small working (or better failing) example.

Closing for now. Please reopen if this is still applicable.

Edit: Finally, I got the point. This is now fixed. Thanks again.