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.
-Rwest/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.
-:[i|o] (more ...)
Swap 1st and 2nd column on input and/or output.
-^
Print a short message about the syntax of the command, then exits. (NOTE: on Windows use just -)
+
Print an extensive usage (help) message, including the explanation of any module-specific option (but not the GMT common options), then exits
-?
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

See Also

gmt, psbasemap, psxy

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