1PSCOUPE(1) GMT PSCOUPE(1)
2
6 pscoupe - Plot cross-sections of focal mechanisms
7
9 pscoupe [ files ] -Jparameters
10 -Rregion -Aparameters [ -B[p|s]parameters ] [ -Ecolor ] [
11 -Fmode[args] ] [ -Gcolor ] [ -K ] [ -L[pen] ] [ -M ] [ -N ] [ -O
12 ] [ -Q ] [ -S<symbol><scale>[/d] ] [ -Tn ] [ -U[stamp] ] [
13 -V[level] ] [ -Wpen ] [ -Xx_offset ] [ -Yy_offset ] [ -Zcpt ] [
14 -dinodata ] [ -eregexp ] [ -hheaders ] [ -iflags ] [ -ttransp ] [
15 -:[i|o] ]
16 Note: No space is allowed between the option flag and the associated
18 arguments.
19
21 pscoupe reads data values from files [or standard input] and generates
22 PostScript code that will plot symbols, lines or polygons on a
23 cross-section. Focal mechanisms may be specified and require additional
24 columns of data. The PostScript code is written to standard output.
25 Unless -Q is used, new file is created with the new coordinates (x, y)
27 and the mechanism (from lower focal half-sphere for horizontal plane,
28 to half-sphere behind a vertical plane). When the plane is not horizon‐
29 tal, - north direction becomes upwards steepest descent direction of
30 the plane (u) - east direction becomes strike direction of the plane
31 (s) - down direction (= north^east) becomes u^s Axis angles are defined
32 in the same way as in horizontal plane in the new system. Moment ten‐
33 sor (initially in r, t, f system that is up, south, east) is defined in
34 (-u^s, -u, s) system.
35
37 table One or more ASCII (or binary, see -bi[ncols][type]) data table
38 file(s) holding a number of data columns. If no tables are given
39 then we read from standard input.
40 -Jparameters (more ...)
42 Select map projection.
43 -Rwest/east/south/north[/zmin/zmax][+r][+uunit]
45 west, east, south, and north specify the region of interest, and
46 you may specify them in decimal degrees or in
47 [±]dd:mm[:ss.xxx][W|E|S|N] format Append +r if lower left and
48 upper right map coordinates are given instead of w/e/s/n. The
49 two shorthands -Rg and -Rd stand for global domain (0/360 and
50 -180/+180 in longitude respectively, with -90/+90 in latitude).
51 Alternatively for grid creation, give Rcodelon/lat/nx/ny, where
52 code is a 2-character combination of L, C, R (for left, center,
53 or right) and T, M, B for top, middle, or bottom. e.g., BL for
54 lower left. This indicates which point on a rectangular region
55 the lon/lat coordinate refers to, and the grid dimensions nx and
56 ny with grid spacings via -I is used to create the corresponding
57 region. Alternatively, specify the name of an existing grid
58 file and the -R settings (and grid spacing, if applicable) are
59 copied from the grid. Appending +uunit expects projected (Carte‐
60 sian) coordinates compatible with chosen -J and we inversely
61 project to determine actual rectangular geographic region. For
62 perspective view (-p), optionally append /zmin/zmax. In case of
63 perspective view (-p), a z-range (zmin, zmax) can be appended to
64 indicate the third dimension. This needs to be done only when
65 using the -Jz option, not when using only the -p option. In the
66 latter case a perspective view of the plane is plotted, with no
67 third dimension. If frame is defined from cross-section parame‐
68 ters (see -A this option is not taken into account, but must be
69 present.
70 -A selects the cross-section.
72 -Aalon1/lat1/lon2/lat2/dip/p_width/dmin/dmax[f]
74 lon and lat are the longitude and latitude of points 1
75 and 2 limiting the length of the cross-section. dip is
76 the dip of the plane on which the cross-section is made.
77 p_width is the width of the cross-section on each side of
78 a vertical plane or above and under an oblique plane.
79 dmin and dmax are the distances min and max from horizon‐
80 tal plane, along steepest descent direction. Add f to get
81 the frame from the cross-section parameters.
82 -Ablon1/lat1/strike/p_length/dip/p_width/dmin/dmax[f]
84 lon1 and lat1 are the longitude and latitude of the
85 beginning of the cross-section. strike is the azimuth of
86 the direction of the cross-section. p_length is the
87 length along which the cross-section is made. The other
88 parameters are the same as for -Aa option.
89 -Acx1/y1/x2/y2/dip/p_width/dmin/dmax[f]
91 The same as -Aa option with x and y cartesian coordi‐
92 nates.
93 -Adx1/y1/strike/p_length/dip/p_width/dmin/dmax[f]
95 The same as -Ab option with x and y cartesian coordi‐
96 nates.
97 -S selects the meaning of the columns in the data file and the fig‐
99 ure to be plotted.
100 -Sascale[/fontsize[/offset[u]]]
102 Focal mechanisms in Aki and Richards convention. scale
103 adjusts the scaling of the radius of the "beach ball",
104 which will be proportional to the magnitude. The scale is
105 the size for magnitude = 5 in PROJ_LENGTH_UNIT (unless c,
106 i, or p is appended to indicate that the size information
107 is in units of cm, inches, meters, or points, respec‐
108 tively). Use the -T option to render the beach ball
109 transparent by drawing only the nodal planes and the cir‐
110 cumference. The color or shade of the compressive quad‐
111 rants can be specified with the -G option. The color or
112 shade of the extensive quadrants can be specified with
113 the -E option. Parameters are expected to be in the fol‐
114 lowing columns:
115 1,2: longitude, latitude of event (-: option inter‐
117 changes order)
118 3: depth of event in kilometers
120 4,5,6: strike, dip and rake
122 7: magnitude
124 8,9: not used; can be 0 0; allows use of the psmeca
126 file format
127 10: text string to appear above the beach ball
129 (default) or under (add u).
130 -Scscale
132 Focal mechanisms in Harvard CMT convention. scale adjusts
133 the scaling of the radius of the "beach ball", which will
134 be proportional to the magnitude. The scale is the size
135 for magnitude = 5 (that is M0 = 4E+23 dynes-cm.) in
136 PROJ_LENGTH_UNIT (unless c, i, or p is appended to indi‐
137 cate that the size information is in units of cm, inches,
138 meters, or points, respectively). Use the -T option to
139 render the beach ball transparent by drawing only the
140 nodal planes and the circumference. The color or shade of
141 the compressive quadrants can be specified with the -G
142 option. The color or shade of the extensive quadrants can
143 be specified with the -E option. Parameters are expected
144 to be in the following columns:
145 1,2: longitude, latitude of event (-: option inter‐
147 changes order)
148 3: depth of event in kilometers
150 4,5,6: strike, dip, and slip of plane 1
152 7,8,9: strike, dip, and slip of plane 2
154 10,11: mantissa and exponent of moment in dyne-cm (if
156 magnitude is uses instead of scalar moment,
157 magnitude is in column 10 and 0 must be in col‐
158 umn 11)
159 12,13: not used; can be 0 0; allows use of the psmeca
161 file format
162 14: text string to appear above the beach ball
164 (default) or under (add u).
165 -Spscale[/fontsize[/offset[u]]]
167 Focal mechanisms given with partial data on both planes.
168 scale adjusts the scaling of the radius of the "beach
169 ball", which will be proportional to the magnitude. The
170 scale is the size for magnitude = 5 in PROJ_LENGTH_UNIT
171 (unless c, i, or p is appended to indicate that the size
172 information is in units of cm, inches, meters, or points,
173 respectively). The color or shade of the compressive
174 quadrants can be specified with the -G option. The color
175 or shade of the extensive quadrants can be specified with
176 the -E option. Parameters are expected to be in the fol‐
177 lowing columns:
178 1,2: longitude, latitude of event (-: option inter‐
180 changes order)
181 3: depth
183 4,5: strike, dip of plane 1
185 6: strike of plane 2
187 7: must be -1/+1 for a normal/inverse fault
189 8: magnitude
191 9,10: not used; can be 0 0; allows use of the psmeca
193 file format
194 11: text string to appear above the beach ball
196 (default) or under (add u).
197 -Sm|d|zscale[/fontsize[/offset[u]]]
199 Seismic moment tensor (Harvard CMT, with zero trace).
200 scale adjusts the scaling of the radius of the "beach
201 ball", which will be proportional to the magnitude. The
202 scale is the size for magnitude = 5 (that is seismic
203 scalar moment = 4E+23 dynes-cm) in PROJ_LENGTH_UNIT
204 (unless c, i, or p is appended to indicate that the size
205 information is in units of cm, inches, meters, or points,
206 respectively). (-T0 option overlays best double couple
207 transparently.)
208 -Sdscale[/fontsize[/offset[u]]]
210 to plot the only double couple part of moment tensor.
211 -Szscale[/fontsize[/offset[u]]]
213 to plot anisotropic part of moment tensor (zero trace).
214 The color or shade of the compressive quadrants can be
215 specified with the -G option. The color or shade of the
216 extensive quadrants can be specified with the -E option.
217 Parameters are expected to be in the following columns:
218 1,2: longitude, latitude of event (-: option inter‐
220 changes order)
221 3: depth of event in kilometers
223 4,5,6,7,8,9:
225 mrr, mtt, mff, mrt, mrf, mtf in 10*exponent
226 dynes-cm
227 10: exponent
229 11,12: Not used; can be 0 0; allows use of the psmeca
231 file format
232 13: Text string to appear above the beach ball
234 (default) or under (add u).
235 -Sxscale[/fontsize[/offset[u]]]
237 Principal axis. scale adjusts the scaling of the radius
238 of the "beach ball", which will be proportional to the
239 magnitude. The scale is the size for magnitude = 5 (that
240 is seismic scalar moment = 4*10e+23 dynes-cm) in
241 PROJ_LENGTH_UNIT (unless c, i, or p is appended to indi‐
242 cate that the size information is in units of cm, inches,
243 meters, or points, respectively). (-T0 option overlays
244 best double couple transparently.)
245 -Syscale[/fontsize[/offset[u]]]
247 to plot the only double couple part of moment tensor.
248 -Stscale[/fontsize[/offset[u]]]
250 to plot anisotropic part of moment tensor (zero trace).
251 The color or shade of the compressive quadrants can be
252 specified with the -G option. The color or shade of the
253 extensive quadrants can be specified with the -E option.
254 Parameters are expected to be in the following columns:
255 1,2: longitude, latitude of event (-: option inter‐
257 changes order)
258 3: depth of event in kilometers
260 4,5,6,7,8,9,10,11,12:
262 value (in 10*exponent dynes-cm), azimuth,
263 plunge of the T, N, and P axes.
264 13: exponent
266 14,15: longitude, latitude at which to place beach
268 ball. Entries in these columns are necessary
269 with the -C option. Using 0,0 in columns 9 and
270 10 will plot the beach ball at the longitude,
271 latitude given in columns 1 and 2. The -:
272 option will interchange the order of columns
273 (1,2) and (9,10).
274 16: Text string to appear above the beach ball
276 (optional).
277
279 -B[p|s]parameters (more ...)
280 Set map boundary frame and axes attributes.
281 -Ecolor
283 Sets color or fill pattern for extensive quadrants [Default is
284 white].
285 -Fmode[args]
287 Sets one or more attributes; repeatable. The various combina‐
288 tions are
289 -Fssymbol[size[/fontsize[/offset[u]]]
291 selects a symbol instead of mechanism. Choose from the follow‐
292 ing: (c) circle, (d) diamond, (i) itriangle, (s) square, (t)
293 triangle, (x) cross. size is the symbol size in PROJ_LENGTH_UNIT
294 (unless c, i, or p is appended to indicate that the size infor‐
295 mation is in units of cm, inches, meters, or points, respec‐
296 tively). If size must be read, it must be in column 4 and the
297 text string will start in column 5. Parameters are expected to
298 be in the following columns:
299 1,2: longitude, latitude of event (-: option interchanges
301 order)
302 3: depth of event in kilometers
304 4: Text string to appear above the beach ball (default)
306 or under (add u).
307 -Fa[size][/P_symbol[/T_symbol]]
309 Computes and plots P and T axes with symbols. Optionally specify
310 size and (separate) P and T axis symbols from the following: (c)
311 circle, (d) diamond, (h) hexagon, (i) inverse triangle, (p)
312 point, (s) square, (t) triangle, (x) cross. [Default: 6p/cc]
313 -Fecolor
315 Sets the color or fill pattern for the T axis symbol. [Default
316 as set by -E]
317 -Fgcolor
319 Sets the color or fill pattern for the P axis symbol. [Default
320 as set by -G]
321 -Fp[pen]
323 Draws the P axis outline using current pen (see -W), or sets
324 pen attributes.
325 -Fr[color]
327 Draw a box behind the label (if any). [Default fill is white]
328 -Ft[pen]
330 Draws the T axis outline using current pen (see -W), or sets pen
331 attributes.
332 -Gcolor
334 Sets color or fill pattern for compressional quadrants [Default
335 is black].
336 -K (more ...)
338 Do not finalize the PostScript plot.
339 -L[pen]
341 Draws the "beach ball" outline using current pen (see -W) or
342 sets pen attributes.
343 -M Same size for any magnitude.
345 -N Does not skip symbols that fall outside map border [Default
347 plots points inside border only].
348 -O (more ...)
350 Append to existing PostScript plot.
351 -P (more ...)
353 Select "Portrait" plot orientation.
354 -Q Suppress the production of files with cross-section and mecha‐
356 nism information.
357 -T[num_of_planes]
359 Plots the nodal planes and outlines the bubble which is trans‐
360 parent. If num_of_planes is 0: both nodal planes are plotted;
361 1: only the first nodal plane is plotted; 2: only the second
362 nodal plane is plotted [Default: 0].
363 -U[[just]/dx/dy/][c|label] (more ...)
365 Draw GMT time stamp logo on plot.
366 -V[level] (more ...)
368 Select verbosity level [c].
369 -W[-|+][pen][attr] (more ...)
371 set pen attributes for text string or default pen attributes for
372 fault plane edges. [Defaults: width = default, color = black,
373 style = solid].
374 -X[a|c|f|r][x-shift[u]]
376 -Y[a|c|f|r][y-shift[u]] (more ...)
378 Shift plot origin.
379 -Zcpt Give a CPT and let compressive part color be determined by the
381 z-value in the third column.
382 -dinodata (more ...)
384 Replace input columns that equal nodata with NaN.
385 -e[~]"pattern" | -e[~]/regexp/[i] (more ...)
387 Only accept data records that match the given pattern.
388 -h[i|o][n][+c][+d][+rremark][+rtitle] (more ...)
390 Skip or produce header record(s).
391 -icols[+l][+sscale][+ooffset][,...] (more ...)
393 Select input columns and transformations (0 is first column).
394 -t[transp] (more ...)
396 Set PDF transparency level in percent.
397 -:[i|o] (more ...)
399 Swap 1st and 2nd column on input and/or output.
400 -^ or just -
402 Print a short message about the syntax of the command, then
403 exits (NOTE: on Windows just use -).
404 -+ or just +
406 Print an extensive usage (help) message, including the explana‐
407 tion of any module-specific option (but not the GMT common
408 options), then exits.
409 -? or no arguments
411 Print a complete usage (help) message, including the explanation
412 of all options, then exits.
413
415 psmeca, psvelo, pspolar, gmt, psbasemap, psxy
416
418 Bomford, G., Geodesy, 4th ed., Oxford University Press, 1980.
419 Aki, K. and P. Richards, Quantitative Seismology, Freeman, 1980.
421 F. A. Dahlen and Jeroen Tromp, Theoretical Seismology, Princeton, 1998,
423 p.167. Definition of scalar moment.
424 Cliff Frohlich, Cliff's Nodes Concerning Plotting Nodal Lines for P, Sh
426 and Sv
427 Seismological Research Letters, Volume 67, Number 1, January-February,
429 1996
430 Thorne Lay, Terry C. Wallace, Modern Global Seismology, Academic Press,
432 1995, p.384.
433 W.H. Press, S.A. Teukolsky, W.T. Vetterling, B.P. Flannery, Numerical
435 Recipes in C, Cambridge University press (routine jacobi)
436
438 Genevieve Patau, Laboratory of Seismogenesis
439 <http://www.ipgp.fr/rech/sismogenese/>, Institut de Physique du Globe
440 de Paris, Departement de Sismologie, Paris, France
441
443 2019, P. Wessel, W. H. F. Smith, R. Scharroo, J. Luis, and F. Wobbe
4445.4.5 Feb 24, 2019 PSCOUPE(1)