NAME
pspolar - Plot polarities on the inferior focal half-sphere on maps
SYNOPSIS
pspolar files -Jparameters -Rwest/east/south/north -Dlongitude/latitude
-Msize -S<symbol><size> [ -Btickinfo ] [ -Clongitude/lati-
tude[/dash_width/pointsize] ] [ -Fred/green/blue ] [ -Gfill ] [ -gfill
] [ -H ] [ -h ] [ -K ] [ -L ] [ -N ] [ -O ] [ -P ] [ -sHalf-
size[v[[<v_width/h_length/h_width/shape>]][g[<r/g/b>]][l] [ -Tan-
gle/form/justify/fontsize ] [ -t<pen> ] [ -U[label] ] [ -V ] [ -Wpen ]
[ -Xx-shift ] [ -Yy-shift ] [ -ccopies ]
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,4
station_code, azimuth, take-off angle, polarity
polarity:
- compression can be c,C,u,U,+
- rarefaction can be d,D,r,R,-
- not defined is anything else
ARGUMENTS
files List one or more file-names. If no files are given, pspolar will
read standard input.
-J Selects the map projection. Scale is UNIT/degree, 1:xxxxx, or
width in UNIT (upper case modifier). UNIT is cm, inch, or m,
depending on the MEASURE_UNIT setting in .gmtdefaults4, but this
can be overridden on the command line by appending c, i, or m to
the scale/width value. For map height, max dimension, or min
dimension, append h, +, or - to the width, respectively.
CYLINDRICAL PROJECTIONS:
-Jclon0/lat0/scale (Cassini)
-Jjlon0/scale (Miller)
-Jmscale (Mercator - Greenwich and Equator as origin)
-Jmlon0/lat0/scale (Mercator - Give meridian and standard paral-
lel)
-Joalon0/lat0/azimuth/scale (Oblique Mercator - point and
azimuth)
-Joblon0/lat0/lon1/lat1/scale (Oblique Mercator - two points)
-Joclon0/lat0/lonp/latp/scale (Oblique Mercator - point and
pole)
-Jqlon0/scale (Equidistant Cylindrical Projection (Plate Car-
ree))
-Jtlon0/scale (TM - Transverse Mercator, with Equator as y = 0)
-Jtlon0/lat0/scale (TM - Transverse Mercator, set origin)
CONIC PROJECTIONS:
-Jblon0/lat0/lat1/lat2/scale (Albers)
-Jdlon0/lat0/lat1/lat2/scale (Equidistant)
-Jllon0/lat0/lat1/lat2/scale (Lambert)
MISCELLANEOUS PROJECTIONS:
-Jhlon0/scale (Hammer)
-Jilon0/scale (Sinusoidal)
-Jk[f|s]lon0/scale (Eckert IV (f) and VI (s))
-Jnlon0/scale (Robinson)
-Jrlon0/scale (Winkel Tripel)
-Jvlon0/scale (Van der Grinten)
-Jwlon0/scale (Mollweide)
NON-GEOGRAPHICAL PROJECTIONS:
-Jp[a]scale[/origin] (polar (theta,r) coordinates, optional a
for azimuths and offset theta [0])
-Jxx-scale[l|ppow][/y-scale[l|ppow]][d] (Linear, log, and power
scaling)
More details can be found in the psbasemap man pages.
-R xmin, xmax, ymin, and ymax specify the Region of interest. For
geographic regions, these limits correspond to west, east,
south, and north 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 wesn. The
two shorthands -Rg -Rd stand for global domain (0/360 or
-180/+180 in longitude respectively, with -90/+90 in latitude).
For calendar time coordinates you may either give relative time
(relative to the selected TIME_EPOCH and in the selected
TIME_UNIT; append t to -JX|x), or absolute time of the form
[date]T[clock] (append T to -JX|x). At least one of date and
clock must be present; the T is always required. The date string
must be of the form [-]yyyy[-mm[-dd]] (Gregorian calendar) or
yyyy[-Www[-d]] (ISO week calendar), while the clock string must
be of the form hh:mm:ss[.xxx]. The use of delimiters and their
type and positions must be as indicated (however, input/output
and plotting formats are flexible).
-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.
If -Efill is the same as -Ffill, use -e to outline.
-e[pen]
Outline symbols in extensive quadrants using pen or the default
pen (see -W).
-Ffill Sets background color of the beach ball. Default is no fill.
-f[pen]
Outline the beach ball using pen or the default pen (see -W).
-Gfill Selects filling of symbols for stations in compressional quad-
rants. Set the shade (0-255) or color (r/g/b) [Default is
0/0/0].
-g[pen]
Outline symbols in compressional quadrants using pen or the
default pen (see -W).
-H Input file(s) has Header record(s). Number of header records can
be changed by editing your .gmtdefaults4 file. If used, GMT
default is 1 header record. Use -Hi if only input data should
have header records [Default will write out header records if
the input data have them].
-h Use special format derived from HYPO71 output
-K More PostScript code will be appended later [Default terminates
the plot system].
-N Does NOT skip symbols that fall outside map border [Default
plots points inside border only].
-O Selects Overlay plot mode [Default initializes a new plot sys-
tem].
-P Selects Portrait plotting mode [GMT Default is Landscape, see
gmtdefaults to change this].
-shalf-size/[V[<v_width/h_length/h_width/shape>]][G<r/g/b>][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.
-Tangle/form/justify/fontsize in points
To write station code. [Default is 0.0/0/5/12].
-t<pen>
Set pen color to write station code. Default uses the default
pen (see -W).
-U Draw Unix System time stamp on plot. User may specify where the
lower left corner of the stamp should fall on the page relative
to lower left corner of plot. Optionally, append a label, or c
lute coordinates; the default (r) will reset plot origin. Give
c to center plot using current page size.
-: Toggles between (longitude,latitude) and (latitude,longitude)
input and/or output. [Default is (longitude,latitude)]. Append
i to select input only or o to select output only. [Default
affects both].
-c Specifies the number of plot copies. [Default is 1].
EXAMPLES
pspolar -R239/240/34/35.2 -JM8 -N -Sc0.4 -H1 -D239.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
pspolar -R239/240/34/35.2 -JM8 -N -Sc0.4 -H1 -D239.5/34.5 -M5 -h
<<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(GMTMANSECTION), psbasemap(GMTMANSECTION), psxy(GMTMANSECTION)
REFERENCES
Bomford, G., Geodesy, 4th ed., Oxford University Press, 1980.
Aki, K. and P. Richards, Quantitative Seismology, Freeman, 1980.
AUTHORS
Genevieve Patau
Seismology Dept.
Institut de Physique du Globe de Paris
(patau@ipgp.jussieu.fr)
GMT4.0 1 Oct 2004 PSPOLAR(GMTMANSECTION)
Man(1) output converted with
man2html