GenericMappingTools / GenericMappingTools/gmt

`pssac` produces different polarity by Y-size in J-option

Open
#6,448 3 comments 0 reactions 0 assignees View on GitHub
bug
Dominant language
C
Stars
979
Forks
414
Avg merge
17h 26m
Merged PRs (30d)
54

Description

I have updated `gmt6` recently and found a very useful module `pssac` (`gmt sac` as a modern name). It reads sac files directly and make a record section while `pswiggle` needs an additional conversion from the sac file to xy-file. But, I like to notice that a feature of `pssac` can mislead `gmt` users. `pssac` represents the seismograms with different polarity when Y-size in J-option is positive and negative.

For example, "-J10c/5c" and "-J10c/-5c" make two reversed plots. This behavior is not natural because one expects that the upward direction of the paper is positive, in my opinion.

Let us make a sac file.

```
sac < a.ps
```
![a](https://user-images.githubusercontent.com/82625620/158328826-761bfddf-3b8d-4624-a670-7e54f8e8d1a1.png)
*Normal polarity (a.ps)*

```
gmt pssac seismogram.sac -R10/12/-1/1 -JX10c/-10c -P -M1 -Ba -Xc -Yc > b.ps
```
![b](https://user-images.githubusercontent.com/82625620/158328837-8140c832-504c-439b-bfec-e3400ddec603.png)
*Reversed polarity with the negative Y-size in J-option (b.ps)*

The similar module `pswiggle` does not behave like this way.

**Suggestion**
The polarity of the seismograms should be fixed as positive upward in the paper.

Contributor guide

Open the contributing guide

Assessment

This issue has not been assessed yet.

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.