pspolar - Plot polarities on the inferior focal half-sphere on maps
pspolar [ table ] -Dlon/lat -Jparameters -Rregion -Msize -S<symbol><size> [ -B[p|s]parameters ] [ -Clon/lat[/dash_width/pointsize] ] [ -Ecolor ] [ -Fcolor ] [ -Gcolor ] [ -K ] [ -L ] [ -N ] [ -O ] [ -Qmode[args] ] [ -Tangle/form/justify/fontsize ] [ -U[just/dx/dy/][c|label] ] [ -V[level] ] [ -Wpen ] [ -Xx_offset ] [ -Yy_offset ] [ -ccopies ] [ -di<nodata> ] [ -h<headers> ] [ -i<flags> ] [ -:[i|o] ]
Note: No space is allowed between the option flag and the associated arguments.
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
Sets one or more attributes; repeatable. The various combinations are
-X[a|c|f|r][x-shift[u]]
Several modifiers may be appended to the vector-producing options to specify the placement of vector heads, their shapes, and the justification of the vector. Below, left and right refers to the side of the vector line when viewed from the start point to the end point of the segment:
+aangle sets the angle of the vector head apex [30].
+b places a vector head at the beginning of the vector path [none]. Optionally, append t for a terminal line, c for a circle, or a for arrow [Default]. Further append l|r to only draw the left or right side of this head [both sides].
+e places a vector head at the end of the vector path [none]. Optionally, append t for a terminal line, c for a circle, or a for arrow [Default]. Further append l|r to only draw the left or right side of this head [both sides].
+g-|fill turns off vector head fill (if -) or sets the vector head fill [Default fill is used, which may be no fill].
+l draws half-arrows, using only the left side of specified heads [both sides].
+m places a vector head at the mid-point the vector path [none]. Append f or r for forward or reverse direction of the vector [forward]. Optionally, append t for a terminal line, c for a circle, or a for arrow head [Default]. Further append l|r to only draw the left or right side of this head [both sides]. Cannot be combined with +b or +e.
+nnorm scales down vector attributes (pen thickness, head size) with decreasing length, where vectors shorter than norm will have their attributes scaled by length/norm [arrow attributes remains invariant to length].
+oplon/plat specifies the oblique pole for the great or small circles. Only needed for great circles if +q is given.
+p[-][pen] sets the vector pen attributes. If pen has a leading - then the head outline is not drawn. [Default pen is used, and head outline is drawn]
+q means the input angle, length data instead represent the start and stop opening angles of the arc segment relative to the given point.
+r draws half-arrows, using only the right side of specified heads [both sides].
t[b|e]trim will shift the beginning or end point (or both) along the vector segment by the given trim; append suitable unit. If the modifiers b|e are not used then trim may be two values separated by a slash, which is used to specify different trims for the two ends. Positive trims will shorted the vector while negative trims will lengthen it [no trim].
In addition, all but circular vectors may take these modifiers:
+jjust determines how the input x,y point relates to the vector. Choose from beginning [default], end, or center.
+s means the input angle, length is instead the x, y coordinates of the vector end point.
Finally, Cartesian vectors may take these modifiers:
+zscale[unit] expects input dx,dy vector components and uses the scale to convert to polar coordinates with length in given unit.
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
Bomford, G., Geodesy, 4th ed., Oxford University Press, 1980.
Aki, K. and P. Richards, Quantitative Seismology, Freeman, 1980.
Genevieve Patau, Laboratory of Seismogenesis <http://www.ipgp.fr/rech/sismogenese/>, Institut de Physique du Globe de Paris, Departement de Sismologie, Paris, France