Task: Compute a Single Travel Time
To compute the travel time for a single phase and a given repicentral distance,
do the following:
- read the velocity model from an .nd file using model=mkreadnd(pfad);
- 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
-
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.