psmeca bug, -T option

[jpesicek]

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

• Operating system: macOS
• GMT version (gmt --version): 6.4.0

[PaulWessel]

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]

Need to find some time to debug it.

[jpesicek]

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]

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

[PaulWessel]

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

[seisman]

OK, will look at it again next week.

[jpesicek]

Just checking in on the status of this bug, thanks

[jpesicek]

another note that when I use gmt4, the lines are plotted correctly