Usage: taup setsac [--debug] [--evdpkm] [--geodetic] [--help] [--verbose]
                   [--az=<azimuth>] [--baz=<backAzimuth>]
                   [--equitorialradius=m] [--invflattening=f]
                   [--planet=<planet>] [--prop=<toolPropsFilename>] [--evt=lat
                   lon]... [--sta=lat lon]... [--geodist=<geoDistTypes>...]...
                   [[--mod=<modelName>] [-h=depth[,depth...]]...
                   [--stadepth=depth[,depth...]]... [--scat=depth deg]...]
                   [[-p=phase[,phase...]]... [--phasefile=<phaseFiles>]...]
                   [sacfile...]
Set headers to travel times of phases using depth and distance from SAC files.
      [sacfile...]           SAC files to process

Options:

      --az=<azimuth>         azimuth in degrees, source to receiver
      --baz=<backAzimuth>    backazimuth in degrees, receiver to source
      --debug                enable debugging output
      --equitorialradius=m   Equitorial radius in meters for distance
                               calculations when geodetic or geocentric,
                               defaults to WGS84 ~ 6378137 meters.
      --evdpkm               sac depth header is in km, default is meters
      --evt, --event=lat lon event latitude and longitude.  Creates a distance
                               if station is also given.
      --geodetic             use geodetic latitude for distance calculations,
                               which implies an ellipticity. Default is
                               spherical. Note this only affects calculation of
                               distance from lat/lon pairs, all travel time
                               calculations are done in a purely spherical
                               model.
      --geodist=<geoDistTypes>...
                             Type of distance calculation to use for lat,lon
                               distance calculation, one of spherical,
                               geocentric, geodetic. Default is spherical. Note
                               this only affects calculation of distance from
                               lat/lon pairs, all travel time calculations are
                               done in a purely spherical model.
      --help                 display this help message
      --invflattening, --geodeticflattening=f
                             Inverse Elliptical flattening for distance
                               calculations when geocentric or geodetic,
                               defaults to WGS84 ~ 298.257. The distance
                               calculation uses 1/x.
      --planet=<planet>      Geodesic for planets, for distance calculations
                               when --geodetic, one of: EARTH, MOON, MARS,
                               VENUS, MERCURY, CERES, JUPITER, SATURN, URANUS,
                               NEPTUNE.
      --prop=<toolPropsFilename>
                             load defaults from properties file
      --sta, --station=lat lon
                             station latitude and longitude. Creates a distance
                               if event is also given.
      --verbose              enable verbose output
Model Args
  -h, --evdepth, --sourcedepth=depth[,depth...]
                             source depth in km
      --mod, --model=<modelName>
                             use velocity model "modelName" for calculations.
                             Default is iasp91. Other builtin models include
                               prem, ak135, ak135fcont, and ak135favg.
      --scat, --scatter=depth deg
                             scattering depth and distance in degrees, which
                               may be negative. Only effects phases with 'o' or
                               'O' in the phase name.
      --stadepth, --receiverdepth=depth[,depth...]
                             the receiver depth in km for stations not at the
                               surface
Phase Names
  -p, --ph, --phase=phase[,phase...]
                             seismic phase names
      --phasefile=<phaseFiles>
                             read list of phase names from file
