1PSSAC(1) GMT PSSAC(1)
2
3
4
6 pssac - Plot seismograms in SAC format on maps
7
9 pssac [ saclist|SACfiles ] -Jparameters
10 -Rregion [ -B[p|s]parameters ] [ -C[t0/t1] ] [ -Ddx[/dy] ] [
11 -Ea|b|k|d|n[n]|u[n] ] [ -F[i][q][r] ] [
12 -G[p|n][+gfill][+zzero][+tt0/t1] ] [ -K ] [ -Msize[u][/alpha] ] [ -O
13 ] [ -P ] [ -Q ] [ -S[i]scale[unit] ] [ -T[+tn][+rre‐
14 duce_vel][+sshift] ] [ -U[stamp] ] [ -V[level] ] [ -Wpen ] [
15 -Xx_offset ] [ -Yy_offset ] [ -hheaders ] [ -ttransp ]
16
17 Note: No space is allowed between the option flag and the associated
18 arguments.
19
21 pssac reads SACfiles in SAC format or reads filenames and controlling
22 parameters from saclist [or standard input] and generates PostScript
23 that will plot seismograms on a map. The PostScript code is written to
24 standard output.
25
27 SACfiles
28 SAC files to plot on a map. Only evenly spaced SAC data is sup‐
29 ported.
30
31 saclist
32 One ASCII data table file holding a number of data columns. If
33 saclist is not given then we read from standard input. Parame‐
34 ters are expected to be in the following columns:
35 filename [X Y [pen]]
36
37 filename is the name of SAC file to plot. X and Y are the posi‐
38 tion of seismograms to plot on a map. On linear plots, the
39 default X is the begin time of SAC file, which will be adjusted
40 if -T option is used, the default Y is determined by -E option.
41 On geographic plots, the default X and Y are station longitude
42 and latitude specified in SAC header. The X and Y given here
43 will override the position determined by command line options.
44 pen, if given, will override the pen from -W option for current
45 SAC file only.
46
47 -Jparameters (more ...)
48 Select map projection.
49
50 -Rwest/east/south/north[/zmin/zmax][+r][+uunit]
51 west, east, south, and north specify the region of interest, and
52 you may specify them in decimal degrees or in
53 [±]dd:mm[:ss.xxx][W|E|S|N] format Append +r if lower left and
54 upper right map coordinates are given instead of w/e/s/n. The
55 two shorthands -Rg and -Rd stand for global domain (0/360 and
56 -180/+180 in longitude respectively, with -90/+90 in latitude).
57 Alternatively for grid creation, give Rcodelon/lat/nx/ny, where
58 code is a 2-character combination of L, C, R (for left, center,
59 or right) and T, M, B for top, middle, or bottom. e.g., BL for
60 lower left. This indicates which point on a rectangular region
61 the lon/lat coordinate refers to, and the grid dimensions nx and
62 ny with grid spacings via -I is used to create the corresponding
63 region. Alternatively, specify the name of an existing grid
64 file and the -R settings (and grid spacing, if applicable) are
65 copied from the grid. Appending +uunit expects projected (Carte‐
66 sian) coordinates compatible with chosen -J and we inversely
67 project to determine actual rectangular geographic region. For
68 perspective view (-p), optionally append /zmin/zmax. In case of
69 perspective view (-p), a z-range (zmin, zmax) can be appended to
70 indicate the third dimension. This needs to be done only when
71 using the -Jz option, not when using only the -p option. In the
72 latter case a perspective view of the plane is plotted, with no
73 third dimension.
74
76 -B[p|s]parameters (more ...)
77 Set map boundary frame and axes attributes.
78
79 -C[t0/t1]
80 Read and plot seismograms in timewindow between t0 and t1 only.
81 t0 and t1 are relative to a reference time specified by -T. If
82 -T option is not specified, use the reference time (kzdate and
83 kztime) defined in SAC header instead. If only -C is used,
84 t0/t1 is determined as xmin/xmax from -R option.
85
86 -Ddx[/dy]
87 Offset seismogram positions by the given mount dx/dy [Default is
88 no offset]. If dy is not given it is set equal to dx.
89
90 -Ea|b|k|d|n[n]|u[n]
91 Choose profile type (the type of Y axis).
92 a: azimuth profile.
93
94 b: back-azimuth profile.
95
96 k: epicentral distance (in km) profile.
97
98 d: epicentral distance (in degree) profile.
99
100 n: trace number profile. The Y position of first trace is
101 numbered as n [Default n is 0].
102
103 u: user defined profile. The Y positions are determined by
104 SAC header variable usern, default using user0.
105
106 -F[i][q][r]
107 Data preprocess before plotting.
108 i: integral
109
110 q: square
111
112 r: remove mean value
113
114 i|q|r can repeat multiple times. For example, -Frii will convert
115 acceleration to displacement. The order of i|q|r controls the
116 order of the data processing.
117
118 -G[p|n][+gfill][+zzero][+tt0/t1]
119 Paint positive or negative portion of traces. If only -G is
120 used, default to fill the positive portion black.
121 p|n controls the painting of positive portion or negative
122 portion. Repeat -G option to specify fills for positive and
123 negative portions, respectively.
124
125 +gfill: color to fill
126
127 +tt0/t1: paint traces between t0 and t1 only. The reference
128 time of t0 and t1 is determined by -T option.
129
130 +zzero: define zero line. From zero to top is positive por‐
131 tion, from zero to bottom is negative portion.
132
133 -K (more ...)
134 Do not finalize the PostScript plot.
135
136 -Msize[u][/alpha]
137 Vertical scaling.
138 size[u]: scale all traces size[u] on a map. The default unit
139 is PROJ_LENGTH_UNIT. The scaling factor is defined as yscale
140 = size*(north-south)/(depmax-depmin)/map_height.
141
142 size/alpha:
143 alpha < 0, use the same scaling factor for all traces.
144 The scaling factor will scale the first trace to size[u].
145
146 alpha = 0, multiply all traces by size. No unit is
147 allowed.
148
149 alpha > 0, multiply all traces by size*r^alpha, r is the
150 distance range in km.
151
152 -O (more ...)
153 Append to existing PostScript plot.
154
155 -P (more ...)
156 Select "Portrait" plot orientation.
157
158 -Q Plot traces vertically.
159
160 -S[i]scale[unit]
161 Sets time scale in seconds per <unit> while plotting on geo‐
162 graphic plots. Append c, i, or p to indicate cm, inch or points
163 as the unit. Use PROJ_LENGTH_UNIT if unit is omitted. Use -Sis‐
164 caleunit to give the reciprocal scale, i.e. cm per second or
165 inch per second.
166
167 -T[+tn][+rreduce_vel][+sshift]
168 Time alignment and shift.
169 +ttmark: align all trace along time mark. tmark are -5(b),
170 -4(e), -3(o), -2(a), 0-9(t0-t9).
171
172 +rreduce_vel: reduce velocity in km/s.
173
174 +sshift: shift all traces by shift seconds.
175
176 -U[[just]/dx/dy/][c|label] (more ...)
177 Draw GMT time stamp logo on plot.
178
179 -V[level] (more ...)
180 Select verbosity level [c].
181
182 -Wpen Set pen attributes for all traces unless overruled by pen speci‐
183 fied in saclist. [Defaults: width = default, color = black,
184 style = solid].
185
186 -X[a|c|f|r][x-shift[u]]
187
188 -Y[a|c|f|r][y-shift[u]] (more ...)
189 Shift plot origin.
190
191 -h[i|o][n][+c][+d][+rremark][+rtitle] (more ...)
192 Skip or produce header record(s).
193
194 -t[transp] (more ...)
195 Set PDF transparency level in percent.
196
197 -^ or just -
198 Print a short message about the syntax of the command, then
199 exits (NOTE: on Windows just use -).
200
201 -+ or just +
202 Print an extensive usage (help) message, including the explana‐
203 tion of any module-specific option (but not the GMT common
204 options), then exits.
205
206 -? or no arguments
207 Print a complete usage (help) message, including the explanation
208 of all options, then exits.
209
211 To plot a single seismogram seis.SAC (generated by SAC command funcgen
212 seismogram) and paint positive portion black and negative portion red:
213
214 gmt pssac seis.SAC -JX10c/5c -R9/20/-2/2 -Baf -Fr -Gp+gblack -Gn+gred > single.ps
215
216 To plot several seismograms (generated by SAC command datagen sub tele
217 *.z) on a distance profile:
218
219 gmt pssac *.z -R200/1600/12/45 -JX15c/5c -Bx200+l'T(s)' -By5+lDegree -BWSen \
220 -Ed -M1.5c -W0.5p,red > distance_profile.ps
221
222 To plot seismograms (generated by SAC command datagen sub tele *.z) on
223 a geographic map:
224
225 gmt pssac *.z -JM15c -R-120/-40/35/65 -Baf -M1i -S300c -K > map.ps
226 saclst stlo stla f *.z | gmt psxy -J -R -St0.4c -Gblack -i1,2 -O >> map.ps
227
229 psmeca, pspolar, pscoupe, psvelo, gmt, psbasemap, psxy
230
232 Refer to SAC User Manual for more details on SAC format and SAC header
233 variables.
234
236 Dongdong Tian, School of Earth and Space Sciences, University of Sci‐
237 ence and Technology of China, Hefei, Anhui, China
238
240 2019, P. Wessel, W. H. F. Smith, R. Scharroo, J. Luis, and F. Wobbe
241
242
243
244
2455.4.5 Feb 24, 2019 PSSAC(1)