earthquake

[meiabdalla]

when you use earthquake deformation you mention in the title "is magnitude correct"?
is the deformation correct because if I use another okada Wrapper it gives another pattern of deformation?
I use (https://github.com/tbenthompson/okada_wrapper)

[meiabdalla]

SyInterferoPy/lib/syinterferopy_functions.py

Line 240 in 627037b

centroid_depth = np.mean((kwargs['bottom_depth'] - kwargs['top_depth']))

I think this line should be
centroid_depth = np.mean((kwargs['bottom_depth'], kwargs['top_depth']))

[matthew-gaddes]

Hi there, Thanks for suggesting this, and sorry for the slow reply. I developed this code to primarily work with volcanoes and included the Okada wrapper for sills and dykes. It was easy to extend this to earthquakes, but when I did this I did not think the magnitudes were correct, hence the comment.

I have tried the fix you suggested and, whilst this improves things, I still don't think the magnitudes are correct (though they are closer). I will look at it in more detail next week.

Thanks, Matt

[meiabdalla]

thank you for your replay, I checked the results from your code with the results from https://github.com/geniusinaction/okapy/blob/master/okada_ascending_descending.ipynb with okada wrapper . The results were the same except for the different epsilon ratios to the two codes. the epsilon ratio in okada wrapper is alpha = (lmda + mu) / (lmda + 2 * mu) # elastic constant used by Okada, which is strange to me.