TauP Wavefront

TauP Wavefront is similar to TauP Path, but plots the wavefront at timesteps instead of the ray paths. It also uses a model like TauP Time and generates depth versus angular distances from the epicenter for the phases, but done at time slices instead of depth slices. The output is in GMT~cite{GMT} psxy format, and is placed into the file taup_wavefront.gmt. The colums are distance, depth, time and ray param, although only the first two are used by the GMT script. If you specify the -gmt flag then this is a complete script with the appropriate psxy command prepended, so if you have GMT installed, you can just:

taup wavefront --mod iasp91 -h 550 -p S,ScS,sS,sScS --gmt

begin{verbatim} taup wavefront -mod iasp91 -h 550 -ph S,ScS,sS,sScS –gmt sh taup_wavefront.gmt end{verbatim}

and you have a plot of the wavefronts in taup_wavefront.pdf.

Or use –svg to generate a SVG plot

taup wavefront --mod iasp91 -h 550 -p S,ScS,sS,sScS --svg

The usage is:

Usage: taup wavefront [arguments]
  or, for purists, java edu.sc.seis.TauP.TauP_Wavefront [arguments]

Arguments are:
-ph phase list     -- comma separated phase list
-pf phasefile      -- file containing phases

-mod[el] modelname -- use velocity model "modelname" for calculations
                      Default is iasp91.

-h depth           -- source depth in km

--stadepth depth   -- receiver depth in km
--scat[ter] depth deg   -- scattering depth and distance

--gmt             -- outputs path as a complete GMT script.
--svg             -- outputs path as a complete SVG file.
--mapwidth        -- sets map width for GMT script.
--rays  num      -- number of raypaths/distances to sample.
--timestep  num  -- steps in time (seconds) for output.
--timefiles      -- outputs each time into a separate .ps file within the gmt script.
--negdist        -- outputs negative distance as well so wavefronts are in both halves.

-o [stdout|outfile]         -- output is redirected to stdout or to the "outfile" file
--prop [propfile]   -- set configuration properties
--debug             -- enable debugging output
--verbose           -- enable verbose output
--version           -- print the version
--help              -- print this out, but you already know that!