# History file 'p1-ss.h' generates the Born seismogram from the RF # format file computed by history file 'p1-born.h' and visualizes it # together with the ray-theory seismogram computed by history file # 'p1-ss2.h'. # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ # Input files required chk.pl: "data/p1i/" "p1-recb.dat" chk.pl: "data/p1i/" "p1-pts.dat" chk.pl: "data/p1i/" "p1-spttc.dat" chk.pl: "data/p1i/" "p1-edge2.dat" chk.pl: "data/p1i/" "p1-edge4.dat" chk.pl: "data/p1i/" "p1-edge5.dat" chk.pl: "data/p1i/" "p1-edge6.dat" # output of 'p1-born.h' chk.pl: "" "rf-born.out" # output of 'p1-ss2.h' chk.pl: "" "p1-rss.gse" # Data for ss.for # Data describing the source time function KSIG=1 # (Gabor signal) SIGT=0. SIGF=10. SIGW=4. SIGPH=0. SIGA=258.174 # (258.174=12*pi*rho*v**3) # Data describing the filtration of the source time function FMIN=1. FLOW=2. FHIGH=19. FMAX=20. # (cosine filter) # Time step and time interval for the Fast Fourier Transform DT=0.001 NFFT=16384 # Data to control optional plotting of program SS SIGPLOT='p1-sig.ps' MPTS=100 # Name of the input RF format file RF='rf-born.out' # Data for sp.for KODESP= SPCHRH=0.40 SPTMIN=1.5 SPTMAX=5.4 SPXLEN=17 SPTDIV=3.9 SPTSUB=10 SPTLEN=25 NORMSP=1 SPAMP=5.5e6 SPTEXT3='BLACK...Born approximation' SPTEXT4='RED...Ray theory' PTS='p1-pts.dat' SPTTC='p1-spttc.dat' REC='p1-recb.dat' # Frequency-domain seismograms SS='p1-ss.gse' SS1='p1-rss.gse' SP1='p1-ss1.ps' SP2='p1-ss2.ps' SP3='p1-ss3.ps' KOLOR=1 KOLOR1=2 SPTDEC=1 ss: # synthetic seismograms sp: # seismogram plotting