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  argu-
     ments.

     -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  compres-
     sional  quadrants.   Set  the shade (0-255) or color (r/g/b)
     [Default is 0/0/0].

     -g[<pen attributes>] Outline symbols in  compressional  qua-
     drants.  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)