1PSPOLAR(1)                            GMT                           PSPOLAR(1)
2
3
4

NAME

6       pspolar - Plot polarities on the inferior focal half-sphere on maps
7

SYNOPSIS

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

DESCRIPTION

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

REQUIRED ARGUMENTS

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

OPTIONAL ARGUMENTS

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

VECTOR ATTRIBUTES

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

EXAMPLES

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

SEE ALSO

285       psmeca, psvelo, pscoupe, gmt, psbasemap, psxy
286

REFERENCES

288       Bomford, G., Geodesy, 4th ed., Oxford University Press, 1980.
289
290       Aki, K. and P. Richards, Quantitative Seismology, Freeman, 1980.
291

AUTHORS

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