Bug #894

wrong beach ball when rake=0

Added by Dongdong over 5 years ago. Updated over 5 years ago.

Status:ClosedStart date:2016-04-27
Priority:NormalDue date:
Assignee:Paul% Done:


Target version:Candidate for next bugfix release
Affected version:5.x-svn Platform:


The mechanism I am trying to plot is from GCMT (http://www.globalcmt.org/CMTsearch.html). Since the url link is too long, I just paste the screenshot.

Commands below can reproduce the bug. This command plots three beach balls with slightly different rake angles.

gmt psmeca -JX15c/7c -R116/127/23/26 -B1 -Sa2c -Fa0.2c/cc > meca.ps << EOF
118.68 24.39 15.0 351 81 -0.001 5.4 0 0 rake=-0.001
120.68 24.39 15.0 351 81 0      5.4 0 0 rake=0
122.68 24.39 15.0 351 81 0.001  5.4 0 0 rake=0.001

The output is:

When rake is exactly equal to 0, the beach ball is wrong. So, I think this may be a bug of psmeca.

PS: I am using GMT 5.2.1 under Linux x86_64.

GCMT.png (29.7 KB) Dongdong, 2016-04-27 22:52

meca.jpg (33.7 KB) Dongdong, 2016-04-28 01:14


#1 Updated by Joaquim over 5 years ago

  • Status changed from New to Resolved

Ok, in r16346 I made the hacky fix of setting
if (rake == 0) rake = 0.00001;

Not ideal but should work and given this code was not written by us a perfect fix would risk to take real long time.

#2 Updated by Paul over 5 years ago

I have added new test script meca_13.sh to remind us. This now passes after JL's fix.

#3 Updated by Paul over 5 years ago

  • Status changed from Resolved to Closed
  • Assignee set to Paul
  • % Done changed from 0 to 100

Closed as fixed - reopen if there are other problems related to this.

#4 Updated by Joaquim over 5 years ago

Paul, but it was you who said that some of our meca tests were failing after this fix. No more?

#5 Updated by Paul over 5 years ago

I have no more failures after reconciling the various troubles, so I am closing it for that reason.

Also available in: Atom PDF