search

GenericMappingTools / gmt

The Generic Mapping Tools
https://www.generic-mapping-tools.org
Other
843 stars 351 forks source link

psmeca dump option #80

Open 79seismo opened 6 years ago

79seismo commented 6 years ago

When you are using psmeca with -Sa flag (Aki and Richards convention where you specify only one fault plane), write out results of the auxiliary plane and P and T axes. These are already plotted but there is no option for outputing them.

welcome[bot] commented 6 years ago

👋 Thanks for opening your first issue here! Please make sure you filled out the template with as much detail as possible. We appreciate that you took the time to contribute!

seisman commented 5 years ago

Although it's straightforward to output focal mechanism information, I think it's better to keep it as a simple plotting tool rather than a focal mechanism calculator.

Feel free to open it if someone still would like to add such feature to psmeca.

79seismo commented 5 years ago

Not sure why it's hard to output the other nodal plane because it's already being calculated within GMT.

seisman commented 5 years ago

No, it's not hard. The problem is, nodal plane information is just a subset of full focal mechanism. If we output nodal plane, why not also output T, N, P axis information, full moment tensor or some other related information.

BTW, the mechanism codes in GMT is now buggy and not fully tested.

79seismo commented 5 years ago

I've used my own code for calculations, and so far, GMT output has been consistent with that.

joa-quim commented 5 years ago

I see output P,T,N as an interesting feature. But it should be made with the IO mechanism to allow sending to external interfaces, not simple fprintf()