1PSPOLAR(1) GMT PSPOLAR(1)
2
3
4
6 pspolar - Plot polarities on the inferior focal half-sphere on maps
7
9 pspolar [ table ] -Dlon/lat -Jparameters
10 -Rregion
11 -Msize -S<symbol><size> [ -B[p|s]parameters ] [
12 -Clon/lat[/dash_width/pointsize] ] [ -Ecolor ] [ -Fcolor ] [ -Gcolor
13 ] [ -K ] [ -L ] [ -N ] [ -O ] [ -Qmode[args] ] [ -Tan‐
14 gle/form/justify/fontsize ] [ -U[stamp] ] [ -V[level] ] [ -Wpen ] [
15 -Xx_offset ] [ -Yy_offset ] [ -dinodata ] [ -eregexp ] [ -hheaders ] [
16 -iflags ] [ -ttransp ] [ -:[i|o] ]
17
18 Note: No space is allowed between the option flag and the associated
19 arguments.
20
22 pspolar reads data values from files [or standard input] and generates
23 PostScript code that will plot stations on focal mechanisms on a map.
24 The PostScript code is written to standard output.
25
26 Parameters are expected to be in the following columns:
27
28 1,2,3: station_code, azimuth, take-off angle
29
30 4: polarity:
31
32 · compression can be c,C,u,U,+
33
34 · rarefaction can be d,D,r,R,-
35
36 · not defined is anything else
37
39 table One or more ASCII (or binary, see -bi[ncols][type]) data table
40 file(s) holding a number of data columns. If no tables are given
41 then we read from standard input.
42
43 -Jparameters (more ...)
44 Select map projection.
45
46 -Rwest/east/south/north[/zmin/zmax][+r][+uunit]
47 west, east, south, and north specify the region of interest, and
48 you may specify them in decimal degrees or in
49 [±]dd:mm[:ss.xxx][W|E|S|N] format Append +r if lower left and
50 upper right map coordinates are given instead of w/e/s/n. The
51 two shorthands -Rg and -Rd stand for global domain (0/360 and
52 -180/+180 in longitude respectively, with -90/+90 in latitude).
53 Alternatively for grid creation, give Rcodelon/lat/nx/ny, where
54 code is a 2-character combination of L, C, R (for left, center,
55 or right) and T, M, B for top, middle, or bottom. e.g., BL for
56 lower left. This indicates which point on a rectangular region
57 the lon/lat coordinate refers to, and the grid dimensions nx and
58 ny with grid spacings via -I is used to create the corresponding
59 region. Alternatively, specify the name of an existing grid
60 file and the -R settings (and grid spacing, if applicable) are
61 copied from the grid. Appending +uunit expects projected (Carte‐
62 sian) coordinates compatible with chosen -J and we inversely
63 project to determine actual rectangular geographic region. For
64 perspective view (-p), optionally append /zmin/zmax. In case of
65 perspective view (-p), a z-range (zmin, zmax) can be appended to
66 indicate the third dimension. This needs to be done only when
67 using the -Jz option, not when using only the -p option. In the
68 latter case a perspective view of the plane is plotted, with no
69 third dimension.
70
71 -Dlongitude/latitude
72 Maps the bubble at given longitude and latitude point.
73
74 -Msize Sets the size of the beach ball to plot polarities in. Size is
75 in default units (unless c, i, or p is appended).
76
77 -S<symbol_type><size>
78 Selects symbol_type and symbol size. Size is in default inits
79 (unless c, i, or p is appended). Choose symbol type from st(a)r,
80 (c)ircle, (d)iamond, (h)exagon, (i)nverted triangle, (p)oint,
81 (s)quare, (t)riangle, (x)cross.
82
84 -B[p|s]parameters (more ...)
85 Set map boundary frame and axes attributes.
86
87 -C Offsets focal mechanisms to the latitude and longitude specified
88 in the last two columns of the input file.
89
90 -Ecolor
91 Selects filling of symbols for stations in extensive quadrants.
92 Set the color [Default is 250]. If -Ecolor is the same as
93 -Fcolor, use -e to outline.
94
95 -Fcolor
96 Sets background color of the beach ball. Default is no fill.
97
98 -Gcolor
99 Selects filling of symbols for stations in compressional quad‐
100 rants. Set the color [Default is black].
101
102 -K (more ...)
103 Do not finalize the PostScript plot.
104
105 -N Does not skip symbols that fall outside map border [Default
106 plots points inside border only].
107
108 -O (more ...)
109 Append to existing PostScript plot.
110
111 -P (more ...)
112 Select "Portrait" plot orientation.
113
114 -Qmode[args]
115 Sets one or more attributes; repeatable. The various combina‐
116 tions are
117
118 -Qe[pen]
119 Outline symbols in extensive quadrants using pen or the
120 default pen (see -W).
121
122 -Qf[pen]
123 Outline the beach ball using pen or the default pen (see
124 -W).
125
126 -Qg[pen]
127 Outline symbols in compressional quadrants using pen or
128 the default pen (see -W).
129
130 -Qh Use special format derived from HYPO71 output
131
132 -Qshalf-size[+vv_size[vecspecs]]
133 Plots S polarity azimuth. S polarity is in last column.
134 Append +v to select a vector and append head size and any
135 vector specifications. If +v is given without arguments
136 then we default to +v0.3i+e+gblack [Default is a line
137 segment]. Give half-size in default units (unless c, i,
138 or p is appended). See Vector Attributes for specifying
139 additional attributes.
140
141 -Qtpen Set pen color to write station code. Default uses the
142 default pen (see -W).
143
144 -Tangle/form/justify/fontsize in points
145 To write station code. [Default is 0.0/0/5/12].
146
147 -U[[just]/dx/dy/][c|label] (more ...)
148 Draw GMT time stamp logo on plot.
149
150 -V[level] (more ...)
151 Select verbosity level [c].
152
153 -W[-|+][pen][attr] (more ...)
154 Set current pen attributes [Defaults: width = default, color =
155 black, style = solid].
156
157 -X[a|c|f|r][x-shift[u]]
158
159 -Y[a|c|f|r][y-shift[u]] (more ...)
160 Shift plot origin.
161
162 -dinodata (more ...)
163 Replace input columns that equal nodata with NaN.
164
165 -e[~]"pattern" | -e[~]/regexp/[i] (more ...)
166 Only accept data records that match the given pattern.
167
168 -icols[+l][+sscale][+ooffset][,...] (more ...)
169 Select input columns and transformations (0 is first column).
170
171 -t[transp] (more ...)
172 Set PDF transparency level in percent.
173
174 -:[i|o] (more ...)
175 Swap 1st and 2nd column on input and/or output.
176
177 -^ or just -
178 Print a short message about the syntax of the command, then
179 exits (NOTE: on Windows just use -).
180
181 -+ or just +
182 Print an extensive usage (help) message, including the explana‐
183 tion of any module-specific option (but not the GMT common
184 options), then exits.
185
186 -? or no arguments
187 Print a complete usage (help) message, including the explanation
188 of all options, then exits.
189
191 Several modifiers may be appended to the vector-producing options to
192 specify the placement of vector heads, their shapes, and the justifica‐
193 tion of the vector. Below, left and right refers to the side of the
194 vector line when viewed from the start point to the end point of the
195 segment:
196 +aangle sets the angle of the vector head apex [30].
197
198 +b places a vector head at the beginning of the vector path [none].
199 Optionally, append t for a terminal line, c for a circle, a for
200 arrow [Default], i for tail, A for plain arrow, and I for plain
201 tail. Further append l|r to only draw the left or right side of
202 this head [both sides].
203
204 +e places a vector head at the end of the vector path [none].
205 Optionally, append t for a terminal line, c for a circle, a for
206 arrow [Default], i for tail, A for plain arrow, and I for plain
207 tail. Further append l|r to only draw the left or right side of
208 this head [both sides].
209
210 +g-|fill turns off vector head fill (if -) or sets the vector head
211 fill [Default fill is used, which may be no fill].
212
213 +hshape sets the shape of the vector head (range -2/2). Default is
214 controlled by MAP_VECTOR_SHAPE [0].
215
216 +l draws half-arrows, using only the left side of specified heads
217 [both sides].
218
219 +m places a vector head at the mid-point the vector path [none].
220 Append f or r for forward or reverse direction of the vector [for‐
221 ward]. Optionally, append t for a terminal line, c for a circle, or
222 a for arrow head [Default]. Further append l|r to only draw the
223 left or right side of this head [both sides]. Cannot be combined
224 with +b or +e.
225
226 +nnorm scales down vector attributes (pen thickness, head size) with
227 decreasing length, where vectors shorter than norm will have their
228 attributes scaled by length/norm [arrow attributes remains invariant
229 to length].
230
231 +oplon/plat specifies the oblique pole for the great or small cir‐
232 cles. Only needed for great circles if +q is given.
233
234 +p[-][pen] sets the vector pen attributes. If pen has a leading -
235 then the head outline is not drawn. [Default pen is used, and head
236 outline is drawn]
237
238 +q means the input angle, length data instead represent the start
239 and stop opening angles of the arc segment relative to the given
240 point.
241
242 +r draws half-arrows, using only the right side of specified heads
243 [both sides].
244
245 +t[b|e]trim will shift the beginning or end point (or both) along
246 the vector segment by the given trim; append suitable unit. If the
247 modifiers b|e are not used then trim may be two values separated by
248 a slash, which is used to specify different trims for the two ends.
249 Positive trims will shorted the vector while negative trims will
250 lengthen it [no trim].
251
252 In addition, all but circular vectors may take these modifiers:
253 +jjust determines how the input x,y point relates to the vector.
254 Choose from beginning [default], end, or center.
255
256 +s means the input angle, length are instead the x, y coordinates of
257 the vector end point.
258
259 Finally, Cartesian vectors may take these modifiers:
260 +zscale[unit] expects input dx,dy vector components and uses the
261 scale to convert to polar coordinates with length in given unit.
262
264 gmt pspolar -R239/240/34/35.2 -JM8c -N -Sc0.4 -h1 -D39.5/34.5 -M5 << END > test.ps
265 #stat azim ih pol
266 0481 11 147 c
267 6185 247 120 d
268 0485 288 114 +
269 0490 223 112 -
270 0487 212 109 .
271 END
272
273 or
274
275 gmt pspolar -R239/240/34/35.2 -JM8c -N -Sc0.4 -h1 -D239.5/34.5 -M5 <<END > test.ps
276 #Date Or. time stat azim ih
277 910223 1 22 0481 11 147 ipu0
278 910223 1 22 6185 247 120 ipd0
279 910223 1 22 0485 288 114 epu0
280 910223 1 22 0490 223 112 epd0
281 910223 1 22 0487 212 109 epu0
282 END
283
285 psmeca, psvelo, pscoupe, gmt, psbasemap, psxy
286
288 Bomford, G., Geodesy, 4th ed., Oxford University Press, 1980.
289
290 Aki, K. and P. Richards, Quantitative Seismology, Freeman, 1980.
291
293 Genevieve Patau, Laboratory of Seismogenesis
294 <http://www.ipgp.fr/rech/sismogenese/>, Institut de Physique du Globe
295 de Paris, Departement de Sismologie, Paris, France
296
298 2019, P. Wessel, W. H. F. Smith, R. Scharroo, J. Luis, and F. Wobbe
299
300
301
302
3035.4.5 Feb 24, 2019 PSPOLAR(1)