I searched through the source code on GitHub (master) and it looks the team has not yet got around to implementing a method for psmeca in base_plotting.py. Am I correct? Or is there another way to draw earthquake focal mechanisms using PyGMT?
If not, it looks like I could try to implement it myself following the examples you have calling other modules (e.g. coast).
Looks like I was right, you guys did a really great job of making the code easily modular. I was able to (at least partially) implement the meca feature. I did this quickly and committed this back to my forked version of the master branch. Have not submitted a pull request because it probably violates your contributers code of conduct - haven’t had the time to review it thoroughly yet.
fig = pygmt.Figure()
fig.basemap(region=region, frame=“af”, projection=“M0/0/8i”)
fig.coast(shorelines=True, N=2, land=‘gray’, water=‘skyblue’)
fig.meca(table=“MTTable.txt”, Sd=“1”, C=" ")