Task: Compute a Single Travel Time

To compute the travel time for a single phase and a given repicentral distance, do the following:

  1. read the velocity model from an .nd file using model=mkreadnd(pfad);
  2. call [tt,p]=mkttime(phase,delta,h,model); where
    phase
    a string containing the phase name, e.g. 'P' or 'SKS'
    delta
    the epicentral distance, in degrees
    h
    the focal depth in km
    model
    a model structure containing the velocity structure
  3. As Result, the routine returns
    tt
    The requested travel time in seconds
    p
    The ray parameter in s/° which is needed to reach the given epicentral distance from the given focal depth.

NOTE that bot tt and p might be vectors in case of phase triplications!


eof.