TauP Path

TauP Path uses a model like TauP Time but generates the angular distances from the epicenter at which the specified rays pierce path that the phases travel. The output is in GMT~cite{GMT} ``psxy’’ format, and is placed into the file taup_path.gmt. 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 path --mod iasp91 -h 550 --deg 74 -p S,ScS,sS,sScS --gmt

begin{verbatim} sh taup_path.gmt end{verbatim}

and you have a plot of the ray paths in TauP path.pdf. To avoid possible plotting errors for phases like Sdiff, the ray paths are interpolated to less than 1 degree increments.

The usage is:

Usage: taup path [arguments]
  or, for purists, java edu.sc.seis.TauP.TauP_Path [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

Distance is given by:

-deg degrees       -- distance in degrees,
-km kilometers     -- distance in kilometers,
                      assumes radius of earth is 6371km,

or by giving the station and event latitude and lonitude,
                      assumes a spherical earth,

-sta[tion] lat lon -- sets the station latitude and longitude
-evt       lat lon -- sets the event latitude and longitude

or by giving the takeoff angle,
--takeoff angle    -- takeoff angle from the source
                      zero is down, 90 horizontal, 180 is up

or by giving the ray parameter,
--shootray param   -- ray parameter from the source in s/deg
                      up or down is determined by the phase


--first            -- only output the first arrival for each phase, no triplications
--withtime        -- include time for each path point
--gmt             -- outputs path as a complete GMT script.
--svg             -- outputs path as a complete SVG file.
--mapwidth        -- sets map width for GMT script.

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