pspolar
pspolar - Plot polarities on the inferior focal half-sphere on maps
Synopsis
pspolar [ table ] -Jparameters
-R[unit]xmin/xmax/ymin/ymax[r]
-Msize -S<symbol><size>
[ -B[p|s]parameters ]
[ -Clon/lat[/dash_width/pointsize] ] [ -Fcolor ]
[ -Gfill ] [ -K ] [ -L ] [ -N ]
[ -O ] [ -Qmode[args] ]
[ -Tangle/form/justify/fontsize ]
[ -U[just/dx/dy/][c|label] ]
[ -V[level] ]
[ -Wpen ]
[ -X[a|c|f|r][x-shift[u]] ]
[ -Y[a|c|f|r][y-shift[u]] ]
[ -ccopies ]
[ -h[i|o][n][+c][+d][+rremark][+rtitle] ]
[ -icols[l][sscale][ooffset][,...] ]
[ -:[i|o] ]
Note: No space is allowed between the option flag and the associated arguments.
Description
pspolar reads data values from files [or standard input] and
generates PostScript code that will plot stations on focal mechanisms
on a map. The PostScript code is written to standard output.
Parameters are expected to be in the following columns:
- 1,2,3:
- station_code, azimuth, take-off angle
- 4:
polarity:
- compression can be c,C,u,U,+
- rarefaction can be d,D,r,R,-
- not defined is anything else
Required Arguments
- table
- One or more ASCII (or binary, see -bi[ncols][type]) data
table file(s) holding a number of data columns. If no tables are given
then we read from standard input.
- -Jparameters (more ...)
- Select map projection.
- -R[unit]west/east/south/north[/zmin/zmax][r]
- west, east, south, and north specify the region of interest,
and you may specify them in decimal degrees or in
[+-]dd:mm[:ss.xxx][W|E|S|N] format. Append r if lower left
and upper right map coordinates are given instead of w/e/s/n. The
two shorthands -Rg and -Rd stand for global domain (0/360
and -180/+180 in longitude respectively, with -90/+90 in latitude).
Alternatively, specify the name of an existing grid file and the
-R settings (and grid spacing, if applicable) are copied from
the grid. Using -Runit expects projected (Cartesian)
coordinates compatible with chosen -J and we inversely project
to determine actual rectangular geographic region.
For perspective view (-p), optionally append /zmin/zmax.
In case of perspective view (-p), a z-range (zmin, zmax)
can be appended to indicate the third dimension. This needs to be
done only when using the -Jz option, not when using only the
-p option. In the latter case a perspective view of the plane is
plotted, with no third dimension.
- -Dlongitude/latitude)
- Maps the bubble at given longitude and latitude point.
- -Msize)
- Sets the size of the beach ball to plot polarities in. Size is in
inch (unless c, i, m, or p is appended).
- -S<symbol_type><size>)
- Selects symbol_type and symbol size. Size is in inch (unless
c, i, m, or p is appended). Choose symbol type from
st(a)r, (c)ircle, (d)iamond, (h)exagon, (i)nverted
triangle, (p)oint, (s)quare, (t)riangle, (x)cross.
Optional Arguments
- -B[p|s]parameters (more ...)
- Set map boundary intervals.
- -C
- Offsets focal mechanisms to the latitude and longitude specified in
the last two columns of the input file.
- -Efill)
- Selects filling of symbols for stations in extensive quadrants. Set
the color [Default is 250]. If -Efill is the same as
-Ffill, use -e to outline.
- -Ffill
- Sets background color of the beach ball. Default is no fill.
- -Gfill
- Selects filling of symbols for stations in compressional quadrants.
Set the color [Default is black].
- -K (more ...)
- Do not finalize the PostScript plot.
- -N
- Does not skip symbols that fall outside map border [Default
plots points inside border only].
- -O (more ...)
- Append to existing PostScript plot.
- -P (more ...)
- Select “Portrait” plot orientation.
- -Qmode[args]
- Sets one or more attributes; repeatable. The various combinations are
- -Qe[pen])
- Outline symbols in extensive quadrants using pen or the default
pen (see -W).
- -Qf[pen])
- Outline the beach ball using pen or the default pen (see -W).
- -Qg[pen])
- Outline symbols in compressional quadrants using pen or the
default pen (see -W).
- -Qh
- Use special format derived from HYPO71 output
- -Qshalf-size/[V[v_width/h_length/h_width/shape]][Gcolor][L]
- Plots S polarity azimuth. S polarity is in last column. It may be a
vector (V option) or a segment. Give
half-size,v_width,h_length,h_width in inch (unless c, i,
m, or p is appended). [L] option is for outline.
- -Qtpen
- Set pen color to write station code. Default uses the default pen
(see -W).
- -Tangle/form/justify/fontsize in points
- To write station code. [Default is 0.0/0/5/12].
- -U[just/dx/dy/][c|label] (more ...)
- Draw GMT time stamp logo on plot.
- -V[level] (more ...)
- Select verbosity level [c].
- -Wpen
- Set current pen attributes [Defaults: width = default, color = black, style = solid].
-X[a|c|f|r][x-shift[u]]
- -Y[a|c|f|r][y-shift[u]] (more ...)
- Shift plot origin.
- -ccopies (more ...)
- Specify number of plot copies [Default is 1].
- -icols[l][sscale][ooffset][,...] (more ...)
- Select input columns (0 is first column).
- -:[i|o] (more ...)
- Swap 1st and 2nd column on input and/or output.
- -^ or just -
- Print a short message about the syntax of the command, then exits (NOTE: on Windows use just -).
- -+ or just +
- Print an extensive usage (help) message, including the explanation of
any module-specific option (but not the GMT common options), then exits.
- -? or no arguments
- Print a complete usage (help) message, including the explanation of
options, then exits.
- --version
- Print GMT version and exit.
- --show-datadir
- Print full path to GMT share directory and exit.
Examples
gmt pspolar -R239/240/34/35.2 -JM8c -N -Sc0.4 -h1 -D39.5/34.5 -M5 << END > test.ps
#stat azim ih pol
0481 11 147 c
6185 247 120 d
0485 288 114 +
0490 223 112 -
0487 212 109 .
END
or
gmt pspolar -R239/240/34/35.2 -JM8c -N -Sc0.4 -h1 -D239.5/34.5 -M5 <<END > test.ps
#Date Or. time stat azim ih
910223 1 22 0481 11 147 ipu0
910223 1 22 6185 247 120 ipd0
910223 1 22 0485 288 114 epu0
910223 1 22 0490 223 112 epd0
910223 1 22 0487 212 109 epu0
END
References
Bomford, G., Geodesy, 4th ed., Oxford University Press, 1980.
Aki, K. and P. Richards, Quantitative Seismology, Freeman, 1980.
Authors
Genevieve Patau, Laboratory of Seismogenesis <http://www.ipgp.fr/rech/sismogenese/>,
Institut de Physique du Globe de Paris, Departement de Sismologie, Paris, France