psvelo(1) GMT psvelo(1)
NAME
psvelo - Plot velocity vectors, crosses, and wedges on maps
SYNOPSIS
psvelo [ table ] -Jparameters
-Rregion [ -Aparameters ] [ -B[p|s]parameters ] [ -Ecolor ] [
-Fcolor ] [ -Gcolor ] [ -K ] [ -L ] [ -N ] [ -O ] [ -P ] [
-Ssymbol/scale/conf/font_size ] [ [ -U[stamp] ] [ -V[level] ] [
-Wpen ] [ -Xx_offset ] [ -Yy_offset ] [ -dinodata ] [ -eregexp ] [
-hheaders ] [ -iflags ] [ -ttransp ] [ -:[i|o] ]
Note: No space is allowed between the option flag and the associated
arguments.
DESCRIPTION
psvelo reads data values from files [or standard input] and generates
PostScript code that will plot velocity arrows on a map. Most options
are the same as for psxy, except -S. The PostScript code is written to
standard output. The previous version (psvelomeca) is now obsolete. It
has been replaced by psvelo and psmeca.
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 a|)
Select map projection.
-Rwest/east/south/north[/zmin/zmax][+r][+uunit]
west, east, south, and north specify the region of interest, and
you may specify them in decimal degrees or in
[A+-]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 for grid creation, give Rcodelon/lat/nx/ny, where
code is a 2-character combination of L, C, R (for left, center,
or right) and T, M, B for top, middle, or bottom. e.g., BL for
lower left. This indicates which point on a rectangular region
the lon/lat coordinate refers to, and the grid dimensions nx and
ny with grid spacings via -I is used to create the corresponding
region. Alternatively, specify the name of an existing grid
file and the -R settings (and grid spacing, if applicable) are
copied from the grid. Appending +uunit expects projected (Carte-
sian) 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.
-S
Selects the meaning of the columns in the data file and the figure to
be plotted.
-Sevelscale/confidence/fontsize.
Velocity ellipses in (N,E) convention. Vscale sets the scaling
of the velocity arrows. This scaling gives inches (unless c, i,
or p is appended). Confidence sets the 2-dimensional confidence
limit for the ellipse, e.g., 0.95 for 95% confidence ellipse.
Fontsize sets the size of the text in points. The ellipse will
be filled with the color or shade specified by the -G option
[default transparent]. The arrow and the circumference of the
ellipse will be drawn with the pen attributes specified by the
-W option. Parameters are expected to be in the following col-
umns:
1,2: longitude, latitude of station (-: option interchanges
order) 3,4: eastward, northward velocity (-: option inter-
changes order) 5,6: uncertainty of eastward, northward veloc-
ities (1-sigma) (-: option interchanges order) 7: correlation
between eastward and northward components 8: name of station
(optional).
-Snbarscale.
Anisotropy bars. Barscale sets the scaling of the bars This
scaling gives inches (unless c, i, or p is appended). Parame-
ters are expected to be in the following columns:
1,2: longitude, latitude of station (-: option interchanges
order) 3,4: eastward, northward components of anisotropy vec-
tor (-: option interchanges order)
-Srvelscale/confidence/fontsize
Velocity ellipses in rotated convention. Vscale sets the scaling
of the velocity arrows. This scaling gives inches (unless c, i,
or p is appended). Confidence sets the 2-dimensional confidence
limit for the ellipse, e.g., 0.95 for 95% confidence ellipse.
Fontsize sets the size of the text in points. The ellipse will
be filled with the color or shade specified by the -G option
[default transparent]. The arrow and the circumference of the
ellipse will be drawn with the pen attributes specified by the
-W option. Parameters are expected to be in the following col-
umns:
1,2: longitude, latitude, of station (-: option interchanges
order) 3,4: eastward, northward velocity (-: option inter-
changes order) 5,6: semi-major, semi-minor axes 7:
counter-clockwise angle, in degrees, from horizontal axis to
major axis of ellipse. 8: name of station (optional)
-Swwedge_scale/wedge_mag.
Rotational wedges. Wedge_scale sets the size of the wedges in
inches (unless c, i, or p is appended). Values are multiplied by
Wedge_mag before plotting. For example, setting Wedge_mag to
1.e7 works well for rotations of the order of 100 nanoradi-
ans/yr. Use -G to set the fill color or shade for the wedge, and
-E to set the color or shade for the uncertainty. Parameters are
expected to be in the following columns:
1,2: longitude, latitude, of station (-: option interchanges
order) 3: rotation in radians 4: rotation uncertainty in
radians
-Sxcross_scale
gives Strain crosses. Cross_scale sets the size of the cross in
inches (unless c, i, or p is appended). Parameters are expected
to be in the following columns:
1,2: longitude, latitude, of station (-: option interchanges
order) 3: eps1, the most extensional eigenvalue of strain
tensor, with extension taken positive. 4: eps2, the most
compressional eigenvalue of strain tensor, with extension
taken positive. 5: azimuth of eps2 in degrees CW from North.
OPTIONAL ARGUMENTS
-Aparameters
Modify vector parameters. For vector heads, append vector head
size [Default is 9p]. See Vector Attributes for specifying
additional attributes.
-B[p|s]parameters (more a|)
Set map boundary frame and axes attributes.
-DSigma_scale
can be used to rescale the uncertainties of velocities (-Se and
-Sr) and rotations (-Sw). Can be combined with the confidence
variable.
-Efill Sets the color or shade used for filling uncertainty wedges
(-Sw) or velocity error ellipses (-Se or -Sr). [If -E is not
specified, the uncertainty regions will be transparent.]
-Ffill Sets the color or shade used for frame and annotation. [Default
is black]
-Gfill Specify color (for symbols/polygons) or pattern (for polygons)
[Default is black]. Optionally, specify -Gpicon_size/pattern,
where pattern gives the number of the image pattern (1-90) OR
the name of a icon-format file. icon_size sets the unit size in
inches. To invert black and white pixels, use -GP instead of
-Gp. See the CookBook for information on individual patterns.
-K (more a|)
Do not finalize the PostScript plot.
-L Draw lines. Ellipses and fault planes will have their outlines
drawn using current pen (see -W).
-N Do NOT skip symbols that fall outside the frame boundary speci-
fied by -R. [Default plots symbols inside frame only].
-O (more a|)
Append to existing PostScript plot.
-P (more a|)
Select aPortraita plot orientation.
-U[[just]/dx/dy/][c|label] (more a|)
Draw GMT time stamp logo on plot.
-V[level] (more a|)
Select verbosity level [c].
-W Set pen attributes for velocity arrows, ellipse circumference
and fault plane edges. [Defaults: width = default, color =
black, style = solid].
-X[a|c|f|r][x-shift[u]]
-Y[a|c|f|r][y-shift[u]] (more a|)
Shift plot origin.
-dinodata (more a|)
Replace input columns that equal nodata with NaN.
-e[~]^<i>apattern^<i>a | -e[~]/regexp/[i] (more a|)
Only accept data records that match the given pattern.
-h[i|o][n][+c][+d][+rremark][+rtitle] (more a|)
Skip or produce header record(s).
-icols[+l][+sscale][+ooffset][,^<i>a|] (more a|)
Select input columns and transformations (0 is first column).
-t[transp] (more a|)
Set PDF transparency level in percent.
-:[i|o] (more a|)
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 just use -).
-+ or just +
Print an extensive usage (help) message, including the explana-
tion 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 all options, then exits.
VECTOR ATTRIBUTES
Several modifiers may be appended to the vector-producing options to
specify the placement of vector heads, their shapes, and the justifica-
tion 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, a for
arrow [Default], i for tail, A for plain arrow, and I for plain
tail. 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, a for
arrow [Default], i for tail, A for plain arrow, and I for plain
tail. 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].
+hshape sets the shape of the vector head (range -2/2). Default is
controlled by MAP_VECTOR_SHAPE [0].
+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 [for-
ward]. 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 cir-
cles. 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 are 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.
EXAMPLES
The following should make big red arrows with green ellipses, outlined
in red. Note that the 39% confidence scaling will give an ellipse which
fits inside a rectangle of dimension Esig by Nsig.
gmt psvelo << END -h2 -R-10/10/-10/10 -W0.25p,red -Ggreen -L -Se0.2/0.39/18 \
-B1g1 -Jx0.4/0.4 -A0.3p -P -V > test.ps
#Long. Lat. Evel Nvel Esig Nsig CorEN SITE
#(deg) (deg) (mm/yr) (mm/yr)
0. -8. 0.0 0.0 4.0 6.0 0.500 4x6
-8. 5. 3.0 3.0 0.0 0.0 0.500 33
0. 0. 4.0 6.0 4.0 6.0 0.500
-5. -5. 6.0 4.0 6.0 4.0 0.500 6x4
5. 0. -6.0 4.0 6.0 4.0 -0.500 -6x4
0. -5. 6.0 -4.0 6.0 4.0 -0.500 6x-4
END
This example should plot some residual rates of rotation in the Western
Transverse Ranges, California. The wedges will be dark gray, with light
gray wedges to represent the 2-sigma uncertainties.
gmt psvelo << END -Sw0.4/1.e7 -W0.75p -Gdarkgray -Elightgray -h1 -D2 -Jm2.2 \
-R240./243./32.5/34.75 -Bf10ma60m/WeSn -P > test.ps
#lon lat spin(rad/yr) spin_sigma (rad/yr)
241.4806 34.2073 5.65E-08 1.17E-08
241.6024 34.4468 -4.85E-08 1.85E-08
241.0952 34.4079 4.46E-09 3.07E-08
241.2542 34.2581 1.28E-07 1.59E-08
242.0593 34.0773 -6.62E-08 1.74E-08
241.0553 34.5369 -2.38E-07 4.27E-08
241.1993 33.1894 -2.99E-10 7.64E-09
241.1084 34.2565 2.17E-08 3.53E-08
END
SEE ALSO
psmeca(1), pspolar(1), pscoupe(1), gmt(1), psbasemap(1), psxy(1)
REFERENCES
Bomford, G., Geodesy, 4th ed., Oxford University Press, 1980.
AUTHORS
Kurt L. Feigl, Department of Geology and Geophysics at University of
Wisconsin-Madison, Madison, Wisconsin, USA
Genevieve Patau, Laboratory of Seismogenesis, Institut de Physique du
Globe de Paris, Departement de Sismologie, Paris, France
COPYRIGHT
2017, P. Wessel, W. H. F. Smith, R. Scharroo, J. Luis, and F. Wobbe
5.4.2 Jun 24, 2017 psvelo(1)
gmt5 5.4.2 - Generated Thu Jun 29 16:02:44 CDT 2017