1PSVELO(1) Generic Mapping Tools PSVELO(1)
2
6 psvelo - Plot velocity vectors, crosses, and wedges on maps
7
9 psvelo files -Jparameters -Rwest/east/south/north[r] [ -B[p|s]parame‐
10 ters ] [ -Dsigma_scale ] [ -Fcolor ] [ -Ecolor ] [ -Gfill ] [
11 -H[i][nrec] ] [ -K ] [ -L ] [ -N ] [ -O ] [ -P ] [ -Ssym‐
12 bol/scale/conf/font_size ] [ -U[just/dx/dy/][c|label] ] [ -V ] [ -Wpen
13 ] [ -X[a|c|r][x-shift[u]] ] [ -Y[a|c|r][y-shift[u]] ] [ -:[i|o] ] [
14 -ccopies ]
15
17 psvelo reads data values from files [or standard input] and generates
18 PostScript code that will plot velocity arrows on a map. Most options
19 are the same as for psxy, except -S. The PostScript code is written to
20 standard output. The previous version (psvelomeca) is now obsolete. It
21 has been replaced by psvelo and psmeca.
22
25 files List one or more file-names. If no files are given, psvelo will
26 read standard input.
27 -J Selects the map projection. Scale is UNIT/degree, 1:xxxxx, or
29 width in UNIT (upper case modifier). UNIT is cm, inch, or m,
30 depending on the MEASURE_UNIT setting in .gmtdefaults4, but this
31 can be overridden on the command line by appending c, i, or m to
32 the scale/width value. When central meridian is optional,
33 default is center of longitude range on -R option. Default
34 standard parallel is the equator. For map height, max dimen‐
35 sion, or min dimension, append h, +, or - to the width, respec‐
36 tively.
37 More details can be found in the psbasemap man pages.
38 CYLINDRICAL PROJECTIONS:
40 -Jclon0/lat0/scale (Cassini)
42 -Jcyl_stere/[lon0/[lat0/]]scale (Cylindrical Stereographic)
43 -Jj[lon0/]scale (Miller)
44 -Jm[lon0/[lat0/]]scale (Mercator)
45 -Jmlon0/lat0/scale (Mercator - Give meridian and standard paral‐
46 lel)
47 -Jo[a]lon0/lat0/azimuth/scale (Oblique Mercator - point and
48 azimuth)
49 -Jo[b]lon0/lat0/lon1/lat1/scale (Oblique Mercator - two points)
50 -Joclon0/lat0/lonp/latp/scale (Oblique Mercator - point and
51 pole)
52 -Jq[lon0/[lat0/]]scale (Cylindrical Equidistant)
53 -Jtlon0/[lat0/]scale (TM - Transverse Mercator)
54 -Juzone/scale (UTM - Universal Transverse Mercator)
55 -Jy[lon0/[lat0/]]scale (Cylindrical Equal-Area)
56 CONIC PROJECTIONS:
58 -Jblon0/lat0/lat1/lat2/scale (Albers)
60 -Jdlon0/lat0/lat1/lat2/scale (Conic Equidistant)
61 -Jllon0/lat0/lat1/lat2/scale (Lambert Conic Conformal)
62 -Jpoly/[lon0/[lat0/]]scale ((American) Polyconic)
63 AZIMUTHAL PROJECTIONS:
65 -Jalon0/lat0[/horizon]/scale (Lambert Azimuthal Equal-Area)
67 -Jelon0/lat0[/horizon]/scale (Azimuthal Equidistant)
68 -Jflon0/lat0[/horizon]/scale (Gnomonic)
69 -Jglon0/lat0[/horizon]/scale (Orthographic)
70 -Jglon0/lat0/altitude/azimuth/tilt/twist/Width/Height/scale
71 (General Perspective).
72 -Jslon0/lat0[/horizon]/scale (General Stereographic)
73 MISCELLANEOUS PROJECTIONS:
75 -Jh[lon0/]scale (Hammer)
77 -Ji[lon0/]scale (Sinusoidal)
78 -Jkf[lon0/]scale (Eckert IV)
79 -Jk[s][lon0/]scale (Eckert VI)
80 -Jn[lon0/]scale (Robinson)
81 -Jr[lon0/]scale (Winkel Tripel)
82 -Jv[lon0/]scale (Van der Grinten)
83 -Jw[lon0/]scale (Mollweide)
84 NON-GEOGRAPHICAL PROJECTIONS:
86 -Jp[a]scale[/origin][r|z] (Polar coordinates (theta,r))
88 -Jxx-scale[d|l|ppow|t|T][/y-scale[d|l|ppow|t|T]] (Linear, log,
89 and power scaling)
90 -R west, east, south, and north specify the Region of interest, and
92 you may specify them in decimal degrees or in
93 [+-]dd:mm[:ss.xxx][W|E|S|N] format. Append r if lower left and
94 upper right map coordinates are given instead of w/e/s/n. The
95 two shorthands -Rg and -Rd stand for global domain (0/360 and
96 -180/+180 in longitude respectively, with -90/+90 in latitude).
97 Alternatively, specify the name of an existing grid file and the
98 -R settings (and grid spacing, if applicable) are copied from
99 the grid.
100 -SSelects the meaning of the columns in the data file and the
102 figure to be plotted.
103 -Sevelscale/confidence/fontsize.
106 Velocity ellipses in (N,E) convention. Vscale sets the scaling
107 of the velocity arrows. This scaling gives inches (unless c, i,
108 m, or p is appended). Confidence sets the 2-dimensional confi‐
109 dence limit for the ellipse, e.g., 0.95 for 95% confidence
110 ellipse. Fontsize sets the size of the text in points. The
111 ellipse will be filled with the color or shade specified by the
112 -G option [default transparent]. The arrow and the circumfer‐
113 ence of the ellipse will be drawn with the pen attributes speci‐
114 fied by the -W option. Parameters are expected to be in the
115 following columns:
116 1,2 longitude, latitude of station (-: option interchanges order)
118 3,4 eastward, northward velocity (-: option interchanges order)
120 5,6 uncertainty of eastward, northward velocities (1-sigma) (-:
122 option interchanges order)
123 7 correlation between eastward and northward components
125 8 name of station (optional).
127 -Snbarscale.
130 Anisotropy bars. Barscale sets the scaling of the bars This
131 scaling gives inches (unless c, i, m, or p is appended). Param‐
132 eters are expected to be in the following columns:
133 1,2 longitude, latitude of station (-: option interchanges order)
135 3,4 eastward, northward components of anisotropy vector (-: option
137 interchanges order)
138 -Srvelscale/confidence/fontsize
141 Velocity ellipses in rotated convention. Vscale sets the scal‐
142 ing of the velocity arrows. This scaling gives inches (unless c,
143 i, m, or p is appended). Confidence sets the 2-dimensional con‐
144 fidence limit for the ellipse, e.g., 0.95 for 95% confidence
145 ellipse. Fontsize sets the size of the text in points. The
146 ellipse will be filled with the color or shade specified by the
147 -G option [default transparent]. The arrow and the circumfer‐
148 ence of the ellipse will be drawn with the pen attributes speci‐
149 fied by the -W option. Parameters are expected to be in the
150 following columns:
151 1,2 longitude, latitude, of station (-: option interchanges order)
153 3,4 eastward, northward velocity (-: option interchanges order)
155 5,6 semi-major, semi-minor axes
157 7 counter-clockwise angle, in degrees, from horizontal axis to
159 major axis of ellipse.
160 8 name of station (optional)
162 -Swwedge_scale/wedge_mag.
165 Rotational wedges. Wedge_scale sets the size of the wedges in
166 inches (unless c, i, m, or p is appended). Values are multi‐
167 plied by Wedge_mag before plotting. For example, setting
168 Wedge_mag to 1.e7 works well for rotations of the order of 100
169 nanoradians/yr. Use -G to set the fill color or shade for the
170 wedge, and -E to set the color or shade for the uncertainty.
171 Parameters are expected to be in the following columns:
172 1,2 longitude, latitude, of station (-: option interchanges order)
174 3 rotation in radians
176 4 rotation uncertainty in radians
178 -Sxcross_scale
181 gives Strain crosses. Cross_scale sets the size of the cross in
182 inches (unless c, i, m, or p is appended). Parameters are
183 expected to be in the following columns:
184 1,2 longitude, latitude, of station (-: option interchanges order)
186 3 eps1, the most extensional eigenvalue of strain tensor, with
188 extension taken positive.
189 4 eps2, the most compressional eigenvalue of strain tensor, with
191 extension taken positive.
192 5 azimuth of eps2 in degrees CW from North.
194
197 No space between the option flag and the associated arguments.
198 -A Arrow_width/Head_length/Head_width Size of arrow in inches.
200 [Default is 0.03/0.12/0.09].
201 -B Sets map boundary annotation and tickmark intervals; see the
203 psbasemap man page for all the details.
204 -C Offset focal mechanisms to the latitude and longitude specified
206 in the last two columns of the input file.
207 -D Sigma_scale can be used to rescale the uncertainties of veloci‐
209 ties (-Se and -Sr) and rotations (-Sw). Can be combined with
210 the confidence variable.
211 -Ffill Sets the color or shade used for frame and annotation. [Default
213 is 0/0/0 (black)]
214 -Efill Sets the color or shade used for filling uncertainty wedges
216 (-Sw) or velocity error ellipses (-Se or -Sr). [If -E is not
217 specified, the uncertainty regions will be transparent.]
218 -Gfill Select filling of ellipses, wedges, and focal mechanisms. By
220 convention, the compressional quadrants of the focal mechanism
221 beach balls are shaded. Set the shade (0-255) or color (r/g/b)
222 [Default is 0/0/0]. Optionally, specify -Gpicon_size/pattern,
223 where pattern gives the number of the image pattern (1-90) OR
224 the name of a icon-format file. icon_size sets the unit size in
225 inches. To invert black and white pixels, use -GP instead of
226 -Gp. See pspatterns for information on individual patterns.
227 -H Input file(s) has header record(s). If used, the default number
229 of header records is N_HEADER_RECS. Use -Hi if only input data
230 should have header records [Default will write out header
231 records if the input data have them]. Blank lines and lines
232 starting with # are always skipped.
233 -K More PostScript code will be appended later [Default terminates
235 the plot system].
236 -L Draw lines. Ellipses and fault planes will have their outlines
238 drawn using current pen (see -W).
239 -N Do NOT skip symbols that fall outside the frame boundary speci‐
241 fied by -R.
242 [Default plots symbols inside frame only].
243 -O Selects Overlay plot mode [Default initializes a new plot sys‐
245 tem].
246 -P Selects Portrait plotting mode [Default is Landscape, see gmtde‐
248 faults to change this].
249 -U Draw Unix System time stamp on plot. By adding just/dx/dy/, the
251 user may specify the justification of the stamp and where the
252 stamp should fall on the page relative to lower left corner of
253 the plot. For example, BL/0/0 will align the lower left corner
254 of the time stamp with the lower left corner of the plot.
255 Optionally, append a label, or c (which will plot the command
256 string.). The GMT parameters UNIX_TIME, UNIX_TIME_POS, and
257 UNIX_TIME_FORMAT can affect the appearance; see the gmtdefaults
258 man page for details. The time string will be in the locale set
259 by the environment variable TZ (generally local time).
260 -V Selects verbose mode, which will send progress reports to stderr
262 [Default runs "silently"].
263 -W Set pen attributes for velocity arrows, ellipse circumference
265 and fault plane edges.
266 [Defaults: width = 1, color = 0/0/0, texture = solid].
267 -X -Y Shift plot origin relative to the current origin by (x-shift,y-
269 shift) and optionally append the length unit (c, i, m, p). You
270 can prepend a to shift the origin back to the original position
271 after plotting, or prepend r [Default] to reset the current
272 origin to the new location. If -O is used then the default (x-
273 shift,y-shift) is (0,0), otherwise it is (r1i, r1i) or (r2.5c,
274 r2.5c). Alternatively, give c to align the center coordinate (x
275 or y) of the plot with the center of the page based on current
276 page size.
277 -: Toggles between (longitude,latitude) and (latitude,longitude)
279 input and/or output. [Default is (longitude,latitude)]. Append
280 i to select input only or o to select output only. [Default
281 affects both].
282 -c Specifies the number of plot copies. [Default is 1].
284
286 The following should make big red arrows with green ellipses, outlined
287 in red. Note that the 39% confidence scaling will give an ellipse
288 which fits inside a rectangle of dimension Esig by Nsig.
289 psvelo << END -H2 -R-10/10/-10/10 -W0.25p,red -Ggreen -L -Se0.2/0.39/18
291 -B1g1 -Jx0.4/0.4 -A0.1/0.3/0.3 -P -V >! test.ps
292 Long. Lat. Evel Nvel Esig Nsig CorEN SITE
293 (deg) (deg) (mm/yr) (mm/yr)
294 0. -8. 0.0 0.0 4.0 6.0 0.500 4x6
295 -8. 5. 3.0 3.0 0.0 0.0 0.500 3x3
296 0. 0. 4.0 6.0 4.0 6.0 0.500
297 -5. -5. 6.0 4.0 6.0 4.0 0.500 6x4
298 5. 0. -6.0 4.0 6.0 4.0 -0.500 -6x4
299 0. -5. 6.0 -4.0 6.0 4.0 -0.500 6x-4
300 END
301 This example should plot some residual rates of rotation in the Western
302 Transverse Ranges, California. The wedges will be dark gray, with
303 light gray wedges to represent the 2-sigma uncertainties.
304 psvelo <<END -Sw0.4/1.e7 -W0.75p -Gdarkgray -Elightgray -H1 -D2 -Jm2.2
306 -R240./243./32.5/34.75 -Bf10ma60m/WeSn -P >! test.ps
307 lon lat spin(rad/yr) spin_sigma (rad/yr)
308 241.4806 34.2073 5.65E-08 1.17E-08
309 241.6024 34.4468 -4.85E-08 1.85E-08
310 241.0952 34.4079 4.46E-09 3.07E-08
311 241.2542 34.2581 1.28E-07 1.59E-08
312 242.0593 34.0773 -6.62E-08 1.74E-08
313 241.0553 34.5369 -2.38E-07 4.27E-08
314 241.1993 33.1894 -2.99E-10 7.64E-09
315 241.1084 34.2565 2.17E-08 3.53E-08
316 END
317
320 GMT(1), psbasemap(1), psxy(1)
321
323 Bomford, G., Geodesy, 4th ed., Oxford University Press, 1980.
324
326 Kurt Feigl
327 CNRS UMR 5562
328 Toulouse, France
329 (Kurt.Feigl@.cnes.fr)
330 Genevieve Patau
331 CNRS UMR 7580
332 Seismology Dept.
333 Institut de Physique du Globe de Paris
334 (patau@ipgp.jussieu.fr)
335GMT 4.5.6 10 Mar 2011 PSVELO(1)