NAME

     psrose - Plot (length, azimuth) as windrose diagram or polar
     histogram (sector or rose diagram).

SYNOPSIS

     psrose file [ -Asector_width[r] ] [ -Btickinfo ] [ - C[mode-
     file]  ]  [  - Eazimuth/elevation ] [ -Fred/green/blue ] [ -
     Gfill ] [ -H ] [ -I ] [ -K ] [ -M[parameters ] [ -O ] [ -P ]
     [   -  Rr0/r1/az_0/az_1   ]   [  - Sradial_scale[n]  ]  [  -
     U[/dx/dy/][label] ] [ -V ] [ -Wpen ] [ -Xx-shift ]  [  - Yy-
     shift ] [ -Zscale ] [ -ccopies ] [ -: ]

DESCRIPTION

     psrose reads (length,azimuth) pairs from file  [or  standard
     input]  and  generates PostScript code that will plot a win-
     drose diagram.  Optionally (with -A), polar  histograms  may
     be  drawn (sector diagram or rose diagram).  Options include
     full circle and half circle plots.  The PostScript  code  is
     written to standard output.

     file Name of data file.  If no file is  given,  psrose  will
          read standard input.

OPTIONS

     No space between the option flag and  the  associated  argu-
     ments.

     -A    Gives the sector width in degrees for sector and  rose
          diagram.  [Default 0 means windrose diagram].  Append r
          to draw rose diagram instead of sector diagram.

     -B    Sets map boundary tickmark  intervals.  See  psbasemap
          for details.  Remember that "x" here is radial distance
          and "y" is azimuth.  The ylabel may be used to  plot  a
          figure caption.

     -C    Plot vectors showing the principal directions given in
          the  modes file.  If no file is given, compute and plot
          mean direction.

     -E    Sets the viewpoint's azimuth and elevation [180/90]

     -G    Selects color (r/g/b, each 0-255)  for  sector  infill
          [Default is black(0/0/0)].

     -H    Input file(s) has Header record(s).  Number of  header
          records  can  be  changed  by editing your .gmtdefaults
          file.  If used, GMT default is 1 header record.

     -I    Inquire.  Computes statistics needed to specify useful
          -R.  No plot is generated.

     -K    More PostScript code will be appended  later  [Default
          terminates the plot system].

          -     M    Specify      new      arrow       attributes
          tailwidth/headlength/headwidth/r/g/b   to   change  the
          appearance of arrows (Only if -C is set).

     -O    Selects Overlay plot mode [Default initializes  a  new
          plot system].

      - P    Selects  Portrait  plotting  mode  [GMT  Default  is
          Landscape, see gmtdefaults to change this].

     -R    Specifies the  'region'  of  interest  in  (r,azimuth)
          space.   r0  is  0,  r1  is  max  length in units.  For
          azimuth, specify -90/90 for half circle plot  or  0/360
          for full circle.

     -S    Specifies radius of circle in inch.  Append n to  nor-
          malize input radii to go from 0 to 1.

     -U    Draw Unix System time stamp on plot.  User may specify
          where the lower left corner of the stamp should fall on
          the page relative to lower left corner of plot in  inch
          [Default   is  (-0.75,-0.75)].   Optionally,  append  a
          label, or c (which will plot the command string.)

     -V    Selects verbose mode, which will send progress reports
          to stderr [Default runs "silently"].

     -W    Set pen attributes for sector outline  or  rose  plot.
          [Default is no outline].

     -X -Y
          Shift  origin  of  plot   by   (x-shift,y-shift)   inch
          [Default is (1,1) for new plots, (0,0) for overlays].

     -Z    Multiply the data radii by scale.  E.g., use  - Z0.001
          to  convert your data from m to km [Default is no scal-
          ing].

     -c    Specifies the number of plot copies. [Default is 1]

     -:    Input file has (azimuth,radius) pairs rather than  the
          expected (radius,azimuth).

EXAMPLES

     To plot a half circle rose diagram of the data in  the  file
     fault_segments.az_r (containing pairs of (azimuth, length in
     meters), using a 10 degree bin sector width, on a circle  of
     radius  = 3 inch, grid going out to radius = 150 km in steps
     of 25 km with a 30 degree sector interval, radial  direction
     annotated   every   50   km,  using  a  light  blue  shading
     (200/240/255) outlined by a solid red pen (width = 3),  draw
     the mean azimuth, and shown in Portrait orientation, try:

     psrose fault_segments.az_r  - R0/150/-90/90  - B50g25:"Fault
     length":/g30:."Rose  diagram":   - S3  -A10r -G200/240/255 -
     W3/255/0/0 -Z0.001 -C -P -: | lpr

     To plot a full circle wind rose diagram of the data  in  the
     file  lines.r_az, on a circle of radius = 2 inch, grid going
     out to radius = 500 units in steps of 100 with a  45  degree
     sector interval, using a solid pen (width = 2), and shown in
     landscape [Default] orientation with UNIX timestamp and com-
     mand line plotted, try:

     psrose lines.az_r -R0/500/0/360 - S2  - Bg100/g45:."Windrose
     diagram": -W2 -Uc | lpr

BUGS

     No default radial scale and grid settings for  polar  histo-
     grams.  User must run psrose -I to find max length in binned
     data set.

SEE ALSO

     gmt, gmtdefaults, pshistogram

REFERENCES

     Wessel, P., and W. H. F. Smith, 1995,  The  Generic  Mapping
     Tools  (GMT)  version  3.0  Technical  Reference & Cookbook,
     SOEST/NOAA.
     Wessel, P., and W. H. F. Smith, 1995,  New  Version  of  the
     Generic Mapping Tools Released, EOS Trans. AGU, 76, p. 329.
     Wessel, P., and W. H. F. Smith, 1995,  New  Version  of  the
     Generic           Mapping           Tools          Released,
     http://www.agu.org/eos_elec/95154e.html, Copyright  1995  by
     the American Geophysical Union.
     Wessel, P., and W. H. F. Smith, 1991,  Free  Software  Helps
     Map and Display Data, EOS Trans. AGU, 72, p. 441.