-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathREADME
120 lines (87 loc) · 3.45 KB
/
README
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
Originally Developed by: Jeroen Tromp
Modified by:
Qinya [email protected]
03/21/02
Hom Nath Gharti
xnmsyn calculates a normal mode synthetics using normal mode
summation for given CMTSOLUTION and STATIONS files.
1. INPUT PARAMTERS:
EXAMPLE---
./database/prem_noocean_model
./database/toroidal_prem_noocean.bin
./database/radial_prem_noocean.bin
./database/spheroidal_prem_noocean.bin
./input/CMTSOLUTION
./input/STATIONS
17.0
1 12000
are respectively: model deck file, toroidal catalogue, radial catalogue,
and spheroidal catalogue, CMTSOLUTION file name, STATIONS file name, shortest period
in summation, dt in output time series, and total length of record in
seconds.
example execution: xnmsyn input/prem_noocean.in
2. SETUP
Make sure you have output/ folder in the current directory.
file STATIONS has the format of
2
KMI IC 25.1233 102.7400 1975.0 0.0
PAS TS 34.1490 -118.1720 0.0 0.0
first line -- number of stations
second - last line -- station_name, network_code stalat, stalon, staele, burial
CMTSOLUTION files are the Harvard format:
PDE 1994 6 9 0 33 45.40 -13.8200 -67.2500 647.1 6.9 6.8 NORTHERN BOLIVIA
event name: M060994A
time shift: 0.0000
half duration: 30.0000
latitude: -13.8200
longitude: -67.2500
depth: 647.1000
Mrr: -0.759000E+28
Mtt: 0.775000E+28
Mpp: -0.160000E+27
Mrt: -0.250300E+29
Mrp: 0.420000E+27
Mtp: -0.248000E+28
3. OUTPUT
OUTPUT will be in output/ with names like:
PAS.TS.LHE.nmsyn PAS.TS.LHN.nmsyn PAS.TS.LHZ.nmsyn
and the old version of synthetics with the same name will
be overwritten, so be sure to rename the old ones if they
are still useful.
You can type
xmgr -graph 0 PAS.TS.LHE.nmsyn -graph 1 PAS.TS.LHN.nmsyn -graph 2 PAS.TS.LHZ.nmsyn -rows 3 to see preliminary result.
REMARKS:
1. If you want to make synthetics with another mode catalogue, you might
want to change the parameters in minos.h, including NR and NL. Make sure
you find out the written format of the catalogue and change correspondingly
in the nmsyn.f.
2. If you want output to be RTZ component instead of NEZ component, you can
change parts of write_record.f, or do the rotation with your own program.
3. A minor problem exists for the output synthetics. The synthetics is not
perfectly zero before P arrivals, and I cannot figure out what causes this.
You are more than welcomed to find out why. ( Jeroen has a program to remove
the offset)
4. There are several mode catalogues available for summation:
radial.aniprem808_40s.bin
toroidal.aniprem808_40s.bin
spheroidal.aniprem808_40s.bin
is a 40s mode catalogue for aniprem808 with top 165 layers stored.
radial.aniprem808_8s.bin
toroidal.aniprem808_8s.bin
spheroidal.aniprem808_8s.bin
is a 8s mode catalogue for aniprem808 with top 165 layers stored.
Use the 40s catalogue directly instead of using tmin=40 for the 8s
catalogue will save you a lot of time.
5. One can compare the output with the output of QmXD in Sun
workstation, using the following script
#!/usr/bin/tcsh -f
setenv MYDBS /net/freawaru/export/home2/jtromp/dbs
/net/freawaru/export/home2/jtromp/bin/QmXD -M \
/net/freawaru/export/home1/jtromp/dbs/catalogues/seigsml_prem_an.750.8s \
-O 8. -P 6000. -s
The output generally agrees farely well except the problem I described
earlier.
6. The dt can be chosen as long as 2*dt < Tmin, larger dt will result
in faster computation.
7. The program was originally written by Jeroen, and I only changed the
input and output format.