#
# History file 'sc-ss-an.h' for computation of anisotropic ray # theory S-wave synthetic seismograms in models SC1_I, SC1_II and ORT # using the program package ANRAY. # The seismograms are compared with results calculated by Fourier # Method. # Input files required # ~~~~~~~~~~~~~~~~~~~~ chk.pl: "anray/data/qi" "anraysc1.qi" chk.pl: "anray/data/qi" "anraysc2.qi" chk.pl: "anray/data/qi" "anrayort.qi" chk.pl: "anray/data/qi" "fresan.scv" chk.pl: "anray/data/qi" "fresan.sct" chk.pl: "anray/data/qi" "fresan.scr" chk.pl: "anray/data/qi" "synfan.sc" # GSE files for FM - results of "data/qi/fm/fm2gse.h" chk.pl: "" "sc1-x-fm.gse" chk.pl: "" "sc1-y-fm.gse" chk.pl: "" "sc1-z-fm.gse" chk.pl: "" "sc2-x-fm.gse" chk.pl: "" "sc2-y-fm.gse" chk.pl: "" "sc2-z-fm.gse" chk.pl: "" "ort-x-fm.gse" chk.pl: "" "ort-y-fm.gse" chk.pl: "" "ort-z-fm.gse" # Data to control seismogram plotting (program SP) KODESP=3 SPCHRH=0.25 SPTMIN=0.35 SPTMAX=0.65 SPTLEN=6.0 SPTDIV=3 SPTSUB=10 SPXMIN=0.00 SPXMAX=1.40 SPXLEN=15.0 SPXDIV=7 SPXSUB=5 SPTTEXT='Time (s)' SPXTEXT='Depth (km)' SPTDEC=2 NORMSP=1 SPAMP=1. CALCOPS='0.50 setlinewidth' GSEWIDTH=8 (writing GSE files in greenss.for) # Data describing the source time function # KSIG=1 # (Gabor signal) # SIGT=0. SIGF=50. SIGW=4. SIGPH=0. SIGA=100. # Data describing the filtration of the source time function # FMIN=0.0 FLOW=5.0 FHIGH=60. FMAX=100. # (cosine filter) # Time and frequency step and interval for the Fast Fourier Transform # DT=0.000488 NFFT=2048 # OF=0.0 NF=101 # DF=1.000576332 # DT=0.000250 NFFT=2048 # OF= NF= DF= # ===================================================================== # Model SC1_I # ===================================================================== # Anisotropic ray theory synthetic seismograms - package ANRAY # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ REC= SPXMAX=1.60 SPXDIV=8 SPTMAX=0.70 SPTDIV=7 SPTSUB=5 anray: 'anraysc1.qi' 1* 'lu1.out' 'lu2.out' / # Vertical component fresan: 'fresan.scv' 1* 'lu2.out' 'lu7.out' ' ' / synfan: 'synfan.sc' 1* 'lu7.out' 'lu8.out' / LU8='lu8.out' SS='sc1-ss3a.gse' anraygse: SS='sc1-z-fm.gse' SPAMP=100000000000000. KOLOR=1 SS1='sc1-ss3a.gse' SPAMP1=1. KOLOR1=4 SP1=' ' SP2=' ' SP3='sc1-ss3a.ps' sp: # Transverse component fresan: 'fresan.sct' 1* 'lu2.out' 'lu7.out' ' ' / synfan: 'synfan.sc' 1* 'lu7.out' 'lu8.out' / SS='sc1-ss2a.gse' anraygse: SS='sc1-y-fm.gse' SPAMP=100000000000000. KOLOR=1 SS1='sc1-ss2a.gse' SPAMP1=1. KOLOR1=4 SP1=' ' SP2='sc1-ss2a.ps' SP3=' ' sp: # Radial component fresan: 'fresan.scr' 1* 'lu2.out' 'lu7.out' ' ' / synfan: 'synfan.sc' 1* 'lu7.out' 'lu8.out' / SS='sc1-ss1a.gse' anraygse: SS='sc1-x-fm.gse' SPAMP=100000000000000. KOLOR=1 SS1='sc1-ss1a.gse' SPAMP1=1. KOLOR1=4 SP1='sc1-ss1a.ps' SP2=' ' SP3=' ' sp: SPAMP=1. SPAMP1= KOLOR= KOLOR1= SS= SS1= SPXMAX=1.40 SPXDIV=7 SPTMAX=0.65 SPTDIV=3 SPTSUB=10 # ===================================================================== # Model SC1_II # ===================================================================== # Anisotropic ray theory synthetic seismograms - package ANRAY # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ REC= SPXMAX=1.60 SPXDIV=8 SPTMAX=0.70 SPTDIV=7 SPTSUB=5 anray: 'anraysc2.qi' 1* 'lu1.out' 'lu2.out' / # Vertical component fresan: 'fresan.scv' 1* 'lu2.out' 'lu7.out' ' ' / synfan: 'synfan.sc' 1* 'lu7.out' 'lu8.out' / LU8='lu8.out' SS='sc2-ss3a.gse' anraygse: SS='sc2-z-fm.gse' SPAMP=100000000000000. KOLOR=1 SS1='sc2-ss3a.gse' SPAMP1=1. KOLOR1=4 SP1=' ' SP2=' ' SP3='sc2-ss3a.ps' sp: # Transverse component fresan: 'fresan.sct' 1* 'lu2.out' 'lu7.out' ' ' / synfan: 'synfan.sc' 1* 'lu7.out' 'lu8.out' / SS='sc2-ss2a.gse' anraygse: SS='sc2-y-fm.gse' SPAMP=100000000000000. KOLOR=1 SS1='sc2-ss2a.gse' SPAMP1=1. KOLOR1=4 SP1=' ' SP2='sc2-ss2a.ps' SP3=' ' sp: # Radial component fresan: 'fresan.scr' 1* 'lu2.out' 'lu7.out' ' ' / synfan: 'synfan.sc' 1* 'lu7.out' 'lu8.out' / SS='sc2-ss1a.gse' anraygse: SS='sc2-x-fm.gse' SPAMP=100000000000000. KOLOR=1 SS1='sc2-ss1a.gse' SPAMP1=1. KOLOR1=4 SP1='sc2-ss1a.ps' SP2=' ' SP3=' ' sp: SPAMP=1. SPAMP1= KOLOR= KOLOR1= SS= SS1= SPXMAX=1.40 SPXDIV=7 SPTMAX=0.65 SPTDIV=3 SPTSUB=10 # ===================================================================== # Model ORT # ===================================================================== # Anisotropic ray theory synthetic seismograms - package ANRAY # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ REC= SPTMIN=0.65 SPTMAX=1.15 SPTLEN=6.0 SPTDIV=5 SPTSUB=10 SPXMIN=0.00 SPXMAX=1.60 SPXLEN=15.0 SPXDIV=8 SPXSUB=5 anray: 'anrayort.qi' 1* 'lu1.out' 'lu2.out' / # Vertical component fresan: 'fresan.scv' 1* 'lu2.out' 'lu7.out' ' ' / synfan: 'synfan.sc' 1* 'lu7.out' 'lu8.out' / LU8='lu8.out' SS='ort-ss3a.gse' anraygse: SS='ort-z-fm.gse' SPAMP=15000000000000. KOLOR=1 SS1='ort-ss3a.gse' SPAMP1=0.15 KOLOR1=4 SP1=' ' SP2=' ' SP3='ort-ss3a.ps' sp: # Transverse component fresan: 'fresan.sct' 1* 'lu2.out' 'lu7.out' ' ' / synfan: 'synfan.sc' 1* 'lu7.out' 'lu8.out' / SS='ort-ss2a.gse' anraygse: SS='ort-y-fm.gse' SPAMP=15000000000000. KOLOR=1 SS1='ort-ss2a.gse' SPAMP1=0.15 KOLOR1=4 SP1=' ' SP2='ort-ss2a.ps' SP3=' ' sp: # Radial component fresan: 'fresan.scr' 1* 'lu2.out' 'lu7.out' ' ' / synfan: 'synfan.sc' 1* 'lu7.out' 'lu8.out' / SS='ort-ss1a.gse' anraygse: SS='ort-x-fm.gse' SPAMP=15000000000000. KOLOR=1 SS1='ort-ss1a.gse' SPAMP1=0.15 KOLOR1=4 SP1='ort-ss1a.ps' SP2=' ' SP3=' ' sp: SPAMP=1. SPAMP1= KOLOR= KOLOR1= SS= SS1= SPTMIN=0.35 SPTMAX=0.65 SPTLEN=6.0 SPTDIV=3 SPTSUB=10 SPXMIN=0.00 SPXMAX=1.40 SPXLEN=15.0 SPXDIV=7 SPXSUB=5 # List of the output files # ~~~~~~~~~~~~~~~~~~~~~~~~ # '???-ss1a.ps','???-ss2a.ps','???-ss3a.ps' ... Anisotropic ray theory # synthetic seismograms compared with FM. # "???" in the names of files identifies the model # sc1 ... model SC1_I # sc2 ... model SC1_II # ort ... model ORT # Colors of seismograms '*ss?a.ps' are # Black - Fourier Method # Blue - Anisotropic ray theory # Figures correspond to Psencik, Farra & Tessmer (2011): # Comparison of the FORT approximation of the coupling ray theory with # the Fourier pseudospectral method # In: Seismic Waves in Complex 3-D Structures, Report 21, # Dep. Geophys., Charles Univ., Prague, pp. 129-165. #