Skip to content

Commit d17ed08

Browse files
author
homnath
committed
- first commit
0 parents  commit d17ed08

Some content is hidden

Large Commits have some content hidden by default. Use the searchbox below for content that may be hidden.

51 files changed

+3199
-0
lines changed

‎.gitignore

+6
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,6 @@
1+
SYNT*/
2+
*.pdf
3+
*.eps
4+
*.bak
5+
xnmsyn*
6+
*.mod

‎CMTSOLUTION

+13
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,13 @@
1+
PDE 2011 3 11 5 46 23.00 38.3200 142.3700 24.4 7.9 7.9 NEAR EAST COAST OF HONSH
2+
event name: 201103110546A
3+
time shift: 69.7600
4+
half duration: 70.0000
5+
latitude: 37.5200
6+
longitude: 143.0500
7+
depth: 20.0000
8+
Mrr: 1.730000e+29
9+
Mtt: -2.810000e+28
10+
Mpp: -1.450000e+29
11+
Mrt: 2.120000e+29
12+
Mrp: 4.550000e+29
13+
Mtp: -6.570000e+28

‎CMTSOLUTION-china

+13
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,13 @@
1+
PDE 1994 6 9 0 33 16.40 -13.8300 -67.5600 637.0 6.9 6.8 NORTHERN BOLIVIA
2+
event name: 060994A
3+
time shift: 0.0000
4+
half duration: 0.0000
5+
latitude: 25.3900
6+
longitude: 101.4000
7+
depth: 33.0000
8+
Mrr: -0.600000e+24
9+
Mtt: -6.290000e+24
10+
Mpp: 6.890000e+24
11+
Mrt: -1.850000e+24
12+
Mrp: 0.120000e+24
13+
Mtp: -4.730000e+24

‎CMTSOLUTIONSL

+14
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,14 @@
1+
PDE 1994 6 9 0 33 16.40 -13.8300 -67.5600 637.0 6.9 6.8 NORTHERN BOLIVIA
2+
event name: 060994A
3+
time shift: 0.0000
4+
half duration: 20.0000
5+
latitude: 89.9900
6+
longitude: 0.0000
7+
depth: 10.0000
8+
Mrr: 0.0000000e00
9+
Mtt: 0.0000000e00
10+
Mpp: 0.0000000e00
11+
Mrt: 0.0000000e00
12+
Mrp: 0.000000e+00
13+
Mtp: -2.500000e+27
14+

‎CMTSOLUTION_Bolivia

+14
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,14 @@
1+
PDE 1994 6 9 0 33 16.40 -13.8300 -67.5600 637.0 6.9 6.8 NORTHERN BOLIVIA
2+
event name: 060994A
3+
time shift: 0.0000
4+
half duration: 60.0000
5+
latitude: 89.9900
6+
longitude: 0.0000
7+
depth: 10.0000
8+
Mrr: -7.600000e+27
9+
Mtt: 7.700000e+27
10+
Mpp: -2.000000e+26
11+
Mrt: -2.500000e+28
12+
Mrp: 4.000000e+26
13+
Mtp: -2.500000e+27
14+

‎CMTSOLUTION_Tohoku

+13
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,13 @@
1+
PDE 2011 3 11 5 46 23.00 38.3200 142.3700 24.4 7.9 7.9 NEAR EAST COAST OF HONSH
2+
event name: 201103110546A
3+
time shift: 69.7600
4+
half duration: 70.0000
5+
latitude: 37.5200
6+
longitude: 143.0500
7+
depth: 20.0000
8+
Mrr: 1.730000e+29
9+
Mtt: -2.810000e+28
10+
Mpp: -1.450000e+29
11+
Mrt: 2.120000e+29
12+
Mrp: 4.550000e+29
13+
Mtp: -6.570000e+28

‎CMTSOLUTION_bol

+13
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,13 @@
1+
PDE 1994 6 9 0 33 16.40 -13.8300 -67.5600 637.0 6.9 6.8 NORTHERN BOLIVIA
2+
event name: 060994A
3+
time shift: 0.0000
4+
half duration: 0.0000
5+
latitude: 89.9900
6+
longitude: 0.0000
7+
depth: 10.0000
8+
Mrr: -7.600000e+27
9+
Mtt: 7.700000e+27
10+
Mpp: -2.000000e+26
11+
Mrt: -2.500000e+28
12+
Mrp: 4.000000e+26
13+
Mtp: -2.500000e+27

‎CMTSOLUTION_peru

+16
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,16 @@
1+
PDE 2001 7 7 9 38 43.50 -17.5400 -72.0800 33.0 6.6 7.3 NEAR COAST OF PERU
2+
event name: 070701F
3+
time shift: 18.3000
4+
half duration: 14.0000
5+
latitude: -17.4500
6+
longitude: -72.4500
7+
depth: 25.0000
8+
Mrr: 1.140000e+27
9+
Mtt: -2.100000e+26
10+
Mpp: -9.300000e+26
11+
Mrt: 6.800000e+26
12+
Mrp: -2.850000e+27
13+
Mtp: 7.700000e+26
14+
15+
16+

‎CMTSOLUTIONcor

+14
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,14 @@
1+
PDE 1994 6 9 0 33 16.40 -13.8300 -67.5600 637.0 6.9 6.8 NORTHERN BOLIVIA
2+
event name: 060994A
3+
time shift: 0.0000
4+
half duration: 60.0000
5+
latitude: 89.9900
6+
longitude: 0.0000
7+
depth: 10.0000
8+
Mrr: -7.600000e+27
9+
Mtt: 7.700000e+27
10+
Mpp: -2.000000e+26
11+
Mrt: -2.500000e+28
12+
Mrp: 4.000000e+26
13+
Mtp: -2.500000e+27
14+

‎Makefile

+11
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,11 @@
1+
FC = gfortran # fortran compiler
2+
FFLAGS = -O3 -funroll-loops -ffixed-line-length-132
3+
4+
xnmsyn: delaz.f90 get_cmt.f90 mnam.f90 plm.f90 timeseries.f90 nmsyn.f90 tdiffer.f90 write_record.f90 get_mode.f90 check_norm.f90 intgrl_disc.f90 find_disc.f90
5+
$(FC) $(FFLAGS) -g -o xnmsyn delaz.f90 get_cmt.f90 mnam.f90 plm.f90 nmsyn.f90 tdiffer.f90 timeseries.f90 write_record.f90 get_mode.f90 check_norm.f90 intgrl_disc.f90 find_disc.f90
6+
7+
xnmsyn2: delaz.f90 get_cmt.f90 mnam.f90 plm.f90 nmsyn2.f90 tdiffer.f90 timeseries.f90 write_record.f90 adr2cmt.f90
8+
$(FC) $(FFLAGS) -g -o xnmsyn2 delaz.f90 get_cmt.f90 mnam.f90 plm.f90 nmsyn2.f90 tdiffer.f90 timeseries.f90 write_record.f90 adr2cmt.f90
9+
10+
clean:
11+
\rm -f *.o xnmsyn *~

‎README

+115
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,115 @@
1+
2+
xnmsyn calculates a normal mode synthetics using normal mode
3+
summation for given CMTSOLUTION and STATIONS files.
4+
5+
1. INPUT PARAMTERS:
6+
7+
EXAMPLE---
8+
/home/lqy/catalogue/data/aniprem808.dk
9+
/home/lqy/catalogue/data/toroidal.aniprem808.bin
10+
/home/lqy/catalogue/data/radial.aniprem808.bin
11+
/home/lqy/catalogue/data/spheroidal.aniprem808.bin
12+
CMTSOLUTION
13+
8.0
14+
1 3600
15+
16+
are respectively: model deck file, toroidal catalogue, radial catalogue,
17+
and spheroidal catalogue, CMTSOLUTION file name, shortest period
18+
in summation, dt in output time series, and total length of record in
19+
seconds.
20+
21+
example execution: xnmsyn < in_8
22+
23+
2. SETUP
24+
25+
Make sure you have STATIONS and SYNT/ in the current directory.
26+
27+
file STATIONS has the format of
28+
29+
2
30+
KMI IC 25.1233 102.7400 1975.0 0.0
31+
PAS TS 34.1490 -118.1720 0.0 0.0
32+
33+
first line -- number of stations
34+
second - last line -- station_name, network_code stalat, stalon, staele, burial
35+
36+
37+
CMTSOLUTION files are the Harvard format:
38+
39+
PDE 1994 6 9 0 33 45.40 -13.8200 -67.2500 647.1 6.9 6.8 NORTHERN BOLIVIA
40+
event name: M060994A
41+
time shift: 0.0000
42+
half duration: 30.0000
43+
latitude: -13.8200
44+
longitude: -67.2500
45+
depth: 647.1000
46+
Mrr: -0.759000E+28
47+
Mtt: 0.775000E+28
48+
Mpp: -0.160000E+27
49+
Mrt: -0.250300E+29
50+
Mrp: 0.420000E+27
51+
Mtp: -0.248000E+28
52+
53+
3. OUTPUT
54+
55+
OUTPUT will be in SYNT/ with names like:
56+
PAS.TS.LHE.nmsyn PAS.TS.LHN.nmsyn PAS.TS.LHZ.nmsyn
57+
and the old version of synthetics with the same name will
58+
be overwritten, so be sure to rename the old ones if they
59+
are still useful.
60+
61+
You can type
62+
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.
63+
64+
65+
REMARKS:
66+
67+
1. If you want to make synthetics with another mode catalogue, you might
68+
want to change the parameters in minos.h, including NR and NL. Make sure
69+
you find out the written format of the catalogue and change correspondingly
70+
in the nmsyn.f.
71+
72+
2. If you want output to be RTZ component instead of NEZ component, you can
73+
change parts of write_record.f, or do the rotation with your own program.
74+
75+
3. A minor problem exists for the output synthetics. The synthetics is not
76+
perfectly zero before P arrivals, and I cannot figure out what causes this.
77+
You are more than welcomed to find out why. ( Jeroen has a program to remove
78+
the offset)
79+
80+
4. There are several mode catalogues available for summation:
81+
82+
radial.aniprem808_40s.bin
83+
toroidal.aniprem808_40s.bin
84+
spheroidal.aniprem808_40s.bin
85+
is a 40s mode catalogue for aniprem808 with top 165 layers stored.
86+
87+
radial.aniprem808_8s.bin
88+
toroidal.aniprem808_8s.bin
89+
spheroidal.aniprem808_8s.bin
90+
is a 8s mode catalogue for aniprem808 with top 165 layers stored.
91+
92+
Use the 40s catalogue directly instead of using tmin=40 for the 8s
93+
catalogue will save you a lot of time.
94+
95+
5. One can compare the output with the output of QmXD in Sun
96+
workstation, using the following script
97+
98+
#!/usr/bin/tcsh -f
99+
setenv MYDBS /net/freawaru/export/home2/jtromp/dbs
100+
/net/freawaru/export/home2/jtromp/bin/QmXD -M \
101+
/net/freawaru/export/home1/jtromp/dbs/catalogues/seigsml_prem_an.750.8s \
102+
-O 8. -P 6000. -s
103+
104+
The output generally agrees farely well except the problem I described
105+
earlier.
106+
107+
6. The dt can be chosen as long as 2*dt < Tmin, larger dt will result
108+
in faster computation.
109+
110+
7. The program was originally written by Jeroen, and I only changed the
111+
input and output format.
112+
113+
114+
Qinya [email protected]
115+
03/21/02

‎STATIONS

+23
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,23 @@
1+
22
2+
X1 Y5 -90 0 0.0 0.0
3+
X11 Y5 -80 0 0.0 0.0
4+
X21 Y5 -70 0 0.0 0.0
5+
X31 Y5 -60 0 0.0 0.0
6+
X41 Y5 -50 0 0.0 0.0
7+
X51 Y5 -40 0 0.0 0.0
8+
X61 Y5 -30 0 0.0 0.0
9+
X71 Y5 -20 0 0.0 0.0
10+
X81 Y5 -10 0 0.0 0.0
11+
X91 Y5 0 0 0.0 0.0
12+
X101 Y5 10 0 0.0 0.0
13+
X111 Y5 20 0 0.0 0.0
14+
X121 Y5 30 0 0.0 0.0
15+
X131 Y5 40 0 0.0 0.0
16+
X141 Y5 50 0 0.0 0.0
17+
X151 Y5 60 0 0.0 0.0
18+
X161 Y5 70 0 0.0 0.0
19+
X171 Y5 80 0 0.0 0.0
20+
X181 Y5 90 0 0.0 0.0
21+
LIGO HN 46.455144 -119.407656 0.0 0.0
22+
LIGO LV 30.562894 -90.774242 0.0 0.0
23+
PRNC GY 40.345763 -74.654775 0.0 0.0

‎STATIONS2

+3
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,3 @@
1+
2
2+
TLY Y5 51.6807 103.6438 0.0 0.0
3+
YSS Y5 46.9583 142.7610 0.0 0.0

‎STATIONSALL

+20
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,20 @@
1+
19
2+
X1 Y5 -90 0 0.0 0.0
3+
X11 Y5 -80 0 0.0 0.0
4+
X21 Y5 -70 0 0.0 0.0
5+
X31 Y5 -60 0 0.0 0.0
6+
X41 Y5 -50 0 0.0 0.0
7+
X51 Y5 -40 0 0.0 0.0
8+
X61 Y5 -30 0 0.0 0.0
9+
X71 Y5 -20 0 0.0 0.0
10+
X81 Y5 -10 0 0.0 0.0
11+
X91 Y5 0 0 0.0 0.0
12+
X101 Y5 10 0 0.0 0.0
13+
X111 Y5 20 0 0.0 0.0
14+
X121 Y5 30 0 0.0 0.0
15+
X131 Y5 40 0 0.0 0.0
16+
X141 Y5 50 0 0.0 0.0
17+
X151 Y5 60 0 0.0 0.0
18+
X161 Y5 70 0 0.0 0.0
19+
X171 Y5 80 0 0.0 0.0
20+
X181 Y5 90 0 0.0 0.0

‎STATIONSOLD

+20
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,20 @@
1+
19
2+
X1 Y5 -90 0 0.0 0.0
3+
X11 Y5 -80 0 0.0 0.0
4+
X21 Y5 -70 0 0.0 0.0
5+
X31 Y5 -60 0 0.0 0.0
6+
X41 Y5 -50 0 0.0 0.0
7+
X51 Y5 -40 0 0.0 0.0
8+
X61 Y5 -30 0 0.0 0.0
9+
X71 Y5 -20 0 0.0 0.0
10+
X81 Y5 -10 0 0.0 0.0
11+
X91 Y5 0 0 0.0 0.0
12+
X101 Y5 10 0 0.0 0.0
13+
X111 Y5 20 0 0.0 0.0
14+
X121 Y5 30 0 0.0 0.0
15+
X131 Y5 40 0 0.0 0.0
16+
X141 Y5 50 0 0.0 0.0
17+
X151 Y5 60 0 0.0 0.0
18+
X161 Y5 70 0 0.0 0.0
19+
X171 Y5 80 0 0.0 0.0
20+
X181 Y5 90 0 0.0 0.0

‎STATIONS_18

+20
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,20 @@
1+
19
2+
X1 Y5 -90 0 0.0 0.0
3+
X11 Y5 -80 0 0.0 0.0
4+
X21 Y5 -70 0 0.0 0.0
5+
X31 Y5 -60 0 0.0 0.0
6+
X41 Y5 -50 0 0.0 0.0
7+
X51 Y5 -40 0 0.0 0.0
8+
X61 Y5 -30 0 0.0 0.0
9+
X71 Y5 -20 0 0.0 0.0
10+
X81 Y5 -10 0 0.0 0.0
11+
X91 Y5 0 0 0.0 0.0
12+
X101 Y5 10 0 0.0 0.0
13+
X111 Y5 20 0 0.0 0.0
14+
X121 Y5 30 0 0.0 0.0
15+
X131 Y5 40 0 0.0 0.0
16+
X141 Y5 50 0 0.0 0.0
17+
X151 Y5 60 0 0.0 0.0
18+
X161 Y5 70 0 0.0 0.0
19+
X171 Y5 80 0 0.0 0.0
20+
X181 Y5 90 0 0.0 0.0

0 commit comments

Comments
 (0)