# 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.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" # 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 NORMSP=1 SPAMP=1. CALCOPS='0.50 setlinewidth' # Anisotropic ray theory synthetic seismograms # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ anray: '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.000488 NFFT=2048 OF=0.0 NF=101 DF=1.000576332 # 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: # 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' sscal: SS1='qi4-y-fm.gse' SS2='qi4-y-fo.gse' SS='qi4-fo2.gse' sscal: SS1='qi4-z-fm.gse' SS2='qi4-z-fo.gse' SS='qi4-fo3.gse' sscal: SS1='qi4-x-fm.gse' SS2='qi4-ssc.gse' SS='qi4-cr1.gse' sscal: SS1='qi4-y-fm.gse' SS2='qi4-ssc.gse' SS='qi4-cr2.gse' sscal: SS1='qi4-z-fm.gse' SS2='qi4-ssc.gse' SS='qi4-cr3.gse' sscal: SS1='qi4-x-fm.gse' SS2='qi4-sscs.gse' SS='qi4-crs1.gse' sscal: SS1='qi4-y-fm.gse' SS2='qi4-sscs.gse' SS='qi4-crs2.gse' sscal: SS1='qi4-z-fm.gse' SS2='qi4-sscs.gse' SS='qi4-crs3.gse' sscal: 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: # 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 # '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