#
# History file 'qi4-ss.h' for computation of anisotropic, isotropic, and # coupling-ray theory S-wave synthetic seismograms in model QI4. # Executables required # ~~~~~~~~~~~~~~~~~~~~ # This history file requires programs FRESAN, SYNFAN and ANRAYGSE # of package ANRAY, and modified version of program ANRAY, please link # program 'anray_44.for' located in directory 'data/qi'. # Input files required # ~~~~~~~~~~~~~~~~~~~~ chk.pl: "data/qi/" "anrayqi4.qi" chk.pl: "data/qi/" "fresan.anv" chk.pl: "data/qi/" "fresan.ant" chk.pl: "data/qi/" "fresan.anr" chk.pl: "data/qi/" "synfan.an" chk.pl: "data/qi/" "qi-mod4.dat" chk.pl: "data/qi/" "qi-src.dat" chk.pl: "data/qi/" "qi-rec0.dat" chk.pl: "data/qi/" "qi-crt.dat" chk.pl: "data/qi/" "qi-rpa.dat" chk.pl: "data/qi/" "qi-rpaa.dat" chk.pl: "data/qi/" "force.dat" chk.pl: "crt/" "writ.dat" # edited Green function combined from S1 and S2 (or SH and SV) rays chk.pl: "data/qi/" "qi4-gr12.dat" chk.pl: "data/qi/" "qi4-grhv.dat" # GSE files for FORT - results of "data/qi/fort/fort2gse.h" chk.pl: "" "qi4-x-fo.gse" chk.pl: "" "qi4-y-fo.gse" chk.pl: "" "qi4-z-fo.gse" # GSE files for FM - results of "data/qi/fm/fm2gse.h" chk.pl: "" "qi4-x-fm.gse" chk.pl: "" "qi4-y-fm.gse" chk.pl: "" "qi4-z-fm.gse" # Data to control seismogram plotting (program SP) KODESP=3 SPCHRH=0.25 SPTMIN=0.40 SPTMAX=0.55 SPTLEN=6.0 SPTDIV=3 SPTSUB=5 SPXMIN=0.0 SPXMAX=0.60 SPXLEN=15.0 SPXDIV=6 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) # Anisotropic ray theory synthetic seismograms # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ anray_44: 'anrayqi4.qi' 1* 'lu1.out' 'lu2.out' / # Vertical component fresan: 'fresan.anv' 1* 'lu2.out' 'lu7.out' ' ' / synfan: 'synfan.an' 1* 'lu7.out' 'lu8.out' / LU8='lu8.out' SS='qi4-ss3a.gse' anraygse: SP1=' ' SP2=' ' SP3='qi4-ss3a.ps' sp: # Transverse component fresan: 'fresan.ant' 1* 'lu2.out' 'lu7.out' ' ' / synfan: 'synfan.an' 1* 'lu7.out' 'lu8.out' / SS='qi4-ss2a.gse' anraygse: SP1=' ' SP2='qi4-ss2a.ps' SP3=' ' sp: # Radial component fresan: 'fresan.anr' 1* 'lu2.out' 'lu7.out' ' ' / synfan: 'synfan.an' 1* 'lu7.out' 'lu8.out' / SS='qi4-ss1a.gse' anraygse: SP1='qi4-ss1a.ps' SP2=' ' SP3=' ' sp: # Isotropic ray theory synthetic seismograms # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ SRC='qi-src.dat' REC='qi-rec0.dat' SOURCE='force.dat' DCRT='qi-crt.dat' CODE='qi-crt.dat' WRIT='writ.dat' INIPAR=3 # 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.000250 NFFT=2048 # OF=0.0 NF=101 # DF=1.000576332 OF= NF= DF= # Ray tracing and synthetic seismograms MODEL='qi-mod4.dat' CRTANI=0 RPAR='qi-rpa.dat' CRTLOG='qi-crt.out' crt: MODEL=' ' GREEN='green.out' SS='qi4-ssi.gse' green: greenss: ss: SP1='qi4-ss1i.ps' SP2='qi4-ss2i.ps' SP3='qi4-ss3i.ps' sp: # Coupling ray theory synthetic seismograms - anisotropic common ray # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ MODEL='qi-mod4.dat' CRTANI=1 RPAR='qi-rpaa.dat' CRTLOG='qi-crta.out' SS='qi4-ssc.gse' crt: green: greenss: ss: SP1='qi4-ss1c.ps' SP2='qi4-ss2c.ps' SP3='qi4-ss3c.ps' sp: # Coupling ray theory synthetic seismograms - anisotropic common ray, second-order perturbations of travel time # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ MODEL='qi-mod4.dat' CRTANI=1 RPAR='qi-rpaa.dat' CRTLOG='qi-crta.out' QIRAY=1 SS='qi4-sscs.gse' crt: green: greenss: ss: SP1='qi4-ss1s.ps' SP2='qi4-ss2s.ps' SP3='qi4-ss3s.ps' sp: # Prevailing-frequency approximation # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ CRTLOG='qi-crtaf.out' SS='qi4-sscf.gse' crt: SINGLF=50. green: SINGLF= greenss: ss: SP1='qi4-ss1f.ps' SP2='qi4-ss2f.ps' SP3='qi4-ss3f.ps' sp: KOLOR=2 (red) KOLOR1=1 (black) SPAMP1=100000000000000. SP1='qi4-ca-1.ps' SS1='qi4-x-fm.gse' SP2=' ' SP3=' ' sp: SP2='qi4-ca-2.ps' SS1='qi4-y-fm.gse' SP1=' ' SP3=' ' sp: SP3='qi4-ca-3.ps' SS1='qi4-z-fm.gse' SP1=' ' SP2=' ' sp: SPAMP1= KOLOR= KOLOR1= SS= SS1= # Coupling ray theory synthetic seismograms, # two S1 and S2 anisotropic rays, # prevailing-frequency approximation # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ MODEL='qi-mod4.dat' CRTANI=1 RPAR='qi-rpaa.dat' QIRAY= KSWAVE=1 SS= CRTLOG='qi-crts1.out' crt: SINGLF=50. GREEN='qi4-gr01.out' green: SINGLF= KSWAVE=2 SS= CRTLOG='qi-crts2.out' crt: SINGLF=50. GREEN='qi4-gr02.out' green: SINGLF= # Green functions 'qi4-gr01.out' and 'qi4-gr02.out' should be manually # edited into file 'qi4-gr12.dat' SS='qi4-sf12.gse' GREEN='qi4-gr12.dat' greenss: ss: GREEN= KSWAVE= KOLOR=2 (red) KOLOR1=1 (black) SPAMP1=100000000000000. SP1='qi4-12-1.ps' SS1='qi4-x-fm.gse' SP2=' ' SP3=' ' sp: SP2='qi4-12-2.ps' SS1='qi4-y-fm.gse' SP1=' ' SP3=' ' sp: SP3='qi4-12-3.ps' SS1='qi4-z-fm.gse' SP1=' ' SP2=' ' sp: SPAMP1= KOLOR= KOLOR1= SS= SS1= # Coupling ray theory synthetic seismograms, # two SH and SV anisotropic rays, # prevailing-frequency approximation # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ MODEL='qi-mod4.dat' CRTANI=1 RPAR='qi-rpaa.dat' QIRAY= TIA1=0.707107 TIA2=0.707107 TIA3=0. DEGREE=2 KSWAVE=1 SS= CRTLOG='qi-crtsh.out' crt: SINGLF=50. GREEN='qi4-grsh.out' green: SINGLF= KSWAVE=2 SS= CRTLOG='qi-crtsv.out' crt: SINGLF=50. GREEN='qi4-grsv.out' green: SINGLF= # Green functions 'qi4-grsh.out' and 'qi4-grsv.out' should be manually # edited into file 'qi4-grhv.dat' SS='qi4-sfhv.gse' GREEN='qi4-grhv.dat' greenss: ss: GREEN= KSWAVE= TIA1= TIA2= TIA3= DEGREE= KOLOR=2 (red) KOLOR1=1 (black) SPAMP1=100000000000000. SP1='qi4-hv-1.ps' SS1='qi4-x-fm.gse' SP2=' ' SP3=' ' sp: SP2='qi4-hv-2.ps' SS1='qi4-y-fm.gse' SP1=' ' SP3=' ' sp: SP3='qi4-hv-3.ps' SS1='qi4-z-fm.gse' SP1=' ' SP2=' ' sp: SPAMP1= KOLOR= KOLOR1= SS= SS1= # Coupling, isotropic and anisotropic seismograms in a single plot # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ # SS='qi4-ssc.gse' SRC= REC= # SS2='qi4-ssi.gse' # SS1='qi4-ss1a.gse' SP1='qi4-ss1.ps' SP2=' ' SP3=' ' # sp: # SS1='qi4-ss2a.gse' SP1=' ' SP2='qi4-ss2.ps' SP3=' ' # sp: # SS1='qi4-ss3a.gse' SP1=' ' SP2=' ' SP3='qi4-ss3.ps' # sp: # FM plus FORT plus both coupling ray theory seismograms # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ KOLOR=3 (green) KOLOR1=2 (red) KOLOR2=4 (blue) KOLOR3=1 (black) SS='qi4-ssc.gse' SRC= REC= SS1='qi4-x-fo.gse' SS2='qi4-sscs.gse' SS3='qi4-x-fm.gse' SPAMP3=100000000000000. SP1='qi4-ss1.ps' SP2=' ' SP3=' ' sp: SS1='qi4-y-fo.gse' SS3='qi4-y-fm.gse' SP1=' ' SP2='qi4-ss2.ps' SP3=' ' sp: SS1='qi4-z-fo.gse' SS3='qi4-z-fm.gse' SP1=' ' SP2=' ' SP3='qi4-ss3.ps' sp: KOLOR= KOLOR1= KOLOR2= KOLOR3= SS= SS1= SS2= SS3= SPAMP= SPAMP1= SPAMP2= SPAMP3= # Subtracting the seismograms # ~~~~~~~~~~~~~~~~~~~~~~~~~~~ DT=0. SS1='qi4-x-fm.gse' COEF1=-100000000000000. SS2='qi4-x-fo.gse' COEF2=1. SS='qi4-fo1.gse' gsecal: SS1='qi4-y-fm.gse' SS2='qi4-y-fo.gse' SS='qi4-fo2.gse' gsecal: SS1='qi4-z-fm.gse' SS2='qi4-z-fo.gse' SS='qi4-fo3.gse' gsecal: SS1='qi4-x-fm.gse' SS2='qi4-ssc.gse' SS='qi4-cr1.gse' gsecal: SS1='qi4-y-fm.gse' SS2='qi4-ssc.gse' SS='qi4-cr2.gse' gsecal: SS1='qi4-z-fm.gse' SS2='qi4-ssc.gse' SS='qi4-cr3.gse' gsecal: SS1='qi4-x-fm.gse' SS2='qi4-sscs.gse' SS='qi4-crs1.gse' gsecal: SS1='qi4-y-fm.gse' SS2='qi4-sscs.gse' SS='qi4-crs2.gse' gsecal: SS1='qi4-z-fm.gse' SS2='qi4-sscs.gse' SS='qi4-crs3.gse' gsecal: SPAMP=2. KOLOR=3 (green) KOLOR1=2 (red) KOLOR2=4 (blue) KOLOR3= SS='qi4-cr1.gse' SS1='qi4-fo1.gse' SS2='qi4-crs1.gse' SP1='qi4-sd1.ps' SP2=' ' SP3=' ' sp: SS='qi4-cr2.gse' SS1='qi4-fo2.gse' SS2='qi4-crs2.gse' SP1=' ' SP2='qi4-sd2.ps' SP3=' ' sp: SS='qi4-cr3.gse' SS1='qi4-fo3.gse' SS2='qi4-crs3.gse' SP1=' ' SP2=' ' SP3='qi4-sd3.ps' sp: KOLOR= KOLOR1= KOLOR2= KOLOR3= SS= SS1= SS2= SS3= SPAMP= SPAMP1= SPAMP2= SPAMP3= # FM plus coupling ray theory plus prevailing-frequency approximation # seismograms # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ KOLOR=2 red KOLOR1=3 green KOLOR2=1 black SS='qi4-sscf.gse' # prevailing-frequency approximation SS1='qi4-sscs.gse' # coupling ray theory SS2='qi4-x-fm.gse' SPAMP2=100000000000000. # FM SP1='qi4-sf1.ps' SP2=' ' SP3=' ' sp: SS2='qi4-y-fm.gse' SP1=' ' SP2='qi4-sf2.ps' SP3=' ' sp: SS2='qi4-z-fm.gse' SP1=' ' SP2=' ' SP3='qi4-sf3.ps' sp: # List of the output files # ~~~~~~~~~~~~~~~~~~~~~~~~ # 'qi4-ss1?.ps','qi4-ss2?.ps','qi4-ss3?.ps' ... Synthetic seismograms # where ? identifies the calculation # a ... anisotropic ray theory # i ... isotropic ray theory # c ... coupling ray theory along anisotropic common ray # s ... coupling ray theory along anisotropic common ray # with second-order travel time perturbation # f ... prevailing-frequency approximtion of the coupling # ray theory # 'qi4-ca-1.ps','qi4-ca-2.ps','qi4-ca-3.ps' ... Prevailing-frequency # approximation seismograms calculated along the common anisotropic # rays with second-order travel time perturbation (red) compared # with FM (black). # Figures used in Klimes & Bulant (2014) # Prevailing-frequency approximation of the coupling ray theory # for S waves along the SH and SV reference rays in a transversely # isotropic medium, Seismic Waves in Complex 3-D Structures, # 24, 165-177, fig. 4. # 'qi4-12-1.ps','qi4-12-2.ps','qi4-12-3.ps' ... Prevailing-frequency # approximation seismograms calculated along S1 and S2 # anisotropic-ray-theory rays (red) compared with FM (black). # Figures used in Klimes & Bulant (2014) # Prevailing-frequency approximation of the coupling ray theory # for S waves along the SH and SV reference rays in a transversely # isotropic medium, Seismic Waves in Complex 3-D Structures, # 24, 165-177, fig. 5. # 'qi4-hv-1.ps','qi4-hv-2.ps','qi4-hv-3.ps' ... Prevailing-frequency # approximation seismograms calculated along SH and SV # anisotropic-ray-theory rays (red) compared with FM (black). # Figures used in Klimes & Bulant (2016), SGG. # 'qi4-sf1.ps','qi4-sf2.ps','qi4-sf3.ps' ... Prevailing-frequency # approximation seismograms compared with standard coupling ray # theory and with FM. # Colors of seismograms are # Black - Fourier Method # Red - Prevailing-frequency approximation of the coupling ray # theory along anisotropic common rays # Green - CRT with second-order travel time perturbation # 'qi4-ss1.ps','qi4-ss2.ps','qi4-ss3.ps' ... Coupling ray theory # seismograms compared with FORT and with FM. # 'qi4-sd1.ps','qi4-sd2.ps','qi4-sd3.ps' ... Differences of the coupling # ray theory and FORT seismograms from FM seismograms. # Colors of seismograms are # Black - Fourier Method # Red - First Order Ray Tracing # Green - Coupling Ray Theory along anisotropic common rays # Blue - CRT with second-order travel time perturbation #