Provided by: gmt-common_5.4.5+dfsg-1_all bug


       pssac - Plot seismograms in SAC format on maps


       pssac [ saclist|SACfiles ]  -Jparameters
        -Rregion  [  -B[p|s]parameters ] [  -C[t0/t1] ] [  -Ddx[/dy] ] [  -Ea|b|k|d|n[n]|u[n] ] [
       -F[i][q][r] ] [  -G[p|n][+gfill][+zzero][+tt0/t1] ] [  -K ] [  -Msize[u][/alpha] ] [  -O ]
       [  -P ] [  -Q ] [  -S[i]scale[unit] ] [  -T[+tn][+rreduce_vel][+sshift] ] [  -U[stamp] ] [
       -V[level] ] [  -Wpen ] [  -Xx_offset ] [  -Yy_offset ] [ -hheaders ] [ -ttransp ]

       Note: No space is allowed between the option flag and the associated arguments.


       pssac reads SACfiles in SAC format or reads  filenames  and  controlling  parameters  from
       saclist  [or standard input] and generates PostScript that will plot seismograms on a map.
       The PostScript code is written to standard output.


              SAC files to plot on a map. Only evenly spaced SAC data is supported.

              One ASCII data table file holding a number of data columns. If saclist is not given
              then  we  read from standard input.  Parameters are expected to be in the following
                 filename [X Y [pen]]

              filename is the name of SAC file to plot.  X and Y are the position of  seismograms
              to  plot  on  a map.  On linear plots, the default X is the begin time of SAC file,
              which will be adjusted if -T option is used, the default  Y  is  determined  by  -E
              option.   On  geographic  plots,  the  default  X  and  Y are station longitude and
              latitude specified in SAC header.  The  X  and  Y  given  here  will  override  the
              position  determined by command line options.  pen, if given, will override the pen
              from -W option for current SAC file only.

       -Jparameters (more ...)
              Select map projection.

              west, east, south, and north specify the region of interest, and  you  may  specify
              them  in decimal degrees or in [±]dd:mm[][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 (Cartesian) 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


       -B[p|s]parameters (more ...)
              Set map boundary frame and axes attributes.

              Read  and  plot  seismograms  in  timewindow between t0 and t1 only.  t0 and t1 are
              relative to a reference time specified by -T.  If -T option is not  specified,  use
              the  reference  time (kzdate and kztime) defined in SAC header instead.  If only -C
              is used, t0/t1 is determined as xmin/xmax from -R option.

              Offset seismogram positions by the given mount dx/dy [Default is no offset].  If dy
              is not given it is set equal to dx.

              Choose profile type (the type of Y axis).
                 a: azimuth profile.

                 b: back-azimuth profile.

                 k: epicentral distance (in km) profile.

                 d: epicentral distance (in degree) profile.

                 n: trace number profile. The Y position of first trace is numbered as n [Default
                 n is 0].

                 u: user defined profile. The Y positions are determined by SAC  header  variable
                 usern, default using user0.

              Data preprocess before plotting.
                 i: integral

                 q: square

                 r: remove mean value

              i|q|r  can  repeat  multiple times. For example, -Frii will convert acceleration to
              displacement.  The order of i|q|r controls the order of the data processing.

              Paint positive or negative portion of traces.  If only -G is used, default to  fill
              the positive portion black.
                 p|n  controls  the  painting of positive portion or negative portion.  Repeat -G
                 option to specify fills for positive and negative portions, respectively.

                 +gfill: color to fill

                 +tt0/t1: paint traces between t0 and t1 only. The reference time of t0 and t1 is
                 determined by -T option.

                 +zzero:  define  zero  line.  From zero to top is positive portion, from zero to
                 bottom is negative portion.

       -K (more ...)
              Do not finalize the PostScript plot.

              Vertical scaling.
                 size[u]:  scale  all  traces  size[u]  on   a   map.   The   default   unit   is
                 PROJ_LENGTH_UNIT.     The    scaling    factor    is   defined   as   yscale   =

                     alpha < 0, use the same scaling factor for all traces.  The  scaling  factor
                     will scale the first trace to size[u].

                     alpha = 0, multiply all traces by size. No unit is allowed.

                     alpha  >  0, multiply all traces by size*r^alpha, r is the distance range in

       -O (more ...)
              Append to existing PostScript plot.

       -P (more ...)
              Select "Portrait" plot orientation.

       -Q     Plot traces vertically.

              Sets time scale in seconds per <unit> while plotting on geographic  plots.   Append
              c, i, or p to indicate cm, inch or points as the unit. Use PROJ_LENGTH_UNIT if unit
              is omitted.  Use -Siscaleunit to give the reciprocal scale, i.e. cm per  second  or
              inch per second.

              Time alignment and shift.
                 +ttmark:  align all trace along time mark. tmark are -5(b), -4(e), -3(o), -2(a),

                 +rreduce_vel: reduce velocity in km/s.

                 +sshift: shift all traces by shift seconds.

       -U[[just]/dx/dy/][c|label] (more ...)
              Draw GMT time stamp logo on plot.

       -V[level] (more ...)
              Select verbosity level [c].

       -Wpen  Set pen attributes for all traces unless overruled by  pen  specified  in  saclist.
              [Defaults: width = default, color = black, style = solid].


       -Y[a|c|f|r][y-shift[u]] (more ...)
              Shift plot origin.

       -h[i|o][n][+c][+d][+rremark][+rtitle] (more ...)
              Skip or produce header record(s).

       -t[transp] (more ...)
              Set PDF transparency level in percent.

       -^ 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  explanation  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.


       To plot a single seismogram seis.SAC (generated by SAC  command  funcgen  seismogram)  and
       paint positive portion black and negative portion red:

          gmt pssac seis.SAC -JX10c/5c -R9/20/-2/2 -Baf -Fr -Gp+gblack -Gn+gred >

       To  plot several seismograms (generated by SAC command datagen sub tele *.z) on a distance

          gmt pssac *.z -R200/1600/12/45 -JX15c/5c -Bx200+l'T(s)' -By5+lDegree -BWSen \
              -Ed -M1.5c -W0.5p,red >

       To plot seismograms (generated by SAC command datagen sub tele *.z) on a geographic map:

          gmt pssac *.z -JM15c -R-120/-40/35/65 -Baf -M1i -S300c -K >
          saclst stlo stla f *.z | gmt psxy -J -R -St0.4c -Gblack -i1,2 -O >>


       psmeca, pspolar, pscoupe, psvelo, gmt, psbasemap, psxy


       Refer to SAC User Manual for more details on SAC format and SAC header variables.


       Dongdong Tian, School of Earth and Space Sciences, University of Science and Technology of
       China, Hefei, Anhui, China


       2019, P. Wessel, W. H. F. Smith, R. Scharroo, J. Luis, and F. Wobbe