NAME
psvelo - Plot velocity vectors, crosses, and wedges on maps
SYNOPSIS
psvelo files -Jparameters -Rwest/east/south/north [ -Btick-
info ] [ -Dsigma_scale ] [ -Fred/green/blue ] [
-Ered/green/blue ] [ -Gfill ] [ -H ] [ -K ] [ -L ] [ -N ] [
-O ] [ -P ] [ -Ssymbol/scale/conf/font_size ] [
-U[/dx/dy/][label] ] [ -V ] [ -Wpen ] [ -Xx-shift ] [ -Yy-
shift ] [ -: ] [ -ccopies ]
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.
ARGUMENTS
files List one or more file-names. If no files are given,
psvelo will read standard input.
-J Selects the map projection. Scale is UNIT/degree,
1:xxxxx. or width in UNIT (upper case modifier).
See psbasemap man pages for details.
-R west, east, south, and north specify the Region of
interest. To specify your boundaries in degrees and
minutes [and seconds], use the dd:mm[:ss] format.
-SSelects 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, m, 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 fol-
lowing columns:
1,2 longitude, latitude of station (-: option interchanges
order)
3,4 eastward, northward velocity (-: option interchanges
order)
5,6 uncertainty of eastward, northward velocities (1-sigma)
(-: option interchanges order)
7 correlation between eastward and northward components
8 name of station (optional). The bug which previously
wrote a garbage string if this argument was omitted has
been fixed.
-Snbarscale. Anisotropy bars. Barscale sets the scaling of
the bars This scaling gives inches (unless c, i, m, or p is
appended). Parameters are expected to be in the following
columns:
1,2 longitude, latitude of station (-: option interchanges
order)
3,4 eastward, northward components of anisotropy vector (-:
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, m, 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 fol-
lowing columns:
1,2 longitude, latitude, of station (-: option interchanges
order)
3,4 eastward, northward velocity (-: option interchanges
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, m, 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 nanoradians/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 Strain crosses. Cross_scale sets the size
of the cross in inches (unless c, i, m, 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 ten-
sor, with extension taken positive.
5 azimuth of eps2 in degrees CW from North.
-R west, east, south, and north specify the region of
interest.
To specify your boundaries in degrees and minutes, use the
dd:mm format by appending 'm'.
OPTIONS
No space between the option flag and the associated argu-
ments.
-A Arrow_width/Head_length/Head_width Size of arrow in
inches. [Default is 0.03/0.12/0.09].
-B Sets map boundary tickmark intervals. See psbasemap for
details.
-C Offset focal mechanisms to the latitude and longitude
specified in the last two columns of the input file.
-D Sigma_scale can be used to rescale the uncertainties of
velocities (-Se and -Sr) and rotations (-Sw). Can be com-
bined with the confidence variable.
-Ffill Sets the color or shade used for frame and annota-
tion. [Default is 0/0/0 (black)]
-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 tran-
sparent.]
-Gfill Select filling of ellipses, wedges, and focal mechan-
isms. By convention, the compressional quadrants of the
focal mechanism beach balls are shaded. Set the shade
(0-255) or color (r/g/b) [Default is 0/0/0]. Optionally,
specify -Gpicon_size/pattern, where pattern gives the number
of the image pattern (1-32) 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 pspat-
terns for information on individual patterns.
-H Input file(s) has Header record(s).
-K Means more PostScript code will be appended later.
[Default terminates the plot system]
-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 boudary
specified by -R.
[Default plots symbols inside frame only].
-O Selects Overlay plot mode. [Default initializes a new
plot system]
-P Selects Portrait plotting mode. [Default is Landscape]
-U Draw Unix System time stamp on plot. Append c to plot
the command line or a text string of your choice.
-V Verbose. Reports the name of the current file being
plotted, as well as number of records, and any scaling
information.
-W Set pen attributes for velocity arrows, ellipse cir-
cumference and fault plane edges.
[Defaults: width = 1, color = 0/0/0, texture = solid].
-X -Y Shift origin of plot by (x-shift,y-shift) inches
[Default is (1,1) for new plots, (0,0) for overlays].
-: Toggles between (longitude,latitude) and
(latitude,longitude) input/output. [Default is
(longitude,latitude)]
-c Specifies the number of plot copies. [Default is 1]
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.
psvelo << END -H2 -R-10/10/-10/10 -W1/255/0/0 \
-G0/255/0 -L -Se0.2/0.39/18 -B1g1 -Jx0.4/0.4 \
-A0.1/0.3/0.3 -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 3x3
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
The 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.
psvelo <<END -Sw0.4/1.e7 -W3 -G60 -E210 -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
gmtsystem, psbasemap, psxy.
REFERENCES
Bomford, G., Geodesy, 4th ed., Oxford University Press,
1980.
AUTHORS
Kurt Feigl
CNRS UMR 5562
Toulouse, France
(Kurt.Feigl@.cnes.fr)
Genevieve Patau
CNRS UMR 7580
Seismology Dept.
Institut de Physique du Globe de Paris
(patau@ipgp.jussieu.fr)