.\" Man page generated from reStructuredText.
.
.TH "PSSAC" "1gmt" "May 21, 2019" "5.4.5" "GMT"
.SH NAME
pssac \- Plot seismograms in SAC format on maps
.
.nr rst2man-indent-level 0
.
.de1 rstReportMargin
\\$1 \\n[an-margin]
level \\n[rst2man-indent-level]
level margin: \\n[rst2man-indent\\n[rst2man-indent-level]]
-
\\n[rst2man-indent0]
\\n[rst2man-indent1]
\\n[rst2man-indent2]
..
.de1 INDENT
.\" .rstReportMargin pre:
. RS \\$1
. nr rst2man-indent\\n[rst2man-indent-level] \\n[an-margin]
. nr rst2man-indent-level +1
.\" .rstReportMargin post:
..
.de UNINDENT
. RE
.\" indent \\n[an-margin]
.\" old: \\n[rst2man-indent\\n[rst2man-indent-level]]
.nr rst2man-indent-level -1
.\" new: \\n[rst2man-indent\\n[rst2man-indent-level]]
.in \\n[rst2man-indent\\n[rst2man-indent-level]]u
..
.SH SYNOPSIS
.sp
\fBpssac\fP [ \fIsaclist\fP|\fISACfiles\fP ]  \fB\-J\fP\fIparameters\fP
 \fB\-R\fP\fIregion\fP
[  \fB\-B\fP[\fBp\fP|\fBs\fP]\fIparameters\fP ]
[  \fB\-C\fP[\fIt0/t1\fP] ]
[  \fB\-D\fP\fIdx\fP[/\fIdy\fP] ]
[  \fB\-E\fP\fBa\fP|\fBb\fP|\fBk\fP|\fBd\fP|\fBn\fP[\fIn\fP]|\fBu\fP[\fIn\fP] ]
[  \fB\-F\fP[\fBi\fP][\fBq\fP][\fBr\fP] ]
[  \fB\-G\fP[\fBp\fP|\fBn\fP][\fB+g\fP\fIfill\fP][\fB+z\fP\fIzero\fP][\fB+t\fP\fIt0/t1\fP] ]
[  \fB\-K\fP ]
[  \fB\-M\fP\fIsize\fP[\fIu\fP][/\fIalpha\fP] ]
[  \fB\-O\fP ]
[  \fB\-P\fP ]
[  \fB\-Q\fP ]
[  \fB\-S\fP[\fBi\fP]\fIscale\fP[\fIunit\fP] ]
[  \fB\-T\fP[\fB+t\fP\fIn\fP][\fB+r\fP\fIreduce_vel\fP][\fB+s\fP\fIshift\fP] ]
[  \fB\-U\fP[\fIstamp\fP] ]
[  \fB\-V\fP[\fIlevel\fP] ]
[  \fB\-W\fP\fIpen\fP ]
[  \fB\-X\fP\fIx_offset\fP ]
[  \fB\-Y\fP\fIy_offset\fP ]
[ \fB\-h\fPheaders ]
[ \fB\-t\fPtransp ]
.sp
\fBNote:\fP No space is allowed between the option flag and the associated arguments.
.SH DESCRIPTION
.sp
\fBpssac\fP reads \fISACfiles\fP in SAC format or reads filenames and controlling parameters
from \fIsaclist\fP [or standard input] and generates PostScript that will plot seismograms on a map.
The PostScript code is written to standard output.
.SH REQUIRED ARGUMENTS
.INDENT 0.0
.TP
.B \fISACfiles\fP
SAC files to plot on a map. Only evenly spaced SAC data is supported.
.TP
.B \fIsaclist\fP
One ASCII data table file holding a number of data columns. If \fIsaclist\fP is not given then we read from standard input.
Parameters are expected to be in the following columns:
.INDENT 7.0
.INDENT 3.5
\fIfilename\fP [\fIX\fP \fIY\fP [\fIpen\fP]]
.UNINDENT
.UNINDENT
.sp
\fIfilename\fP is the name of SAC file to plot.
\fIX\fP and \fIY\fP are the position of seismograms to plot on a map.
On linear plots, the default \fIX\fP is the begin time of SAC file, which will be adjusted if \fB\-T\fP option is used,
the default \fIY\fP is determined by \fB\-E\fP option.
On geographic plots, the default \fIX\fP and \fIY\fP are station longitude and latitude specified in SAC header.
The \fIX\fP and \fIY\fP given here will override the position determined by command line options.
\fIpen\fP, if given, will override the pen from \fB\-W\fP option for current SAC file only.
.UNINDENT
.INDENT 0.0
.TP
\fB\-J\fP\fIparameters\fP (more ...)
Select map projection.  
.UNINDENT
.INDENT 0.0
.TP
\fB\-R\fP\fIwest\fP/\fIeast\fP/\fIsouth\fP/\fInorth\fP[/\fIzmin\fP/\fIzmax\fP][\fB+r\fP][\fB+u\fP\fIunit\fP]
\fIwest\fP, \fIeast\fP, \fIsouth\fP, and \fInorth\fP specify the region of interest,
and you may specify them in decimal degrees or in
[±]dd:mm[:ss.xxx][\fBW\fP|\fBE\fP|\fBS\fP|\fBN\fP] format
Append \fB+r\fP if lower left
and upper right map coordinates are given instead of w/e/s/n. The
two shorthands \fB\-Rg\fP and \fB\-Rd\fP stand for global domain (0/360
and \-180/+180 in longitude respectively, with \-90/+90 in latitude).
Alternatively for grid creation, give \fBR\fP\fIcode\fP\fIlon\fP/\fIlat\fP/\fInx\fP/\fIny\fP, where
\fIcode\fP is a 2\-character combination of L, C, R (for left, center, or right)
and T, M, B for top, middle, or bottom. e.g., BL for lower left.  This
indicates which point on a rectangular region the \fIlon\fP/\fIlat\fP coordinate
refers to, and the grid dimensions \fInx\fP and \fIny\fP with grid spacings via
\fB\-I\fP is used to create the corresponding region.
Alternatively, specify the name of an existing grid file and the
\fB\-R\fP settings (and grid spacing, if applicable) are copied from
the grid. Appending \fB+u\fP\fIunit\fP expects projected (Cartesian)
coordinates compatible with chosen \fB\-J\fP and we inversely project
to determine actual rectangular geographic region.
For perspective view (\fB\-p\fP), optionally append /\fIzmin\fP/\fIzmax\fP\&.
In case of perspective view (\fB\-p\fP), a z\-range (\fIzmin\fP, \fIzmax\fP)
can be appended to indicate the third dimension. This needs to be
done only when using the \fB\-Jz\fP option, not when using only the
\fB\-p\fP option. In the latter case a perspective view of the plane is
plotted, with no third dimension.  
.UNINDENT
.SH OPTIONAL ARGUMENTS
.INDENT 0.0
.TP
\fB\-B\fP[\fBp\fP|\fBs\fP]\fIparameters\fP (more ...)
Set map boundary frame and axes attributes.
.UNINDENT
.INDENT 0.0
.TP
\fB\-C\fP[\fIt0/t1\fP]
Read and plot seismograms in timewindow between \fIt0\fP and \fIt1\fP only.
\fIt0\fP and \fIt1\fP are relative to a reference time specified by \fB\-T\fP\&.
If \fB\-T\fP option is not specified, use the reference time (kzdate and kztime) defined in SAC header instead.
If only \fB\-C\fP is used, \fIt0/t1\fP is determined as \fIxmin/xmax\fP from \fB\-R\fP option.
.UNINDENT
.INDENT 0.0
.TP
\fB\-D\fP\fIdx\fP[\fB/\fP\fIdy\fP]
Offset seismogram positions by the given mount \fIdx/dy\fP [Default is no offset].
If \fIdy\fP is not given it is set equal to \fIdx\fP\&.
.UNINDENT
.INDENT 0.0
.TP
\fB\-Ea\fP|\fBb\fP|\fBk\fP|\fBd\fP|\fBn\fP[\fIn\fP]|\fBu\fP[\fIn\fP]
Choose profile type (the type of Y axis).
.INDENT 7.0
.INDENT 3.5
\fBa\fP: azimuth profile.
.sp
\fBb\fP: back\-azimuth profile.
.sp
\fBk\fP: epicentral distance (in km) profile.
.sp
\fBd\fP: epicentral distance (in degree) profile.
.sp
\fBn\fP: trace number profile. The \fIY\fP position of first trace is numbered as \fIn\fP [Default \fIn\fP is 0].
.sp
\fBu\fP: user defined profile. The \fIY\fP positions are determined by SAC header variable user\fIn\fP, default using user0.
.UNINDENT
.UNINDENT
.UNINDENT
.INDENT 0.0
.TP
\fB\-F\fP[\fBi\fP][\fBq\fP][\fBr\fP]
Data preprocess before plotting.
.INDENT 7.0
.INDENT 3.5
\fBi\fP: integral
.sp
\fBq\fP: square
.sp
\fBr\fP: remove mean value
.UNINDENT
.UNINDENT
.sp
\fBi|q|r\fP can repeat multiple times. For example, \fB\-Frii\fP will convert acceleration to displacement.
The order of \fBi|q|r\fP controls the order of the data processing.
.UNINDENT
.INDENT 0.0
.TP
\fB\-G\fP[\fBp\fP|\fBn\fP][\fB+g\fP\fIfill\fP][\fB+z\fP\fIzero\fP][\fB+t\fP\fIt0/t1\fP]
Paint positive or negative portion of traces.
If only \fB\-G\fP is used, default to fill the positive portion black.
.INDENT 7.0
.INDENT 3.5
\fBp\fP|\fBn\fP controls the painting of positive \fBp\fPortion or \fBn\fPegative portion.
Repeat \fB\-G\fP option to specify fills for positive and negative portions, respectively.
.sp
\fB+g\fP\fIfill\fP: color to fill
.sp
\fB+t\fP\fIt0/t1\fP: paint traces between t0 and t1 only. The reference time of \fIt0\fP and \fIt1\fP is determined by \fB\-T\fP option.
.sp
\fB+z\fP\fIzero\fP: define zero line. From \fIzero\fP to top is positive portion, from \fIzero\fP to bottom is negative portion.
.UNINDENT
.UNINDENT
.UNINDENT
.INDENT 0.0
.TP
\fB\-K\fP (more ...)
Do not finalize the PostScript plot.
.UNINDENT
.INDENT 0.0
.TP
\fB\-M\fP\fIsize\fP[\fIu\fP][\fI/alpha\fP]
Vertical scaling.
.INDENT 7.0
.INDENT 3.5
\fIsize\fP[\fIu\fP]: scale all traces \fIsize\fP[\fIu\fP] on a map. The default unit is PROJ_LENGTH_UNIT.
The scaling factor is defined as \fByscale = size*(north\-south)/(depmax\-depmin)/map_height\fP\&.
.sp
\fIsize\fP/\fIalpha\fP:
.INDENT 0.0
.INDENT 3.5
\fIalpha\fP < 0, use the same scaling factor for all traces. The scaling factor will scale the first trace to \fIsize\fP[\fIu\fP].
.sp
\fIalpha\fP = 0, multiply all traces by \fIsize\fP\&. No unit is allowed.
.sp
\fIalpha\fP > 0, multiply all traces by \fBsize*r^alpha\fP, r is the distance range in km.
.UNINDENT
.UNINDENT
.UNINDENT
.UNINDENT
.UNINDENT
.INDENT 0.0
.TP
\fB\-O\fP (more ...)
Append to existing PostScript plot.
.UNINDENT
.INDENT 0.0
.TP
\fB\-P\fP (more ...)
Select "Portrait" plot orientation.
.UNINDENT
.INDENT 0.0
.TP
\fB\-Q\fP
Plot traces vertically.
.UNINDENT
.INDENT 0.0
.TP
\fB\-S\fP[\fBi\fP]\fIscale\fP[\fIunit\fP]
Sets time scale in seconds per <unit> while plotting on geographic plots.
Append c, i, or p to indicate cm, inch or points as the unit. Use PROJ_LENGTH_UNIT if \fIunit\fP is omitted.
Use \-S\fBi\fP\fIscale\fP\fIunit\fP to give the reciprocal scale, i.e. cm per second or inch per second.
.UNINDENT
.INDENT 0.0
.TP
\fB\-T\fP[\fB+t\fP\fIn\fP][\fB+r\fP\fIreduce_vel\fP][\fB+s\fP\fIshift\fP]
Time alignment and shift.
.INDENT 7.0
.INDENT 3.5
\fB+t\fP\fItmark\fP: align all trace along time mark. \fItmark\fP are \-5(b), \-4(e), \-3(o), \-2(a), 0\-9(t0\-t9).
.sp
\fB+r\fP\fIreduce_vel\fP: reduce velocity in km/s.
.sp
\fB+s\fP\fIshift\fP: shift all traces by \fIshift\fP seconds.
.UNINDENT
.UNINDENT
.UNINDENT
.INDENT 0.0
.TP
\fB\-U\fP[[\fIjust\fP]/\fIdx\fP/\fIdy\fP/][\fBc\fP|\fIlabel\fP] (more ...)
Draw GMT time stamp logo on plot.
.UNINDENT
.INDENT 0.0
.TP
\fB\-V\fP[\fIlevel\fP] (more ...)
Select verbosity level [c].  
.UNINDENT
.INDENT 0.0
.TP
\fB\-W\fP\fIpen\fP
Set pen attributes for all traces unless overruled by \fIpen\fP specified in \fIsaclist\fP\&.
[Defaults: width = default, color = black, style = solid].
.UNINDENT
.sp
\fB\-X\fP[\fBa\fP|\fBc\fP|\fBf\fP|\fBr\fP][\fIx\-shift\fP[\fBu\fP]]
.INDENT 0.0
.TP
\fB\-Y\fP[\fBa\fP|\fBc\fP|\fBf\fP|\fBr\fP][\fIy\-shift\fP[\fBu\fP]] (more ...)
Shift plot origin.
.UNINDENT
.INDENT 0.0
.TP
\fB\-h\fP[\fBi\fP|\fBo\fP][\fIn\fP][\fB+c\fP][\fB+d\fP][\fB+r\fP\fIremark\fP][\fB+r\fP\fItitle\fP] (more ...)
Skip or produce header record(s).  
.UNINDENT
.INDENT 0.0
.TP
\fB\-t\fP[\fItransp\fP] (more ...)
Set PDF transparency level in percent.
.UNINDENT
.INDENT 0.0
.TP
\fB\-^\fP or just \fB\-\fP
Print a short message about the syntax of the command, then exits (NOTE: on Windows just use \fB\-\fP).
.TP
\fB\-+\fP or just \fB+\fP
Print an extensive usage (help) message, including the explanation of
any module\-specific option (but not the GMT common options), then exits.
.TP
\fB\-?\fP or no arguments
Print a complete usage (help) message, including the explanation of all options, then exits.
.UNINDENT
.SH EXAMPLES
.sp
To plot a single seismogram seis.SAC (generated by SAC command \fBfuncgen seismogram\fP)
and paint positive portion black and negative portion red:
.INDENT 0.0
.INDENT 3.5
.sp
.nf
.ft C
gmt pssac seis.SAC \-JX10c/5c \-R9/20/\-2/2 \-Baf \-Fr \-Gp+gblack \-Gn+gred > single.ps
.ft P
.fi
.UNINDENT
.UNINDENT
.sp
To plot several seismograms (generated by SAC command \fBdatagen sub tele *.z\fP) on a distance profile:
.INDENT 0.0
.INDENT 3.5
.sp
.nf
.ft C
gmt pssac *.z \-R200/1600/12/45 \-JX15c/5c \-Bx200+l\(aqT(s)\(aq \-By5+lDegree \-BWSen \e
    \-Ed \-M1.5c \-W0.5p,red > distance_profile.ps
.ft P
.fi
.UNINDENT
.UNINDENT
.sp
To plot seismograms (generated by SAC command \fBdatagen sub tele *.z\fP) on a geographic map:
.INDENT 0.0
.INDENT 3.5
.sp
.nf
.ft C
gmt pssac *.z \-JM15c \-R\-120/\-40/35/65 \-Baf \-M1i \-S300c \-K > map.ps
saclst stlo stla f *.z | gmt psxy \-J \-R \-St0.4c \-Gblack \-i1,2 \-O >> map.ps
.ft P
.fi
.UNINDENT
.UNINDENT
.SH SEE ALSO
.sp
psmeca,
pspolar,
pscoupe,
psvelo,
gmt, psbasemap, psxy
.SH REFERENCES
.sp
Refer to \fI\%SAC User Manual\fP for more details
on SAC format and SAC header variables.
.SH AUTHORS
.sp
Dongdong Tian, School of Earth and Space Sciences,
University of Science and Technology of China, Hefei, Anhui, China
.SH COPYRIGHT
2019, P. Wessel, W. H. F. Smith, R. Scharroo, J. Luis, and F. Wobbe
.\" Generated by docutils manpage writer.
.