Usage: taup setsac [--debug] [--evdpkm] [--geodetic] [--help] [--verbose]
                   [--geodeticflattening=f] [--prop=<toolPropsFilename>]
                   [[--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:

      --debug        enable debugging output
      --evdpkm       sac depth header is in km, default is meters
      --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.
      --geodeticflattening=f
                     Inverse Elliptical flattening for distance calculations
                       when --geodetic, defaults to WGS84 ~ 298.257. The
                       distance calculation uses 1/x.
      --help         display this help message
      --prop=<toolPropsFilename>
                     load defaults from properties file
      --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
