1PSVELO(1) GMT PSVELO(1)
2
6 psvelo - Plot velocity vectors, crosses, and wedges on maps
7
9 psvelo [ table ] -Jparameters
10 -Rregion [ -Aparameters ] [ -B[p|s]parameters ] [ -Ecolor ] [
11 -Fcolor ] [ -Gcolor ] [ -K ] [ -L ] [ -N ] [ -O ] [ -P ] [
12 -Ssymbol/scale/conf/font_size ] [ [ -U[stamp] ] [ -V[level] ] [
13 -Wpen ] [ -Xx_offset ] [ -Yy_offset ] [ -dinodata ] [ -eregexp ] [
14 -hheaders ] [ -iflags ] [ -ttransp ] [ -:[i|o] ]
15 Note: No space is allowed between the option flag and the associated
17 arguments.
18
20 psvelo reads data values from files [or standard input] and generates
21 PostScript code that will plot velocity arrows on a map. Most options
22 are the same as for psxy, except -S. The PostScript code is written to
23 standard output. The previous version (psvelomeca) is now obsolete. It
24 has been replaced by psvelo and psmeca.
25
27 table One or more ASCII (or binary, see -bi[ncols][type]) data table
28 file(s) holding a number of data columns. If no tables are given
29 then we read from standard input.
30 -Jparameters (more ...)
32 Select map projection.
33 -Rwest/east/south/north[/zmin/zmax][+r][+uunit]
35 west, east, south, and north specify the region of interest, and
36 you may specify them in decimal degrees or in
37 [±]dd:mm[:ss.xxx][W|E|S|N] format Append +r if lower left and
38 upper right map coordinates are given instead of w/e/s/n. The
39 two shorthands -Rg and -Rd stand for global domain (0/360 and
40 -180/+180 in longitude respectively, with -90/+90 in latitude).
41 Alternatively for grid creation, give Rcodelon/lat/nx/ny, where
42 code is a 2-character combination of L, C, R (for left, center,
43 or right) and T, M, B for top, middle, or bottom. e.g., BL for
44 lower left. This indicates which point on a rectangular region
45 the lon/lat coordinate refers to, and the grid dimensions nx and
46 ny with grid spacings via -I is used to create the corresponding
47 region. Alternatively, specify the name of an existing grid
48 file and the -R settings (and grid spacing, if applicable) are
49 copied from the grid. Appending +uunit expects projected (Carte‐
50 sian) coordinates compatible with chosen -J and we inversely
51 project to determine actual rectangular geographic region. For
52 perspective view (-p), optionally append /zmin/zmax. In case of
53 perspective view (-p), a z-range (zmin, zmax) can be appended to
54 indicate the third dimension. This needs to be done only when
55 using the -Jz option, not when using only the -p option. In the
56 latter case a perspective view of the plane is plotted, with no
57 third dimension.
58 -S
60 Selects the meaning of the columns in the data file and the figure to
62 be plotted.
63 -Sevelscale/confidence/fontsize.
64 Velocity ellipses in (N,E) convention. Vscale sets the scaling
65 of the velocity arrows. This scaling gives inches (unless c, i,
66 or p is appended). Confidence sets the 2-dimensional confidence
67 limit for the ellipse, e.g., 0.95 for 95% confidence ellipse.
68 Fontsize sets the size of the text in points. The ellipse will
69 be filled with the color or shade specified by the -G option
70 [default transparent]. The arrow and the circumference of the
71 ellipse will be drawn with the pen attributes specified by the
72 -W option. Parameters are expected to be in the following col‐
73 umns:
74 1,2: longitude, latitude of station (-: option interchanges
75 order) 3,4: eastward, northward velocity (-: option inter‐
76 changes order) 5,6: uncertainty of eastward, northward veloc‐
77 ities (1-sigma) (-: option interchanges order) 7: correlation
78 between eastward and northward components 8: name of station
79 (optional).
80 -Snbarscale.
82 Anisotropy bars. Barscale sets the scaling of the bars This
83 scaling gives inches (unless c, i, or p is appended). Parame‐
84 ters are expected to be in the following columns:
85 1,2: longitude, latitude of station (-: option interchanges
86 order) 3,4: eastward, northward components of anisotropy vec‐
87 tor (-: option interchanges order)
88 -Srvelscale/confidence/fontsize
90 Velocity ellipses in rotated convention. Vscale sets the scaling
91 of the velocity arrows. This scaling gives inches (unless c, i,
92 or p is appended). Confidence sets the 2-dimensional confidence
93 limit for the ellipse, e.g., 0.95 for 95% confidence ellipse.
94 Fontsize sets the size of the text in points. The ellipse will
95 be filled with the color or shade specified by the -G option
96 [default transparent]. The arrow and the circumference of the
97 ellipse will be drawn with the pen attributes specified by the
98 -W option. Parameters are expected to be in the following col‐
99 umns:
100 1,2: longitude, latitude, of station (-: option interchanges
101 order) 3,4: eastward, northward velocity (-: option inter‐
102 changes order) 5,6: semi-major, semi-minor axes 7:
103 counter-clockwise angle, in degrees, from horizontal axis to
104 major axis of ellipse. 8: name of station (optional)
105 -Swwedge_scale/wedge_mag.
107 Rotational wedges. Wedge_scale sets the size of the wedges in
108 inches (unless c, i, or p is appended). Values are multiplied by
109 Wedge_mag before plotting. For example, setting Wedge_mag to
110 1.e7 works well for rotations of the order of 100 nanoradi‐
111 ans/yr. Use -G to set the fill color or shade for the wedge, and
112 -E to set the color or shade for the uncertainty. Parameters are
113 expected to be in the following columns:
114 1,2: longitude, latitude, of station (-: option interchanges
115 order) 3: rotation in radians 4: rotation uncertainty in
116 radians
117 -Sxcross_scale
119 gives Strain crosses. Cross_scale sets the size of the cross in
120 inches (unless c, i, or p is appended). Parameters are expected
121 to be in the following columns:
122 1,2: longitude, latitude, of station (-: option interchanges
123 order) 3: eps1, the most extensional eigenvalue of strain
124 tensor, with extension taken positive. 4: eps2, the most
125 compressional eigenvalue of strain tensor, with extension
126 taken positive. 5: azimuth of eps2 in degrees CW from North.
127
129 -Aparameters
130 Modify vector parameters. For vector heads, append vector head
131 size [Default is 9p]. See Vector Attributes for specifying
132 additional attributes.
133 -B[p|s]parameters (more ...)
135 Set map boundary frame and axes attributes.
136 -DSigma_scale
138 can be used to rescale the uncertainties of velocities (-Se and
139 -Sr) and rotations (-Sw). Can be combined with the confidence
140 variable.
141 -Efill Sets the color or shade used for filling uncertainty wedges
143 (-Sw) or velocity error ellipses (-Se or -Sr). [If -E is not
144 specified, the uncertainty regions will be transparent.]
145 -Ffill Sets the color or shade used for frame and annotation. [Default
147 is black]
148 -Gfill Specify color (for symbols/polygons) or pattern (for polygons)
150 [Default is black]. Optionally, specify -Gpicon_size/pattern,
151 where pattern gives the number of the image pattern (1-90) OR
152 the name of a icon-format file. icon_size sets the unit size in
153 inches. To invert black and white pixels, use -GP instead of
154 -Gp. See the CookBook for information on individual patterns.
155 -K (more ...)
157 Do not finalize the PostScript plot.
158 -L Draw lines. Ellipses and fault planes will have their outlines
160 drawn using current pen (see -W).
161 -N Do NOT skip symbols that fall outside the frame boundary speci‐
163 fied by -R. [Default plots symbols inside frame only].
164 -O (more ...)
166 Append to existing PostScript plot.
167 -P (more ...)
169 Select "Portrait" plot orientation.
170 -U[[just]/dx/dy/][c|label] (more ...)
172 Draw GMT time stamp logo on plot.
173 -V[level] (more ...)
175 Select verbosity level [c].
176 -W Set pen attributes for velocity arrows, ellipse circumference
178 and fault plane edges. [Defaults: width = default, color =
179 black, style = solid].
180 -X[a|c|f|r][x-shift[u]]
182 -Y[a|c|f|r][y-shift[u]] (more ...)
184 Shift plot origin.
185 -dinodata (more ...)
187 Replace input columns that equal nodata with NaN.
188 -e[~]"pattern" | -e[~]/regexp/[i] (more ...)
190 Only accept data records that match the given pattern.
191 -h[i|o][n][+c][+d][+rremark][+rtitle] (more ...)
193 Skip or produce header record(s).
194 -icols[+l][+sscale][+ooffset][,...] (more ...)
196 Select input columns and transformations (0 is first column).
197 -t[transp] (more ...)
199 Set PDF transparency level in percent.
200 -:[i|o] (more ...)
202 Swap 1st and 2nd column on input and/or output.
203 -^ or just -
205 Print a short message about the syntax of the command, then
206 exits (NOTE: on Windows just use -).
207 -+ or just +
209 Print an extensive usage (help) message, including the explana‐
210 tion of any module-specific option (but not the GMT common
211 options), then exits.
212 -? or no arguments
214 Print a complete usage (help) message, including the explanation
215 of all options, then exits.
216
218 Several modifiers may be appended to the vector-producing options to
219 specify the placement of vector heads, their shapes, and the justifica‐
220 tion of the vector. Below, left and right refers to the side of the
221 vector line when viewed from the start point to the end point of the
222 segment:
223 +aangle sets the angle of the vector head apex [30].
224 +b places a vector head at the beginning of the vector path [none].
226 Optionally, append t for a terminal line, c for a circle, a for
227 arrow [Default], i for tail, A for plain arrow, and I for plain
228 tail. Further append l|r to only draw the left or right side of
229 this head [both sides].
230 +e places a vector head at the end of the vector path [none].
232 Optionally, append t for a terminal line, c for a circle, a for
233 arrow [Default], i for tail, A for plain arrow, and I for plain
234 tail. Further append l|r to only draw the left or right side of
235 this head [both sides].
236 +g-|fill turns off vector head fill (if -) or sets the vector head
238 fill [Default fill is used, which may be no fill].
239 +hshape sets the shape of the vector head (range -2/2). Default is
241 controlled by MAP_VECTOR_SHAPE [0].
242 +l draws half-arrows, using only the left side of specified heads
244 [both sides].
245 +m places a vector head at the mid-point the vector path [none].
247 Append f or r for forward or reverse direction of the vector [for‐
248 ward]. Optionally, append t for a terminal line, c for a circle, or
249 a for arrow head [Default]. Further append l|r to only draw the
250 left or right side of this head [both sides]. Cannot be combined
251 with +b or +e.
252 +nnorm scales down vector attributes (pen thickness, head size) with
254 decreasing length, where vectors shorter than norm will have their
255 attributes scaled by length/norm [arrow attributes remains invariant
256 to length].
257 +oplon/plat specifies the oblique pole for the great or small cir‐
259 cles. Only needed for great circles if +q is given.
260 +p[-][pen] sets the vector pen attributes. If pen has a leading -
262 then the head outline is not drawn. [Default pen is used, and head
263 outline is drawn]
264 +q means the input angle, length data instead represent the start
266 and stop opening angles of the arc segment relative to the given
267 point.
268 +r draws half-arrows, using only the right side of specified heads
270 [both sides].
271 +t[b|e]trim will shift the beginning or end point (or both) along
273 the vector segment by the given trim; append suitable unit. If the
274 modifiers b|e are not used then trim may be two values separated by
275 a slash, which is used to specify different trims for the two ends.
276 Positive trims will shorted the vector while negative trims will
277 lengthen it [no trim].
278 In addition, all but circular vectors may take these modifiers:
280 +jjust determines how the input x,y point relates to the vector.
281 Choose from beginning [default], end, or center.
282 +s means the input angle, length are instead the x, y coordinates of
284 the vector end point.
285 Finally, Cartesian vectors may take these modifiers:
287 +zscale[unit] expects input dx,dy vector components and uses the
288 scale to convert to polar coordinates with length in given unit.
289
291 The following should make big red arrows with green ellipses, outlined
292 in red. Note that the 39% confidence scaling will give an ellipse which
293 fits inside a rectangle of dimension Esig by Nsig.
294 gmt psvelo << END -h2 -R-10/10/-10/10 -W0.25p,red -Ggreen -L -Se0.2/0.39/18 \
296 -B1g1 -Jx0.4/0.4 -A0.3p -P -V > test.ps
297 #Long. Lat. Evel Nvel Esig Nsig CorEN SITE
298 #(deg) (deg) (mm/yr) (mm/yr)
299 0. -8. 0.0 0.0 4.0 6.0 0.500 4x6
300 -8. 5. 3.0 3.0 0.0 0.0 0.500 3x3
301 0. 0. 4.0 6.0 4.0 6.0 0.500
302 -5. -5. 6.0 4.0 6.0 4.0 0.500 6x4
303 5. 0. -6.0 4.0 6.0 4.0 -0.500 -6x4
304 0. -5. 6.0 -4.0 6.0 4.0 -0.500 6x-4
305 END
306 This example should plot some residual rates of rotation in the Western
308 Transverse Ranges, California. The wedges will be dark gray, with light
309 gray wedges to represent the 2-sigma uncertainties.
310 gmt psvelo << END -Sw0.4/1.e7 -W0.75p -Gdarkgray -Elightgray -h1 -D2 -Jm2.2 \
312 -R240./243./32.5/34.75 -Bf10ma60m/WeSn -P > test.ps
313 #lon lat spin(rad/yr) spin_sigma (rad/yr)
314 241.4806 34.2073 5.65E-08 1.17E-08
315 241.6024 34.4468 -4.85E-08 1.85E-08
316 241.0952 34.4079 4.46E-09 3.07E-08
317 241.2542 34.2581 1.28E-07 1.59E-08
318 242.0593 34.0773 -6.62E-08 1.74E-08
319 241.0553 34.5369 -2.38E-07 4.27E-08
320 241.1993 33.1894 -2.99E-10 7.64E-09
321 241.1084 34.2565 2.17E-08 3.53E-08
322 END
323
325 psmeca, pspolar, pscoupe, gmt, psbasemap, psxy
326
328 Bomford, G., Geodesy, 4th ed., Oxford University Press, 1980.
329
331 Kurt L. Feigl, Department of Geology and Geophysics at University of
332 Wisconsin-Madison, Madison, Wisconsin, USA
333 Genevieve Patau, Laboratory of Seismogenesis, Institut de Physique du
335 Globe de Paris, Departement de Sismologie, Paris, France
336
338 2019, P. Wessel, W. H. F. Smith, R. Scharroo, J. Luis, and F. Wobbe
3395.4.5 Feb 24, 2019 PSVELO(1)