manpagez: man pages & more
man psvelo(1)
Home | html | info | man
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
© manpagez.com 2000-2024
Individual documents may contain additional copyright information.