[CIG-SEISMO] Mineos DEMO3 test produces different arrival times than in DEMO2
Thorsten Becker
twb at usc.edu
Fri Mar 21 17:38:00 PDT 2008
Hi,
This is probably a simple user error, but I cannot work it out right
now, nor did I find anything on the mailing list. So, thanks for your
comments !
I'm trying to run the Mineos (SVN version) DEMO3 test case and get
results different from those in DEMO2 for some stations. Specifically,
I compile the code on Linux, with intel 10.1, -O -static options, then
run "RUN_MINEOS prem_noocean" in DEMO3, get the first attachment as
output. Then, I convert with "cucss2sac -a", and plot the traces with
GMT script "ptrace" as attached. For some stations, the arrival times
are fine (e.g. BJT), but for others the DEMO3 output is shifted a lot
(e.g. ALE).
>From looking at the scripts, it seems that all input is the same up to
the syndat step, which produces the offset traces.
Thanks much for any help in advance!
Cheers
Thorsten
--
Thorsten Becker Department of Earth Sciences
University of Southern California 213.740.8365
Los Angeles CA 90089-0740 geodynamics.usc.edu/~becker
-------------- next part --------------
Step 1: minos_bran runs for S modes .....................
============== Program minos_bran ==================
input model file:
../models/prem_noocean.txt
output file:
prem_noocean_S
eigenfunction file (output):
eprem_noocean_S
enter eps and wgrav
1.000000000000000E-010 1.00000000000000
enter jcom (1=rad;2=tor;3=sph;4=ictor)
3
enter lmin,lmax,wmin,wmax,nmin,nmax
2 8000 0.000000000000000E+000 200.000000000000
0 0
Step 2: minos_bran runs for T modes .....................
============== Program minos_bran ==================
input model file:
../models/prem_noocean.txt
output file:
prem_noocean_T
eigenfunction file (output):
eprem_noocean_T
enter eps and wgrav
1.000000000000000E-010 1.00000000000000
enter jcom (1=rad;2=tor;3=sph;4=ictor)
2
enter lmin,lmax,wmin,wmax,nmin,nmax
2 8000 0.000000000000000E+000 200.000000000000
0 0
Step 3: eigen for S .....................................
============= Program eigcon ====================
spheroidals (3) or toroidals (2) or radial (1) or
inner core toroidals (4) modes
3
enter name of model file
../models/prem_noocean.txt
enter max depth [km] :
1000.000
enter name of minos_bran output text file
prem_noocean_S
minos_bran output binary unformatted file name
eprem_noocean_S
enter path/dbase_name or dbase_name to store eigenfunctions:
test_S
====================================================
eigcon: n,nstart,nrad = 185 119 67
Step 4: eigen for T .....................................
============= Program eigcon ====================
spheroidals (3) or toroidals (2) or radial (1) or
inner core toroidals (4) modes
2
enter name of model file
../models/prem_noocean.txt
enter max depth [km] :
1000.000
enter name of minos_bran output text file
prem_noocean_T
minos_bran output binary unformatted file name
eprem_noocean_T
enter path/dbase_name or dbase_name to store eigenfunctions:
test_T
====================================================
eigcon: n,nstart,nrad = 185 119 67
Step 5: green functions evaluation .........................
============= Program green ====================
enter path to db with sta & stachan:
short
enter name of file within list of nmodes db:
db_list
enter input CMT file name:
china_cmt_event
min and max frequencies to be considered (mHz) :
10.00000 260.0000
enter # pts in greens fns .le. 30000 :
8000
enter Green functions output db file name:
green
====================================================
green: Event: B011400F 2000 14 23:37:10.800 lat = 25.390, lon = 101.400
green: source depth = 33.0 km
green: step = 1.000 sec, nsamples = 8000
green: # sph. modes in band = 2596 must be .le. 300000
green: # tor. modes in band = 2352 must be .le. 300000
green: Input dbname : short
green: Station: ALE 82.5033 -62.3500 , Channels: 3
green: Channel: # 1 ALE LHZ 0.0000000E+00 0.0000000E+00
green: Channel: # 2 ALE LHN 0.0000000E+00 90.00000
green: Channel: # 3 ALE LHE 90.00000 90.00000
green: Epicentral Distance : 72.020
green: Azimuth of Source : 15.433
green: 1 ALE LHZ 2000 14 23:37:10.800 1.000 48000
green: 2 ALE LHN 2000 14 23:37:10.800 1.000 48000
green: 3 ALE LHE 2000 14 23:37:10.800 1.000 48000
green: Station: ANMO 34.9502 -106.4602 , Channels: 3
green: Channel: # 1 ANMO LHZ 0.0000000E+00 0.0000000E+00
green: Channel: # 2 ANMO LH2 10.00000 90.00000
green: Channel: # 3 ANMO LH1 280.0000 90.00000
green: Epicentral Distance : 114.439
green: Azimuth of Source : -27.665
green: 4 ANMO LHZ 2000 14 23:37:10.800 1.000 48000
green: 5 ANMO LH2 2000 14 23:37:10.800 1.000 48000
green: 6 ANMO LH1 2000 14 23:37:10.800 1.000 48000
green: Station: BAK 35.3441 -119.1043 , Channels: 3
green: Channel: # 1 BAK BHZ 0.0000000E+00 0.0000000E+00
green: Channel: # 2 BAK BHN 0.0000000E+00 90.00000
green: Channel: # 3 BAK BHE 90.00000 90.00000
green: Epicentral Distance : 108.462
green: Azimuth of Source : -38.270
green: 7 BAK BHZ 2000 14 23:37:10.800 1.000 48000
green: 8 BAK BHN 2000 14 23:37:10.800 1.000 48000
green: 9 BAK BHE 2000 14 23:37:10.800 1.000 48000
green: Station: BDFB -15.6418 -48.0148 , Channels: 3
green: Channel: # 1 BDFB LHZ 0.0000000E+00 0.0000000E+00
green: Channel: # 2 BDFB LHN 0.0000000E+00 90.00000
green: Channel: # 3 BDFB LHE 90.00000 90.00000
green: Epicentral Distance : 149.822
green: Azimuth of Source : 66.283
green: 10 BDFB LHZ 2000 14 23:37:10.800 1.000 48000
green: 11 BDFB LHN 2000 14 23:37:10.800 1.000 48000
green: 12 BDFB LHE 2000 14 23:37:10.800 1.000 48000
green: Station: BILL 68.0651 166.4524 , Channels: 3
green: Channel: # 1 BILL LHZ 0.0000000E+00 0.0000000E+00
green: Channel: # 2 BILL LHN 0.0000000E+00 90.00000
green: Channel: # 3 BILL LHE 90.00000 90.00000
green: Epicentral Distance : 57.417
green: Azimuth of Source : -103.266
green: 13 BILL LHZ 2000 14 23:37:10.800 1.000 48000
green: 14 BILL LHN 2000 14 23:37:10.800 1.000 48000
green: 15 BILL LHE 2000 14 23:37:10.800 1.000 48000
green: Station: BJT 40.0183 116.1679 , Channels: 3
green: Channel: # 1 BJT LHZ 0.0000000E+00 0.0000000E+00
green: Channel: # 2 BJT LHN 0.0000000E+00 90.00000
green: Channel: # 3 BJT LHE 90.00000 90.00000
green: Epicentral Distance : 19.123
green: Azimuth of Source : -135.267
green: 16 BJT LHZ 2000 14 23:37:10.800 1.000 48000
green: 17 BJT LHN 2000 14 23:37:10.800 1.000 48000
green: 18 BJT LHE 2000 14 23:37:10.800 1.000 48000
green: Station: BRVK 53.0581 70.2828 , Channels: 3
green: Channel: # 1 BRVK LHZ 0.0000000E+00 0.0000000E+00
green: Channel: # 2 BRVK LHN 0.0000000E+00 90.00000
green: Channel: # 3 BRVK LHE 90.00000 90.00000
green: Epicentral Distance : 36.158
green: Azimuth of Source : 127.603
green: 19 BRVK LHZ 2000 14 23:37:10.800 1.000 48000
green: 20 BRVK LHN 2000 14 23:37:10.800 1.000 48000
green: 21 BRVK LHE 2000 14 23:37:10.800 1.000 48000
green: Station: CASY -66.2792 110.5364 , Channels: 3
green: Channel: # 1 CASY LHZ 0.0000000E+00 0.0000000E+00
green: Channel: # 2 CASY LHN 180.0000 90.00000
green: Channel: # 3 CASY LHE 270.0000 90.00000
green: Epicentral Distance : 91.644
green: Azimuth of Source : -8.261
green: 22 CASY LHZ 2000 14 23:37:10.800 1.000 48000
green: 23 CASY LHN 2000 14 23:37:10.800 1.000 48000
green: 24 CASY LHE 2000 14 23:37:10.800 1.000 48000
green: Station: CCM 38.0557 -91.2446 , Channels: 3
green: Channel: # 1 CCM LHZ 0.0000000E+00 0.0000000E+00
green: Channel: # 2 CCM LHN 0.0000000E+00 90.00000
green: Channel: # 3 CCM LHE 90.00000 90.00000
green: Epicentral Distance : 115.783
green: Azimuth of Source : -12.703
green: 25 CCM LHZ 2000 14 23:37:10.800 1.000 48000
green: 26 CCM LHN 2000 14 23:37:10.800 1.000 48000
green: 27 CCM LHE 2000 14 23:37:10.800 1.000 48000
green: Station: TLY 51.6807 103.6438 , Channels: 3
green: Channel: # 1 TLY LHZ 0.0000000E+00 0.0000000E+00
green: Channel: # 2 TLY LHN 0.0000000E+00 90.00000
green: Channel: # 3 TLY LHE 90.00000 90.00000
green: Epicentral Distance : 26.308
green: Azimuth of Source : -175.417
green: 28 TLY LHZ 2000 14 23:37:10.800 1.000 48000
green: 29 TLY LHN 2000 14 23:37:10.800 1.000 48000
green: 30 TLY LHE 2000 14 23:37:10.800 1.000 48000
Step 6: synthetic seismogram construction ..................
============== Program syndat ==================
enter input CMT file name:
china_cmt_event
enter tensor type: 0 - moment, 1 - nodal plane 1, 2 - nodal plane 2
0
enter input dbname
green
enter output dbname
Syndat
syndat: CMT solution:
syndat: Event: B011400F 2000 14 23: 37: 10.8, lat/lon = 25.39 101.40,
syndat: depth = 33.0 km, step = 1.000 sec, Duration = -2.3 sec,
syndat: M0 = 0.833E+25,
syndat: comp: -0.600E+24 -0.629E+25 0.689E+25 -0.185E+25 0.120E+24 -0.473E+25
syndat: plane1: 27. 78. -6. plane2: 118. 84.-168.
syndat: Selected moment tensor components:
syndat: -0.600E+24 -0.629E+25 0.689E+25 -0.185E+25 0.120E+24 -0.473E+25
syndat: Synthetic waveforms:
syndat: 1 ALE LHZ 1.000 8000
syndat: 2 ALE LHN 1.000 8000
syndat: 3 ALE LHE 1.000 8000
syndat: 4 ANMO LHZ 1.000 8000
syndat: 5 ANMO LH2 1.000 8000
syndat: 6 ANMO LH1 1.000 8000
syndat: 7 BAK BHZ 1.000 8000
syndat: 8 BAK BHN 1.000 8000
syndat: 9 BAK BHE 1.000 8000
syndat: 10 BDFB LHZ 1.000 8000
syndat: 11 BDFB LHN 1.000 8000
syndat: 12 BDFB LHE 1.000 8000
syndat: 13 BILL LHZ 1.000 8000
syndat: 14 BILL LHN 1.000 8000
syndat: 15 BILL LHE 1.000 8000
syndat: 16 BJT LHZ 1.000 8000
syndat: 17 BJT LHN 1.000 8000
syndat: 18 BJT LHE 1.000 8000
syndat: 19 BRVK LHZ 1.000 8000
syndat: 20 BRVK LHN 1.000 8000
syndat: 21 BRVK LHE 1.000 8000
syndat: 22 CASY LHZ 1.000 8000
syndat: 23 CASY LHN 1.000 8000
syndat: 24 CASY LHE 1.000 8000
syndat: 25 CCM LHZ 1.000 8000
syndat: 26 CCM LHN 1.000 8000
syndat: 27 CCM LHE 1.000 8000
syndat: 28 TLY LHZ 1.000 8000
syndat: 29 TLY LHN 1.000 8000
syndat: 30 TLY LHE 1.000 8000
-------------- next part --------------
A non-text attachment was scrubbed...
Name: ptrace
Type: application/x-shellscript
Size: 1021 bytes
Desc: not available
Url : http://geodynamics.org/pipermail/cig-seismo/attachments/20080321/496ec161/attachment-0001.bin
-------------- next part --------------
A non-text attachment was scrubbed...
Name: BJT.pdf
Type: application/pdf
Size: 64096 bytes
Desc: not available
Url : http://geodynamics.org/pipermail/cig-seismo/attachments/20080321/496ec161/attachment-0002.pdf
-------------- next part --------------
A non-text attachment was scrubbed...
Name: ALE.pdf
Type: application/pdf
Size: 61845 bytes
Desc: not available
Url : http://geodynamics.org/pipermail/cig-seismo/attachments/20080321/496ec161/attachment-0003.pdf
More information about the CIG-SEISMO
mailing list