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