search

GenericMappingTools / gmt

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

psmeca bug, -T option #8053

Open jpesicek opened 10 months ago

jpesicek commented 10 months ago

Description of the problem

when plotting focal mechanisms with psmeca and the -T option, nodal planes sometimes plot incompletely and/or outside the circle. One example below:

Full script that generated the error

echo 2 2 0 151.71 75.00 180.00 5 | gmt psmeca -JX10 -R0/4/0/4 -Sa10 -N -: -L -T > plot.ps

Actual outcome

line plots outside circle

Expected outcome

if you remove the -T option, the plot behaves as expected

System information

PaulWessel commented 10 months ago

Confirmed. Looks like a problem with exactly 180. Both 179.xxx and 180.xxx seems to work. @seisman, any thoughts on why 180 "blows up" somewhat? The French code must do something like divide by zero (sin (180))?

seisman commented 10 months ago

Need to find some time to debug it.

jpesicek commented 10 months ago

one more:

echo 2 2 0 80.00 90.00 -140.00 5 | gmt psmeca -JX10 -R0/4/0/4 -Sa10 -N -: -L -T > plot2.ps

PaulWessel commented 10 months ago

This time it is the "90.00" that cause trouble, it seems.

PaulWessel commented 8 months ago

Ping @seisman on this again? Gotta be something dumb like sin(180) and then divide by it?

seisman commented 8 months ago

OK, will look at it again next week.