pssac1gmt

Plot seismograms in sac format on maps

Usage

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]

Overview

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.

Required arguments

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

saclist

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 columns:

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.

-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 [±]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 (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 dimension.

Optional arguments

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

-C[t0/t1]: 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.

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

-Ea|b|k|d|n[n]|u[n]: 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.

-F[i][q][r]

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.

-G[p|n][+gfill][+zzero][+tt0/t1]: 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.

-Msize[u][/alpha]: 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 = size*(north-south)/(depmax-depmin)/map_height.

size/alpha:

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 km.

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

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

-Q: Plot traces vertically.

-S[i]scale[unit]: 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.

-T[+tn][+rreduce_vel][+sshift]: Time alignment and shift.

+ttmark: align all trace along time mark. tmark are -5(b), -4(e), -3(o), -2(a), 0-9(t0-t9).

+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].

-X[a|c|f|r][x-shift[u]]

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

Examples

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 > single.ps

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

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

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 > map.ps
saclst stlo stla f *.z | gmt psxy -J -R -St0.4c -Gblack -i1,2 -O >> map.ps

References

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