# History file 's145oe5.h' computes two point rays in the elastic # smoothed SEG/EAGE Salt Model, plots computed rays, synthetic # seismograms and travel times. # The shot corresponds to shot number 145 and receivers correspond # to ocean bottom cable line 5 (E-W direction) in the Next Generation # Seismic Modeling and Imaging project, Phase 1 (House et al., 2004). # The time of computation is approximately 12 hours # on a PC equipped with Athlon XP 2400+(2GHz) and with Linux. # # References: # House, S., Larsen, S., Hoelting, C., Marfurt, K., Wiley, R. (2004): # Next-Generation Seismic Modeling and Imaging project: summary of # elastic modeling results. Exp. Abstr., Soc. Expl. Geophys. # 74th Ann. Meet., pp. 2201-2204. # # Increase MRAY=120000 in rp3d.inc! # ====================================================================== # # Input files required chk.pl: "data/sal/" "sal-mod3.dat" chk.pl: "data/sal/" "s145-src.dat" chk.pl: "data/sal/" "s145oe5r.dat" chk.pl: "data/sal/" "oe-dcr1a.dat" chk.pl: "data/sal/" "oe-dcr1b.dat" chk.pl: "data/sal/" "oe-dcr1c.dat" chk.pl: "data/sal/" "oe-cod1a.dat" chk.pl: "data/sal/" "oe-cod1b.dat" chk.pl: "data/sal/" "oe-cod1c.dat" chk.pl: "data/sal/" "oe-rpa1a.dat" chk.pl: "data/sal/" "oe-rpa1b.dat" chk.pl: "data/sal/" "oe-rpa1c.dat" chk.pl: "data/sal/" "oe-wri1.dat" chk.pl: "data/sal/" "oe-ttc.dat" chk.pl: "crt/" "sourcex.dat" chk.pl: "forms/" "append.pl" chk.pl: "forms/" "del.pl" chk.pl: "forms/" "echo.pl" # File containing the complex-valued seismic moment SOURCE=sourcex.dat # Explosive isotropic source # Data describing the source time function KSIG=2 SIGPLOT='ricker.ps' MPTS=500 SIGT=0.0 SIGF=15. SIGW=2. SIGPH=0. SIGA=1. # Data describing the filtration of the source time function FMIN=0. FLOW=2. FHIGH=120. FMAX=122. # (cosine filter) # Time step and time interval for the Fast Fourier Transform DT=0.001 NFFT=8192 # Data to control seismogram plotting (for details see 'sp.for') KODESP=2 SPCHRH=0.30 SPTMIN=5.0 SPTMAX=0.0 SPTLEN=22.86 SPTDIV=5 SPTSUB=5 SPXMIN=4224. SPXMAX=4224. SPXLEN=15.24 SPXDIV=-280 SPXSUB=1 SPYMAX=10104. SPYMIN=3384. NORMSP=0 SPAMP=3. CALCOPS='0.3 setlinewidth' # Input data for complete ray tracing (see crtin.for) MODEL='sal-mod3.dat' SRC='s145-src.dat' REC='s145oe5r.dat' DCRT='oe-dcr1a.dat' CODE='oe-cod1a.dat' RPAR='oe-rpa1a.dat' WRIT='oe-wri1.dat' CRTLOG='oe5-log.out' INIPAR=3 # Ray tracing of direct wave crt: # Input data for complete ray tracing (see crtin.for) DCRT='oe-dcr1b.dat' CODE='oe-cod1b.dat' RPAR='oe-rpa1b.dat' # Ray tracing of refracted and reflected waves crt: # Input data for complete ray tracing (see crtin.for) DCRT='oe-dcr1c.dat' CODE='oe-cod1c.dat' RPAR='oe-rpa1c.dat' # Ray tracing of refracted and reflected waves crt: # Deleting old files del.pl:'s01.out' del.pl:'s01i.out' del.pl:'s145oe5.mx' # Computed quantities at the given surface and corresponding # quantities at the initial points of rays for all calculated # elementary waves are concatenated. append.pl: 's01.out' 's011.out' append.pl: 's01.out' 's022.out' append.pl: 's01.out' 's033.out' append.pl: 's01.out' 's044.out' append.pl: 's01.out' 's055.out' append.pl: 's01.out' 's026.out' append.pl: 's01.out' 's037.out' append.pl: 's01.out' 's038.out' append.pl: 's01.out' 's029.out' append.pl: 's01.out' 's0210.out' append.pl: 's01.out' 's0211.out' append.pl: 's01.out' 's0212.out' append.pl: 's01.out' 's0213.out' append.pl: 's01.out' 's0214.out' append.pl: 's01.out' 's0215.out' append.pl: 's01.out' 's0216.out' append.pl: 's01.out' 's0217.out' append.pl: 's01.out' 's0218.out' append.pl: 's01.out' 's0219.out' append.pl: 's01.out' 's0220.out' append.pl: 's01.out' 's0221.out' append.pl: 's01.out' 's0222.out' append.pl: 's01.out' 's0323.out' append.pl: 's01.out' 's0224.out' append.pl: 's01.out' 's0225.out' append.pl: 's01.out' 's0226.out' append.pl: 's01.out' 's0227.out' append.pl: 's01.out' 's0228.out' append.pl: 's01.out' 's0229.out' append.pl: 's01.out' 's0230.out' append.pl: 's01i.out' 's011i.out' append.pl: 's01i.out' 's022i.out' append.pl: 's01i.out' 's033i.out' append.pl: 's01i.out' 's044i.out' append.pl: 's01i.out' 's055i.out' append.pl: 's01i.out' 's026i.out' append.pl: 's01i.out' 's037i.out' append.pl: 's01i.out' 's038i.out' append.pl: 's01i.out' 's029i.out' append.pl: 's01i.out' 's0210i.out' append.pl: 's01i.out' 's0211i.out' append.pl: 's01i.out' 's0212i.out' append.pl: 's01i.out' 's0213i.out' append.pl: 's01i.out' 's0214i.out' append.pl: 's01i.out' 's0215i.out' append.pl: 's01i.out' 's0216i.out' append.pl: 's01i.out' 's0217i.out' append.pl: 's01i.out' 's0218i.out' append.pl: 's01i.out' 's0219i.out' append.pl: 's01i.out' 's0220i.out' append.pl: 's01i.out' 's0221i.out' append.pl: 's01i.out' 's0222i.out' append.pl: 's01i.out' 's0323i.out' append.pl: 's01i.out' 's0224i.out' append.pl: 's01i.out' 's0225i.out' append.pl: 's01i.out' 's0226i.out' append.pl: 's01i.out' 's0227i.out' append.pl: 's01i.out' 's0228i.out' append.pl: 's01i.out' 's0229i.out' append.pl: 's01i.out' 's0230i.out' # Frequency-domain seismograms CRTOUT=' ' green: # Green function greenss: # response function SS='oe5-ss.gse' ss: # synthetic seismograms SP1='oe5-sf1.ps' SP2='oe5-sf2.ps' SP3='oe5-sf3.ps' KOLOR=1 sp: # seismogram plotting # Rewriting the two-point rays into a formatted file CRTOUT='crtout.tmp' KALL=1 echo.pl:"'r01.out' ' ' 'r01i.out' ' '" ">crtout.tmp" RAYS='oe5-r1.out' crtray: echo.pl:"'r02.out' ' ' 'r02i.out' ' '" ">crtout.tmp" RAYS='oe5-r2.out' crtray: echo.pl:"'r03.out' ' ' 'r03i.out' ' '" ">crtout.tmp" RAYS='oe5-r3.out' crtray: echo.pl:"'r04.out' ' ' 'r04i.out' ' '" ">crtout.tmp" RAYS='oe5-r4.out' crtray: echo.pl:"'r05.out' ' ' 'r05i.out' ' '" ">crtout.tmp" RAYS='oe5-r5.out' crtray: echo.pl:"'r06.out' ' ' 'r06i.out' ' '" ">crtout.tmp" RAYS='oe5-r6.out' crtray: echo.pl:"'r07.out' ' ' 'r07i.out' ' '" ">crtout.tmp" RAYS='oe5-r7.out' crtray: echo.pl:"'r08.out' ' ' 'r08i.out' ' '" ">crtout.tmp" RAYS='oe5-r8.out' crtray: echo.pl:"'r09.out' ' ' 'r09i.out' ' '" ">crtout.tmp" RAYS='oe5-r9.out' crtray: echo.pl:"'r10.out' ' ' 'r10i.out' ' '" ">crtout.tmp" RAYS='oe5-r10.out' crtray: echo.pl:"'r11.out' ' ' 'r11i.out' ' '" ">crtout.tmp" RAYS='oe5-r11.out' crtray: echo.pl:"'r12.out' ' ' 'r12i.out' ' '" ">crtout.tmp" RAYS='oe5-r12.out' crtray: echo.pl:"'r13.out' ' ' 'r13i.out' ' '" ">crtout.tmp" RAYS='oe5-r13.out' crtray: echo.pl:"'r14.out' ' ' 'r14i.out' ' '" ">crtout.tmp" RAYS='oe5-r14.out' crtray: echo.pl:"'r15.out' ' ' 'r15i.out' ' '" ">crtout.tmp" RAYS='oe5-r15.out' crtray: echo.pl:"'r16.out' ' ' 'r16i.out' ' '" ">crtout.tmp" RAYS='oe5-r16.out' crtray: echo.pl:"'r17.out' ' ' 'r17i.out' ' '" ">crtout.tmp" RAYS='oe5-r17.out' crtray: echo.pl:"'r18.out' ' ' 'r18i.out' ' '" ">crtout.tmp" RAYS='oe5-r18.out' crtray: echo.pl:"'r19.out' ' ' 'r19i.out' ' '" ">crtout.tmp" RAYS='oe5-r19.out' crtray: echo.pl:"'r20.out' ' ' 'r20i.out' ' '" ">crtout.tmp" RAYS='oe5-r20.out' crtray: echo.pl:"'r21.out' ' ' 'r21i.out' ' '" ">crtout.tmp" RAYS='oe5-r21.out' crtray: echo.pl:"'r22.out' ' ' 'r22i.out' ' '" ">crtout.tmp" RAYS='oe5-r22.out' crtray: echo.pl:"'r23.out' ' ' 'r23i.out' ' '" ">crtout.tmp" RAYS='oe5-r23.out' crtray: echo.pl:"'r24.out' ' ' 'r24i.out' ' '" ">crtout.tmp" RAYS='oe5-r24.out' crtray: echo.pl:"'r25.out' ' ' 'r25i.out' ' '" ">crtout.tmp" RAYS='oe5-r25.out' crtray: echo.pl:"'r26.out' ' ' 'r26i.out' ' '" ">crtout.tmp" RAYS='oe5-r26.out' crtray: echo.pl:"'r27.out' ' ' 'r27i.out' ' '" ">crtout.tmp" RAYS='oe5-r27.out' crtray: echo.pl:"'r28.out' ' ' 'r28i.out' ' '" ">crtout.tmp" RAYS='oe5-r28.out' crtray: echo.pl:"'r29.out' ' ' 'r29i.out' ' '" ">crtout.tmp" RAYS='oe5-r29.out' crtray: echo.pl:"'r30.out' ' ' 'r30i.out' ' '" ">crtout.tmp" RAYS='oe5-r30.out' crtray: # Adding two-point rays into GOCAD file VRML='GOCAD' WRL=' ' WRLOUT='oe5-tmp.mx' # GOCAD # Setting blue colour for rays KOLLIN=0 R=0.00 G=0.00 B=1.00 NAME='Rays' PROPERTIES=' ' # just for GOCAD LIN='oe5-r1.out' NAME='Rays1' PROPERTIES=' ' # just for GOCAD linwrl: append.pl: 's145oe5.mx' 'oe5-tmp.mx' # GOCAD LIN='oe5-r2.out' NAME='Rays2' PROPERTIES=' ' # just for GOCAD linwrl: append.pl: 's145oe5.mx' 'oe5-tmp.mx' # GOCAD LIN='oe5-r3.out' NAME='Rays3' PROPERTIES=' ' # just for GOCAD linwrl: append.pl: 's145oe5.mx' 'oe5-tmp.mx' # GOCAD LIN='oe5-r4.out' NAME='Rays4' PROPERTIES=' ' # just for GOCAD linwrl: append.pl: 's145oe5.mx' 'oe5-tmp.mx' # GOCAD LIN='oe5-r5.out' NAME='Rays5' PROPERTIES=' ' # just for GOCAD linwrl: append.pl: 's145oe5.mx' 'oe5-tmp.mx' # GOCAD LIN='oe5-r6.out' NAME='Rays6' PROPERTIES=' ' # just for GOCAD linwrl: append.pl: 's145oe5.mx' 'oe5-tmp.mx' # GOCAD LIN='oe5-r7.out' NAME='Rays7' PROPERTIES=' ' # just for GOCAD linwrl: append.pl: 's145oe5.mx' 'oe5-tmp.mx' # GOCAD LIN='oe5-r8.out' NAME='Rays8' PROPERTIES=' ' # just for GOCAD linwrl: append.pl: 's145oe5.mx' 'oe5-tmp.mx' # GOCAD LIN='oe5-r9.out' NAME='Rays9' PROPERTIES=' ' # just for GOCAD linwrl: append.pl: 's145oe5.mx' 'oe5-tmp.mx' # GOCAD LIN='oe5-r10.out' NAME='Rays10' PROPERTIES=' ' # just for GOCAD linwrl: append.pl: 's145oe5.mx' 'oe5-tmp.mx' # GOCAD LIN='oe5-r11.out' NAME='Rays11' PROPERTIES=' ' # just for GOCAD linwrl: append.pl: 's145oe5.mx' 'oe5-tmp.mx' # GOCAD LIN='oe5-r12.out' NAME='Rays12' PROPERTIES=' ' # just for GOCAD linwrl: append.pl: 's145oe5.mx' 'oe5-tmp.mx' # GOCAD LIN='oe5-r13.out' NAME='Rays13' PROPERTIES=' ' # just for GOCAD linwrl: append.pl: 's145oe5.mx' 'oe5-tmp.mx' # GOCAD LIN='oe5-r14.out' NAME='Rays14' PROPERTIES=' ' # just for GOCAD linwrl: append.pl: 's145oe5.mx' 'oe5-tmp.mx' # GOCAD LIN='oe5-r15.out' NAME='Rays15' PROPERTIES=' ' # just for GOCAD linwrl: append.pl: 's145oe5.mx' 'oe5-tmp.mx' # GOCAD LIN='oe5-r16.out' NAME='Rays16' PROPERTIES=' ' # just for GOCAD linwrl: append.pl: 's145oe5.mx' 'oe5-tmp.mx' # GOCAD LIN='oe5-r17.out' NAME='Rays17' PROPERTIES=' ' # just for GOCAD linwrl: append.pl: 's145oe5.mx' 'oe5-tmp.mx' # GOCAD LIN='oe5-r18.out' NAME='Rays18' PROPERTIES=' ' # just for GOCAD linwrl: append.pl: 's145oe5.mx' 'oe5-tmp.mx' # GOCAD LIN='oe5-r19.out' NAME='Rays19' PROPERTIES=' ' # just for GOCAD linwrl: append.pl: 's145oe5.mx' 'oe5-tmp.mx' # GOCAD LIN='oe5-r20.out' NAME='Rays20' PROPERTIES=' ' # just for GOCAD linwrl: append.pl: 's145oe5.mx' 'oe5-tmp.mx' # GOCAD LIN='oe5-r21.out' NAME='Rays21' PROPERTIES=' ' # just for GOCAD linwrl: append.pl: 's145oe5.mx' 'oe5-tmp.mx' # GOCAD LIN='oe5-r22.out' NAME='Rays22' PROPERTIES=' ' # just for GOCAD linwrl: append.pl: 's145oe5.mx' 'oe5-tmp.mx' # GOCAD LIN='oe5-r23.out' NAME='Rays23' PROPERTIES=' ' # just for GOCAD linwrl: append.pl: 's145oe5.mx' 'oe5-tmp.mx' # GOCAD LIN='oe5-r24.out' NAME='Rays24' PROPERTIES=' ' # just for GOCAD linwrl: append.pl: 's145oe5.mx' 'oe5-tmp.mx' # GOCAD LIN='oe5-r25.out' NAME='Rays25' PROPERTIES=' ' # just for GOCAD linwrl: append.pl: 's145oe5.mx' 'oe5-tmp.mx' # GOCAD LIN='oe5-r26.out' NAME='Rays26' PROPERTIES=' ' # just for GOCAD linwrl: append.pl: 's145oe5.mx' 'oe5-tmp.mx' # GOCAD LIN='oe5-r27.out' NAME='Rays27' PROPERTIES=' ' # just for GOCAD linwrl: append.pl: 's145oe5.mx' 'oe5-tmp.mx' # GOCAD LIN='oe5-r28.out' NAME='Rays28' PROPERTIES=' ' # just for GOCAD linwrl: append.pl: 's145oe5.mx' 'oe5-tmp.mx' # GOCAD LIN='oe5-r29.out' NAME='Rays29' PROPERTIES=' ' # just for GOCAD linwrl: append.pl: 's145oe5.mx' 'oe5-tmp.mx' # GOCAD LIN='oe5-r30.out' NAME='Rays30' PROPERTIES=' ' # just for GOCAD linwrl: append.pl: 's145oe5.mx' 'oe5-tmp.mx' # GOCAD # Adding source and receivers into GOCAD file PTS='s145-src.dat' KOLPTS=0 R=1.00 G=0.00 B=1.00 NAME='Source-145' PROPERTIES=' ' ASIZE=6 # just for GOCAD ptswrl: append.pl: 's145oe5.mx' 'oe5-tmp.mx' # GOCAD PTS='s145oe5r.dat' KOLPTS=0 R=1.00 G=1.00 B=1.00 NAME='Receivers-oe5' PROPERTIES=' ' ASIZE=3 # just for GOCAD ptswrl: append.pl: 's145oe5.mx' 'oe5-tmp.mx' # GOCAD # Conversion of GSE file to SEG-Y files # Input GSE file SS='oe5-ss.gse' # SEG-Y parameters (for details see gse2segy.for) NTRACE=281 NSAMPL=9000 SINTER=0.001 ISFORM=1 ICDFP=1 # Conversion of the first component of GSE file and output SEG-Y file KOMP=1 FISEGY='oe5-ss1.sgy' gse2segy: KOMP=2 FISEGY='oe5-ss2.sgy' gse2segy: KOMP=3 FISEGY='oe5-ss3.sgy' gse2segy: # Travel times # Generating input files CRTOUT for individual elementary waves echo.pl: "' ' 's011.out' 's011i.out' /" "> crt01.out" echo.pl: "' ' 's022.out' 's022i.out' /" "> crt02.out" echo.pl: "' ' 's033.out' 's033i.out' /" "> crt03.out" echo.pl: "' ' 's044.out' 's044i.out' /" "> crt04.out" echo.pl: "' ' 's055.out' 's055i.out' /" "> crt05.out" echo.pl: "' ' 's026.out' 's026i.out' /" "> crt06.out" echo.pl: "' ' 's037.out' 's037i.out' /" "> crt07.out" echo.pl: "' ' 's038.out' 's038i.out' /" "> crt08.out" echo.pl: "' ' 's029.out' 's029i.out' /" "> crt09.out" echo.pl: "' ' 's0210.out' 's0210i.out' /" "> crt10.out" echo.pl: "' ' 's0211.out' 's0211i.out' /" "> crt11.out" echo.pl: "' ' 's0212.out' 's0212i.out' /" "> crt12.out" echo.pl: "' ' 's0213.out' 's0213i.out' /" "> crt13.out" echo.pl: "' ' 's0214.out' 's0214i.out' /" "> crt14.out" echo.pl: "' ' 's0215.out' 's0215i.out' /" "> crt15.out" echo.pl: "' ' 's0216.out' 's0216i.out' /" "> crt16.out" echo.pl: "' ' 's0217.out' 's0217i.out' /" "> crt17.out" echo.pl: "' ' 's0218.out' 's0218i.out' /" "> crt18.out" echo.pl: "' ' 's0219.out' 's0219i.out' /" "> crt19.out" echo.pl: "' ' 's0220.out' 's0220i.out' /" "> crt20.out" echo.pl: "' ' 's0221.out' 's0221i.out' /" "> crt21.out" echo.pl: "' ' 's0222.out' 's0222i.out' /" "> crt22.out" echo.pl: "' ' 's0323.out' 's0323i.out' /" "> crt23.out" echo.pl: "' ' 's0224.out' 's0224i.out' /" "> crt24.out" echo.pl: "' ' 's0225.out' 's0225i.out' /" "> crt25.out" echo.pl: "' ' 's0226.out' 's0226i.out' /" "> crt26.out" echo.pl: "' ' 's0227.out' 's0227i.out' /" "> crt27.out" echo.pl: "' ' 's0228.out' 's0228i.out' /" "> crt28.out" echo.pl: "' ' 's0229.out' 's0229i.out' /" "> crt29.out" echo.pl: "' ' 's0230.out' 's0230i.out' /" "> crt30.out" # Generating files with travel times of individual elementary waves SRC=' ' NQ=4 KALL=0 KREC=1 CRTOUT='crt01.out' PTS='tt01.out' crtpts: CRTOUT='crt02.out' PTS='tt02.out' crtpts: CRTOUT='crt03.out' PTS='tt03.out' crtpts: CRTOUT='crt04.out' PTS='tt04.out' crtpts: CRTOUT='crt05.out' PTS='tt05.out' crtpts: CRTOUT='crt06.out' PTS='tt06.out' crtpts: CRTOUT='crt07.out' PTS='tt07.out' crtpts: CRTOUT='crt08.out' PTS='tt08.out' crtpts: CRTOUT='crt09.out' PTS='tt09.out' crtpts: CRTOUT='crt10.out' PTS='tt10.out' crtpts: CRTOUT='crt11.out' PTS='tt11.out' crtpts: CRTOUT='crt12.out' PTS='tt12.out' crtpts: CRTOUT='crt13.out' PTS='tt13.out' crtpts: CRTOUT='crt14.out' PTS='tt14.out' crtpts: CRTOUT='crt15.out' PTS='tt15.out' crtpts: CRTOUT='crt16.out' PTS='tt16.out' crtpts: CRTOUT='crt17.out' PTS='tt17.out' crtpts: CRTOUT='crt18.out' PTS='tt18.out' crtpts: CRTOUT='crt19.out' PTS='tt19.out' crtpts: CRTOUT='crt20.out' PTS='tt20.out' crtpts: CRTOUT='crt21.out' PTS='tt21.out' crtpts: CRTOUT='crt22.out' PTS='tt22.out' crtpts: CRTOUT='crt23.out' PTS='tt23.out' crtpts: CRTOUT='crt24.out' PTS='tt24.out' crtpts: CRTOUT='crt25.out' PTS='tt25.out' crtpts: CRTOUT='crt26.out' PTS='tt26.out' crtpts: CRTOUT='crt27.out' PTS='tt27.out' crtpts: CRTOUT='crt28.out' PTS='tt28.out' crtpts: CRTOUT='crt29.out' PTS='tt29.out' crtpts: CRTOUT='crt30.out' PTS='tt30.out' crtpts: # Plotting travel times SP1='oe5-tt.ps' SP2=' ' SP3=' ' SS=' ' SS1=' ' SS2=' ' SS3=' ' SS4=' ' SS5=' ' SS6=' ' SS7=' ' SS8=' ' SS9=' ' SS10=' ' SS11=' ' SS12=' ' SS13=' ' SS14=' ' SS15=' ' SS16=' ' SS17=' ' SS18=' ' SS19=' ' SS20=' ' SS21=' ' SS22=' ' SS23=' ' SS24=' ' SS25=' ' SS26=' ' SS27=' ' SS28=' ' SS29=' ' PTS=' ' FTT=' ' SPTTC='oe-ttc.dat' KOLORTD=-1 SPSYMH=0.15 CALCOPS='0.55 setlinewidth' sp: # To save disk space delete temporary output files #del.pl: '*.out' # Important output files # ~~~~~~~~~~~~~~~~~~~~~~ # 'oe5-sf*.ps' ... PostScript files with seismograms # 'oe5-tt.ps' ... PostScript file with travel times # 's145oe5.mx' ... GOCAD file with two-point rays, shot and receivers