30.3 Calculating of local phases: TTLOCAL

Program to calculate travel times for given stations and hypocenter location at local distances. The program uses the HYPOCENTER travel time calculations. By default, the model used is from STATION0.HYP, unless another model is selected in the s-file.

The program operates in two modes:


If iasp.inp is present, and there are no arguments, arrival times for stations in iasp.inp are calculated and put in file iasp.out to be used with MULPLT. The iasp.inp contains the stations in the waveform file. It is thus possible to calculate theoretical arrival times without any picks in the S-file. The iasp.inp is generated by MULPLT when the user asks for synthetic phases (MULPLT IASP option). The phases calculated are first and possibly second arrivals. For short distances there would only be first arrivals like Pg. Although refracted arrivals might be present, they are not calculated since they are unlikely to be seen. For larger distances both refracted P and PG are calculated like PN3 and PG. PN and PB are not specifically calculated but will appear as one of the PNX. When running the program, an output is also give on the text screen:

           1 FNA     230.977905     N4   35.8533058
           1 FNA     230.977905     G    38.0395508
           2 IGT     308.198364     N4   45.9507637
           2 IGT     308.198364     G    50.4460907
           3 MRVN    293.900574     N4   44.1098595
           3 MRVN    293.900574     G    48.1656075
           4 PUK     46.8152657     G    8.61053467
           5 SCTE    248.730225     N4   38.1072121
           5 SCTE    248.730225     G    40.8507156
           6 TIR     104.178276     G    17.6073112
           7 KBN     212.190460     N4   33.3805771
           7 KBN     212.190460     G    35.0128403
           8 PHP     106.046379     G    17.9404621
           9 SGRT    303.668915     N4   45.3243828
           9 SGRT    303.668915     G    49.6978455

This gives the station, epicentral distance, type of phase and travel time for the P-phase. The travel times for the S-phases are calculated using the Vp/VS ratio given in the station file and NOT by using a possible S-velocity model in the station file.

Free standing program for calculation of theoretical arrival times for station given in an S-file To operate in this mode, the command ttlocal sfile is given. The program calculates theoretical arrivals for the phases given in the s-file and write the output file with the original arrival times replaced by the theoretical times. All weights etc are kept. N and B-phases are used if given in S-file and specified in station file. A possible S-velocity model is ignored. Relocating the output file should give same locations as given in input file with zero residuals. In this mode there can be many events in the input s-file. The station file used is the one from the first event. Optionally the arrival times can be have added noise if a second argument is used like e.g. ttlocal collect.out 2.0 where input file is collect out and random noise is added. The noise generator is the Fortran RAND function which gives random numbers from 0-1 so in the example we would get random numbers between 0 and 2. The nose added to the arrival times for P is then between -1 and 1. The noise added to the S is multiplied with Vp/Vs.

The output file is ttlocal.out

Peter Voss : Tue Jun 8 13:38:42 UTC 2021