-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathgtt91.1
executable file
·126 lines (126 loc) · 3.9 KB
/
gtt91.1
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
.TH GTT91
.SH NAME
gtt91 \- construction of computer accessible travel time tables
.SH SYNOPSIS
.B gtt91
.SH DESCRIPTION
.B gtt91
comprises a sequence of routines to set up a representation of
the travel travel time tables for a given velocity model in terms of
the tau-spline procedure of Buland & Chapman (1983), together with
routines for user access (interactive and graphical).
.SS Construction of tau-spline tables:
.br
The two basic stages are:
.br
.I remodl
\- construction of partial tau integrals for different
segments of a velocity model of the Earth.
This requires discretisation of the velocity model at
depths corresponding to the slowness employed in the
representation of the tables. It is at this point in the
process that the velocity model has to be specified and
linked in with the rest of the routines - subsequently all
necessary information on the model is carried in the
direct access table files. The procedure remodl
generates two table files:
.br
.I remodl.hed, remodl.tbl
.br
which are used as the input to the second stage (setbrn)
.LP
.I setbrn
\- construction of the tau-spline coefficients for different
travel time branches corresponding to distinct seismic
phases. The tau-splines are generated for each of the
depths at which the model was sampled. In addition
auxiliary tables are set up to allow interpolation to
arbitrary depth. The procedure setbrn generates two
direct access table files identified by the name of the
velocity model
.br
e.g.
.I iasp91.hed, iasp91.tbl
.br
These tables are used in the routine
.I ttimes
which returns
travel times for seismic phases at specified distance
and in other display and location procedures
.LP
A Makefile template is provided to allow the construction of the
tau-spline tables for different velocity models. The specification
of the velocity model used for iasp91 is in terms of polynomials
(up to cubic) in normalised radius for different portions of the
velocity model.
The main body of routines used by the gtt91 system are stored in the
.br
.I $GTT91/libtau.f
.br
and two include files
.br
.I $GTT91/limits.inc and $GTT91/ttlim.inc
.br
are needed.
.SS Access to travel time tables:
.br
The simplest mode of access to the travel time tables is interactive
use of the
.I ttimes
command which allows user specification of
phases and returns the travel times for all specified phases at
a particular epicentral distance, for a given source depth.
.LP
The travel time information can also be accessed in graphical form
(Postscript output suitable for Laserwriter II) using the
zstp command which allows windowing in range and time.
A X-window version (zstx) is under development.
.LP
The interface to the tau-spline tables is readily adapted
to other computational procedures:
Printed travel time tables can be generated in a two stage process:
.br
i) Set up a summary file of travel timesfor many different source
depths and distances using the routine
.br
.I ttable
.br
which allows interactive specification of source depths
and distances, and includes all available phases.
.br
ii) Run the table sort routine
.br
.I wtac
.br
which allows the interrogation of the summary file produced by
.I ttable.
.br
.I wtac
prompts for a phase code and the number of characters to match and
also the range of epicentral distances required. Output of travel time
tables is to a specified file in
.I tbl
format. This
.I tbl
file can then be
printed using
.I ptroff
facilities.
.LP
Other facilities (including information on take-off angles
for source mechanism studies) are under development.
.SH ENVIRONMENT
.IP GTT91
Home directory for sources of gtt91 routines
access required for constructing modified
travel time tables
.SH SEE ALSO
ttimes
.SH AUTHORS
R. Buland,
National Earthquake Information Centre, U.S. Geological Survey
Golden, Colorado
.br
B.L.N. Kennett,
Research School of Earth Sciences, Australian National University
Canberra, Australia