# Coupling ray theory seismograms in weakly anisotropic model QI # Second-order perturbations of travel times # Anisotropic common ray tracing # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ # This history file may be run directly for one of the models # 'qi-mod.dat', 'qi-modh.dat', 'qi-mod2.dat', 'qi-mod4.dat', # 'qi-mod8.dat', or may be called from the history file 'qi-ss4m.h' # for the model 'qi-mod.tmp' - edit the line containing "MODEL=" below. # Executables required # ~~~~~~~~~~~~~~~~~~~~ # This history file executes compiled programs 'greentc.for' and # 'tccomp.for'. Programs 'greentc.for' and 'tccomp.for' are located # in directory 'data/tc'. chk.pl: "data/tc/" "greentc.for" chk.pl: "data/tc/" "tccomp.for" chk.pl: "forms/" "f.pl" f.pl: "greentc" f.pl: "tccomp" # Input files required chk.pl: "data/qi/" "qi-src.dat" chk.pl: "data/qi/" "qi-rec.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" # Model QI, QIH, QI2, QI4 or QI8 # MODEL='qi-mod.dat' chk.pl: "data/qi/" "qi-mod.dat" # MODEL='qi-modh.dat' chk.pl: "data/qi/" "qi-modh.dat" # MODEL='qi-mod2.dat' chk.pl: "data/qi/" "qi-mod2.dat" # MODEL='qi-mod4.dat' chk.pl: "data/qi/" "qi-mod4.dat" # MODEL='qi-mod8.dat' chk.pl: "data/qi/" "qi-mod8.dat" MODEL='qi-mod.tmp' chk.pl: " " "qi-mod.tmp" # Source and receivers SRC='qi-src.dat' REC='qi-rec.dat' SOURCE='force.dat' # Data filenames for CRT (for details see 'crt.for') DCRT='qi-crt.dat' CODE='qi-crt.dat' WRIT='writ.dat' # only two-point rays stored in files # Kind of ray parameters (for details see 'init.for') INIPAR=3 # azimuthal equidistant projection # Parameter to get the same results of ray tracing as in version 5.70 UEBMUL=1. # Degree of the homogeneous Hamiltonian for anisotropic ray tracing DEGREE=-1 # 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 # Data to control optional plotting of program SS MPTS=200 # 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=-.03 SPXMAX=0.60 SPXLEN=15.0 SPXDIV=7 SPXSUB=9 NORMSP=1 SPAMP=0.80 SS1='qi-ss.gse' KOLOR=2 KOLOR1=1 # Isotropic two-point ray tracing # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ CRTANI=0 RPAR='qi-rpa.dat' CRTLOG='qi-crt.out' crt: # Coupling ray theory seismograms # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ QICHM=0 QITT=0 QIDT=0 QIPV=0 QIRAY=0 QILST='qilst.out' GREEN='green.out' SS='qi-ss.gse' # GSE file with synthetic seismograms SIGPLOT='qi-sig.ps' # Optional plot of program SS green: # Green function greenss: # response function ss: # synthetic seismograms GREENTC='greentc.out' greentc: # Quasi-isotropic ray theory seismograms - Christoffel matrix # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ QICHM=1 QITT=0 QIDT=0 QIPV=0 QIRAY=0 QILST=' ' GREEN='greenc.out' SS='qi-ssc.gse' # GSE file with synthetic seismograms SIGPLOT= # no plotting by program SS green: # Green function greenss: # response function ss: # synthetic seismograms # Plotting seismograms SP1='qi-ss1c.ps' SP2='qi-ss2c.ps' SP3='qi-ss3c.ps' sp: # Comparison of plane-wave propagator matrices GREENTC='greentcq.out' GREENTC2='greentc.out' TCCOMP='qi-errc.out' INORM=2 greentc: tccomp: # Quasi-isotropic ray theory seismograms - travel times # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ QICHM=0 QITT=1 QIDT=1 QIPV=0 QIRAY=0 GREEN='greent.out' SS='qi-sst.gse' # GSE file with synthetic seismograms SIGPLOT= # no plotting by program SS green: # Green function greenss: # response function ss: # synthetic seismograms # Plotting seismograms SP1='qi-ss1t.ps' SP2='qi-ss2t.ps' SP3='qi-ss3t.ps' sp: # Comparison of plane-wave propagator matrices GREENTC='greentcq.out' GREENTC2='greentc.out' TCCOMP='qi-errt.out' INORM=2 greentc: tccomp: # Quasi-isotropic ray theory seismograms - polarization vector # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ QICHM=0 QITT=0 QIDT=0 QIPV=1 QIRAY=0 GREEN='greenp.out' SS='qi-ssp.gse' # GSE file with synthetic seismograms SIGPLOT= # no plotting by program SS green: # Green function greenss: # response function ss: # synthetic seismograms # Plotting seismograms SP1='qi-ss1p.ps' SP2='qi-ss2p.ps' SP3='qi-ss3p.ps' sp: # Comparison of plane-wave propagator matrices GREENTC='greentcq.out' GREENTC2='greentc.out' TCCOMP='qi-errp.out' INORM=2 greentc: tccomp: # Quasi-isotropic ray theory seismograms - anisotropic common ray # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ QICHM=0 QITT=0 QIDT=0 QIPV=0 QIRAY=1 GREEN='greenr.out' SS='qi-ssr.gse' # GSE file with synthetic seismograms SIGPLOT= # no plotting by program SS green: # Green function greenss: # response function ss: # synthetic seismograms # Plotting seismograms SS='qi-ss.gse' SS1='qi-ssr.gse' SP1='qi-ss1r.ps' SP2='qi-ss2r.ps' SP3='qi-ss3r.ps' sp: # Comparison of plane-wave propagator matrices GREENTC='greentcq.out' GREENTC2='greentc.out' TCCOMP='qi-errr.out' INORM=2 greentc: tccomp: # Plotting seismograms of 3 different QI approximations together # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ SS='qi-ssc.gse' SS1='qi-ssp.gse' SS2='qi-sst.gse' SS3='qi-ss.gse' KOLOR=2 KOLOR1=4 KOLOR2=3 KOLOR3=1 SP1='qi-ss1.ps' SP2='qi-ss2.ps' SP3='qi-ss3.ps' sp: # Anisotropic common ray two-point ray tracing # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ CRTANI=1 RPAR='qi-rpaa.dat' CRTLOG='qi-crta.out' crt: # Coupling ray theory seismograms along the common anisotropic rays # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ QICHM=0 QITT=0 QIDT=0 QIPV=0 QIRAY=0 QILST='qilsta.out' GREEN='green.out' SS='qi-ssa.gse' # GSE file with synthetic seismograms SIGPLOT= # no plotting by program SS green: # Green function greenss: # response function ss: # synthetic seismograms # Plotting seismograms SS='qi-ss.gse' SS1='qi-ssa.gse' SS2= SS3= KOLOR=2 KOLOR1=1 KOLOR2= KOLOR3= SP1='qi-ss1a.ps' SP2='qi-ss2a.ps' SP3='qi-ss3a.ps' sp: SS='qi-ssr.gse' SS1='qi-ssa.gse' SP1='qi-ss1q.ps' SP2='qi-ss2q.ps' SP3='qi-ss3q.ps' sp: # Comparison of plane-wave propagator matrices GREENTC='greentcq.out' GREENTC2='greentc.out' TCCOMP='qi-erra.out' INORM=2 greentc: tccomp: # Coupling ray theory seismograms along the common anisotropic rays # with the second-order perturbations # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ QICHM=0 QITT=0 QIDT=0 QIPV=0 QIRAY=1 QILST=' ' GREEN='green.out' SS='qi-sss.gse' # GSE file with synthetic seismograms SIGPLOT= # no plotting by program SS green: # Green function greenss: # response function ss: # synthetic seismograms # Plotting seismograms SS='qi-ss.gse' SS1='qi-sss.gse' SS2= SS3= KOLOR=2 KOLOR1=1 KOLOR2= KOLOR3= SP1='qi-ss1v.ps' SP2='qi-ss2v.ps' SP3='qi-ss3v.ps' sp: SS='qi-ssa.gse' SS1='qi-sss.gse' SP1='qi-ss1w.ps' SP2='qi-ss2w.ps' SP3='qi-ss3w.ps' sp: # List of the output files # ~~~~~~~~~~~~~~~~~~~~~~~~ # 'qi-err?.out' ... Norm of the difference from the coupling ray theory # ray propagator matrix, where "?" identifies either the # quasi-isotropic approximation being compared # c ... Christoffel matrix # t ... travel times # p ... polarization vector # r ... anisotropic common ray (perturbation) | second order travel-time perturbation # or the calculation along anisotropic common rays # a ... anisotropic common ray # 'qilst.out' ... Second-order perturbations of travel times. # 'qilsta.out'... Second-order perturbations of travel times calculated # along anisotropic common rays. # 'qi-ss?.gse' ... Synthetic seismograms in GSE format, where "?" # identifies the calculation # ... coupling ray theory along isotropic reference rays # c ... QI approximation of the Christoffel matrix # t ... QI approximation of the travel times # p ... QI approximation of the polarization vector # r ... anisotropic common ray (perturbation) | second order travel-time perturbation # a ... common anisotropic reference rays # s ... common anisotropic reference rays with second order # travel-time perturbation # 'qi-ss1?.ps','qi-ss2?.ps','qi-ss3?.ps' ... PostScript pictures of # synthetic seismograms calculated by different methods. # c ... QI approximation of the Christoffel matrix (Red), # coupling ray theory along isotropic reference rays (Black) # t ... QI approximation of the travel times (Red), # coupling ray theory along isotropic reference rays (Black) # p ... QI approximation of the polarization vector (Red), # coupling ray theory along isotropic reference rays (Black) # r ... anisotropic common ray (perturbation) (Black), | second order perturbation along isotropic rays (Black) # coupling ray theory along isotropic reference rays (Red) # ... QI approximation of the Christoffel matrix (Red), # QI approximation of the travel times (Green), # QI approximation of the polarization vector (Blue), # coupling ray theory along isotropic reference rays (Black) # a ... common anisotropic reference ray (Black), # coupling ray theory along isotropic reference rays (Red) # q ... common anisotropic reference ray (Black), # second order perturbation along isotropic rays (Red) # v ... second order perturbation along common anis. rays (Black), # coupling ray theory along isotropic reference rays (Red) # w ... second order perturbation along common anis. rays (Black), # common anisotropic reference ray (Red)