## NAME

          pspolar - Plot polarities on the inferior focal half-sphere on maps

## SYNOPSIS

          pspolar files -Jparameters -Rwest/east/south/north [
          -Dlongitude/latitude -Msize -Ssymbol/size [ -Btickinfo ] [
          -Clongitude/latitude [/tiret_width/pointsize ] ] [ -Fred/green/blue ] [
          -Gfill ] [ -gfill ] [ -H ] [ -h ] [ -K ] [ -L ] [ -N ] [ -O ] [ -P ] [
          -sHalf-size[v[[<v_width/h_length/h_width/shape>]][g[<r/g/b>]][l] [
          -U[label] ] [ -V ] [ -Wpen ] [ -Xx-shift ] [ -Yy-shift ] [ -ccopies ]

## DESCRIPTION

          pspolar reads data values from files [or standard input] and generates
          PostScript code that will plot stations on focal mechanisms on a map.
          The PostScript code is written to standard output.

## ARGUMENTS

          files List one or more file-names. If no files are given, pspolar will
          read standard input.

          -J Selects the map projection. Scale is cm/degree, 1:xxxxx. or width
          in cm (upper case modifier).

          See psbasemap manpages for details.

          -R west, east, south, and north specify the Region of interest.  To
          specify your boundaries in degrees and minutes, use the dd:mm format.

          -Dlongitude/latitude Maps the bubble at given longitude and latitude
          point.

          -Msize Sets the size of the beach ball to plot polarities in.
          Size is in inch (unless c, i, m, or p is appended).

          -S<symbol type>sizeSelects symbol type and symbol size .
          Size is in inch (unless c, i, m, or p is appended).
          Choose symbol type between
             st(a)r, (c)ircle, (d)iamond, (h)exagon, (i)nvtriangle
          (p)oint, (s)quare, (t)riangle, (x)cross

## OPTIONS

          No space between the option flag and the associated arguments.

          -B Sets map boundary tickmark intervals. See psbasemap for details.

          -C Offsets focal mechanisms to the latitude and longitude specified in
          the last two columns of the input file.

          -Efill Selects filling of symbols for stations in extensive quadrants.
          Set the shade (0-255) or color (r/g/b) [Default is 250/250/250].  If
          Efill is the same as Ffill, use -e to outline.

          -e[<pen attributes>] Outline symbols in extensive quadrants.
          Default pen attributes is W pen attributes.

          -Ffill Sets background color of the beach ball. Default is no fill.

          -f[<pen attributes>] Outline the beach ball. Default pen attributes is
          W pen attributes.

          -Gfill Selects filling of symbols for stations in compressional
          quadrants.  Set the shade (0-255) or color (r/g/b) [Default is 0/0/0].

          -g[<pen attributes>] Outline symbols in compressional quadrants.
          Default pen attributes is W pen attributes.

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

          -h
          Use special format derived from HYPO71 output

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

          -N Does NOT skip symbols that fall outside map border [Default plots
          points inside border only].

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

          -P Selects Portrait plotting mode. [Default is Landscape]

          -shalf-size/[V[<v_width/h_length/h_width/shape>]][G<r/g/b>][L]
          Plots S polarity azimuth.
          S polarity is in last column.  It may be a vector (V option) or a
          segment. Give half-size,v_width,h_length,h_width in inch (unless c, i,
          m, or p is appended).  [L] option is for outline.

          -Tangle/form/justify/fontsize in points
          To write station code. Default is 0.0/0/5/12.

          -t<pen attributes>
          To write station code. Default attributes is W attributes.

          -U Draws 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 current pen attributes
           [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]

          Parameters are expected to be in the following columns

          1,2,3,4
               station_code, azimuth, take-off angle, polarity

               polarity :
               - compression can be c,C,u,U,+
               - rarefaction can be d,D,r,R,-
               - not defined is anything else

## EXAMPLES

          pspolar << END -R239/240/34/35.2 -JM8 -N -Sc0.4 -H1 -D239.5/34.5 -M5
          >! test.ps
          stat azim  ih pol
          0481   11  147   c
          6185  247  120   d
          0485  288  114   +
          0490  223  112   -
          0487  212  109   .
          END
          or
          pspolar << END -R239/240/34/35.2 -JM8 -N -Sc0.4 -H1 -D239.5/34.5 -M5
          -h >! test.ps
          Date   Or. time stat azim  ih
          910223   1 22   0481   11  147   ipu0
          910223   1 22   6185  247  120   ipd0
          910223   1 22   0485  288  114   epu0
          910223   1 22   0490  223  112   epd0
          910223   1 22   0487  212  109   epu0
          END

## SEE ALSO

          gmtsystem, psbasemap, psxy.

## REFERENCES

          Bomford, G., Geodesy, 4th ed., Oxford University Press, 1980.
          Aki, K. and P. Richards, Quantitative Seismolgy, Freeman, 1980.

## AUTHORS

          Genevieve Patau
          Seismology Dept.
          Institut de Physique du Globe de Paris
          (patau@.ipgp.jussieu.fr)