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!