1PSSAC(1)                              GMT                             PSSAC(1)
2
3
4

NAME

6       pssac - Plot seismograms in SAC format on maps
7

SYNOPSIS

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

DESCRIPTION

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

REQUIRED ARGUMENTS

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

OPTIONAL ARGUMENTS

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

EXAMPLES

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

SEE ALSO

229       psmeca, pspolar, pscoupe, psvelo, gmt, psbasemap, psxy
230

REFERENCES

232       Refer to SAC User Manual for more details on SAC format and SAC  header
233       variables.
234

AUTHORS

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)
Impressum