makefile2

国外免费地震资料处理软件包

	vp_OverUnderAniso f1.v f2.v f3.v > f.v
	< fill.H Window n3=1 f3=120 | Real | \
	${OGREY} pclip=99 title="frequency=5" \
	transp=n label1=midpoint label2=offset out=i1.v
	< fill.H Window n3=1 f3=112 | Real | \
	${OGREY} pclip=99 title="frequency=10" \
	transp=n label1=midpoint label2=offset out=i2.v
	< fill.H Window n3=1 f3=103 | Real | \
	${OGREY} pclip=99 title="frequency=15" \
	transp=n label1=midpoint label2=offset out=i3.v
	vp_OverUnderAniso i1.v i2.v i3.v > i.v
	vp_SideBySideAniso f.v i.v > junk.v
	vp_Movie f.v i.v > ${RESDIR}/fslice.v3
	< junk.v vppen vpstyle=n txscale=3 xsize=7.5 ysize=10 \
	> ${RESDIR}/fslice.v

${RESDIR}/gslice.v ${RESDIR}/gslice.v3: film.H fftd.H
	< fftd.H Window n3=1 f3=120 | Real | \
	${OGREY} pclip=99 title="frequency=5" \
	transp=n label1=midpoint label2=offset out=f1.v
	< fftd.H Window n3=1 f3=112 | Real | \
	${OGREY} pclip=99 title="frequency=10" \
	transp=n label1=midpoint label2=offset out=f2.v
	< fftd.H Window n3=1 f3=103 | Real | \
	${OGREY} pclip=99 title="frequency=15" \
	transp=n label1=midpoint label2=offset out=f3.v
	vp_OverUnderAniso f1.v f2.v f3.v > f.v
	< film.H Window n3=1 f3=120 | Real | \
	${OGREY} pclip=99 title="frequency=5" \
	transp=n label1=midpoint label2=offset out=i1.v
	< film.H Window n3=1 f3=112 | Real | \
	${OGREY} pclip=99 title="frequency=10" \
	transp=n label1=midpoint label2=offset out=i2.v
	< film.H Window n3=1 f3=103 | Real | \
	${OGREY} pclip=99 title="frequency=15" \
	transp=n label1=midpoint label2=offset out=i3.v
	vp_OverUnderAniso i1.v i2.v i3.v > i.v
	vp_SideBySideAniso f.v i.v > junk.v
	vp_Movie f.v i.v > ${RESDIR}/gslice.v3
	< junk.v vppen vpstyle=n txscale=3 xsize=7.5 ysize=10 \
	> ${RESDIR}/gslice.v

${RESDIR}/sslice.v ${RESDIR}/sslice.v3: shotfill.H shotfftd.H
	< shotfftd.H Window n3=1 f3=120 | Real | \
	${OGREY} pclip=99 title="frequency=5" \
	transp=n label1=midpoint label2=offset out=f1.v
	< shotfftd.H Window n3=1 f3=112 | Real | \
	${OGREY} pclip=99 title="frequency=10" \
	transp=n label1=midpoint label2=offset out=f2.v
	< shotfftd.H Window n3=1 f3=103 | Real | \
	${OGREY} pclip=99 title="frequency=15" \
	transp=n label1=midpoint label2=offset out=f3.v
	vp_OverUnderAniso f1.v f2.v f3.v > f.v
	< shotfill.H Window n3=1 f3=120 | Real | \
	${OGREY} pclip=99 title="frequency=5" \
	transp=n label1=midpoint label2=offset out=i1.v
	< shotfill.H Window n3=1 f3=112 | Real | \
	${OGREY} pclip=99 title="frequency=10" \
	transp=n label1=midpoint label2=offset out=i2.v
	< shotfill.H Window n3=1 f3=103 | Real | \
	${OGREY} pclip=99 title="frequency=15" \
	transp=n label1=midpoint label2=offset out=i3.v
	vp_OverUnderAniso i1.v i2.v i3.v > i.v
	vp_SideBySideAniso f.v i.v > junk.v
	vp_Movie f.v i.v > ${RESDIR}/sslice.v3
	< junk.v vppen vpstyle=n txscale=3 xsize=7.5 ysize=10 \
	> ${RESDIR}/sslice.v

${RESDIR}/hslice.v ${RESDIR}/hslice.v3: hfill.H fftd.H
	< fftd.H Window n3=1 f3=120 | Real | \
	${OGREY} pclip=99 title="frequency=5" \
	transp=n label1=midpoint label2=offset out=f1.v
	< fftd.H Window n3=1 f3=112 | Real | \
	${OGREY} pclip=99 title="frequency=10" \
	transp=n label1=midpoint label2=offset out=f2.v
	< fftd.H Window n3=1 f3=103 | Real | \
	${OGREY} pclip=99 title="frequency=15" \
	transp=n label1=midpoint label2=offset out=f3.v
	vp_OverUnderAniso f1.v f2.v f3.v > f.v
	< hfill.H Window n3=1 f3=120 | Real | \
	${OGREY} pclip=99 title="frequency=5" \
	transp=n label1=midpoint label2=offset out=i1.v
	< hfill.H Window n3=1 f3=112 | Real | \
	${OGREY} pclip=99 title="frequency=10" \
	transp=n label1=midpoint label2=offset out=i2.v
	< hfill.H Window n3=1 f3=103 | Real | \
	${OGREY} pclip=99 title="frequency=15" \
	transp=n label1=midpoint label2=offset out=i3.v
	vp_OverUnderAniso i1.v i2.v i3.v > i.v
	vp_SideBySideAniso f.v i.v > junk.v
	vp_Movie f.v i.v > ${RESDIR}/hslice.v3
	< junk.v vppen vpstyle=n txscale=3 xsize=7.5 ysize=10 \
	> ${RESDIR}/hslice.v

${RESDIR}/cslice.v ${RESDIR}/cslice.v3: cfill.H fftd.H
	< fftd.H Window n3=1 f3=120 | Real | \
	${OGREY} pclip=99 title="frequency=5" \
	transp=n label1=midpoint label2=offset out=f1.v
	< fftd.H Window n3=1 f3=112 | Real | \
	${OGREY} pclip=99 title="frequency=10" \
	transp=n label1=midpoint label2=offset out=f2.v
	< fftd.H Window n3=1 f3=103 | Real | \
	${OGREY} pclip=99 title="frequency=15" \
	transp=n label1=midpoint label2=offset out=f3.v
	vp_OverUnderAniso f1.v f2.v f3.v > f.v
	< cfill.H Window n3=1 f3=120 | Real | \
	${OGREY} pclip=99 title="frequency=5" \
	transp=n label1=midpoint label2=offset out=i1.v
	< cfill.H Window n3=1 f3=112 | Real | \
	${OGREY} pclip=99 title="frequency=10" \
	transp=n label1=midpoint label2=offset out=i2.v
	< cfill.H Window n3=1 f3=103 | Real | \
	${OGREY} pclip=99 title="frequency=15" \
	transp=n label1=midpoint label2=offset out=i3.v
	vp_OverUnderAniso i1.v i2.v i3.v > i.v
	vp_SideBySideAniso f.v i.v > junk.v
	vp_Movie f.v i.v > ${RESDIR}/cslice.v3
	< junk.v vppen vpstyle=n txscale=3 xsize=7.5 ysize=10 \
	> ${RESDIR}/cslice.v



fild.H: nmo.cup.H fill.H \
	${BINDIR}/Stretch.x ${BINDIR}/FFT2.x ${BINDIR}/Shift.x 
	< fill.H Transp plane=23 | Transp | \
	${BINDIR}/Shift.x inv=y> transp.H
	< transp.H ${BINDIR}/FFT2.x > finv.H
	< finv.H ${BINDIR}/Stretch.x dens=2 inv=1 > junk1.H
	< nmo.cup.H Window n1=50 > junk2.H
	Cat axis=1 junk2.H junk1.H > $@

fimd.H: nmo.cup.H film.H \
	${BINDIR}/Stretch.x ${BINDIR}/FFT2.x ${BINDIR}/Shift.x 
	< film.H Transp plane=23 | Transp | \
	${BINDIR}/Shift.x inv=y> transp.H
	< transp.H ${BINDIR}/FFT2.x > finv.H
	< finv.H ${BINDIR}/Stretch.x dens=2 inv=1 > junk1.H
	< nmo.cup.H Window n1=50 > junk2.H
	Cat axis=1 junk2.H junk1.H > $@

shotfild.H: shot.cup.H shotfill.H \
	${BINDIR}/Stretch.x ${BINDIR}/FFT2.x ${BINDIR}/Shift.x 
	< shotfill.H Transp plane=23 | Transp | \
	${BINDIR}/Shift.x inv=y> transp.H
	< transp.H ${BINDIR}/FFT2.x > finv.H
	< finv.H ${BINDIR}/Stretch.x dens=2 inv=1 > junk1.H
	< shot.cup.H Window n1=50 > junk2.H
	Cat axis=1 junk2.H junk1.H > $@

shotxd%.H: shot.cup.H shotex%.H \
	${BINDIR}/Stretch.x ${BINDIR}/FFT2.x ${BINDIR}/Shift.x 
	< shotex$*.H Transp ${BINDIR}/Shift.x inv=y> transp.H
	< transp.H ${BINDIR}/FFT2.x > finv.H
	< finv.H ${BINDIR}/Stretch.x dens=2 inv=1 > junk1.H
	< shot.cup.H Window n1=50 n3=1 f3=79 > junk2.H
	Cat axis=1 junk2.H junk1.H > $@

shotlsd.H: shot.cup.H shotls.H \
	${BINDIR}/Stretch.x ${BINDIR}/FFT2.x ${BINDIR}/Shift.x 
	< shotls.H Transp plane=23 | Transp | \
	${BINDIR}/Shift.x inv=y> transp.H
	< transp.H ${BINDIR}/FFT2.x > finv.H
	< finv.H ${BINDIR}/Stretch.x dens=2 inv=1 > junk1.H
	< shot.cup.H Window n1=50 n3=3 f3=78 > junk2.H
	Cat axis=1 junk2.H junk1.H > $@

shotnd%.H: shot.cup.H shotin%.H \
	${BINDIR}/Stretch.x ${BINDIR}/FFT2.x ${BINDIR}/Shift.x 
	< shotin$*.H Transp | ${BINDIR}/Shift.x inv=y> transp.H
	< transp.H ${BINDIR}/FFT2.x > finv.H
	< finv.H ${BINDIR}/Stretch.x dens=2 inv=1 > junk1.H
	< shot.cup.H Window n1=50 n3=1 f3=79 > junk2.H
	Cat axis=1 junk2.H junk1.H > $@

shotbench.H: shot.cup.H
	< shot.cup.H Window n3=1 f3=79 > $@

shotdiff%.H: shotnd%.H shotbench.H
	Add scale=1,-1 shotnd$*.H shotbench.H > diff.H
	Cat shotnd$*.H diff.H > $@

${RESDIR}/shotin.v: shotdiff0.H shotdiff0.5.H
	< shotdiff0.H Byte pclip=99.5 > diff1.A
	< shotdiff0.5.H Byte pclip=99.5 > diff2.A
	< diff1.A Window n3=1 f3=0 | Ta2vplot >/dev/null \
	label1=time label2=offset crowd1=0.85 \
	title="interpolated shot, no reg" out=shot1.v
	< diff1.A Window n3=1 f3=1 | Ta2vplot >/dev/null \
	label1=time label2=offset crowd1=0.85 \
	title="error, no reg" out=diff1.v
	< diff2.A Window n3=1 f3=0 | Ta2vplot >/dev/null \
	label1=time label2=offset crowd1=0.85 \
	title="interpolated shot, reg" out=shot2.v
	< diff2.A Window n3=1 f3=1 | Ta2vplot >/dev/null \
	label1=time label2=offset crowd1=0.85 \
	title="error, reg" out=diff2.v
	vppen vpstyle=n big=n gridnum=2,2 < \
	shot1.v diff1.v shot2.v diff2.v > ${RESDIR}/shotin.v

sfild.H: shotfild.H ${BINDIR}/Shot2cmp.x
	< shotfild.H ${BINDIR}/Shot2cmp.x > $@

hild.H: nmo.cup.H hfill.H \
	${BINDIR}/Stretch.x ${BINDIR}/FFT2.x ${BINDIR}/Shift.x 
	< hfill.H Transp plane=23 | Transp | \
	${BINDIR}/Shift.x inv=y> transp.H
	< transp.H ${BINDIR}/FFT2.x > finv.H
	< finv.H ${BINDIR}/Stretch.x dens=2 inv=1 > junk1.H
	< nmo.cup.H Window n1=50 > junk2.H
	Cat axis=1 junk2.H junk1.H > $@

all.H: nmo.cup.H cuphole.H hild.H fild.H fimd.H
	Merge axis=2 nmo.cup.H cuphole.H hild.H fild.H fimd.H > $@

all.v3: all.H
	< all.H Window f1=50 n1=100 > junk.H
	< junk.H Grey crowd=1 >/dev/null out=$@

${RESDIR}/all.v ${RESDIR}/all.v3: hild.H fild.H 
	< hild.H Window n3=1  f3=0 f1=50 n1=100 | \
	${OGREY} title="Near offset" out=h1.v
	< hild.H Window n3=1 f3=25 f1=50 n1=100 | \
	${OGREY} title="Middle offset" out=h2.v
	< hild.H Window n3=1 f3=50 f1=50 n1=100 | \
	${OGREY} title="Far offset" out=h3.v
	vp_OverUnderAniso h1.v h2.v h3.v > h.v
	< fild.H Window n3=1 f3=0  f1=50 n1=100 | \
	${OGREY} title="Near offset" out=f1.v
	< fild.H Window n3=1 f3=25 f1=50 n1=100 | \
	${OGREY} title="Middle offset" out=f2.v
	< fild.H Window n3=1 f3=50 f1=50 n1=100 | \
	${OGREY} title="Far offset" out=f3.v
	vp_OverUnderAniso f1.v f2.v f3.v > f.v
	vp_Movie h.v f.v > ${RESDIR}/all.v3
	vp_SideBySideAniso h.v f.v > junk.v
	< junk.v vppen \
	vpstyle=n txscale=3 xsize=7.5 ysize=10 > ${RESDIR}/all.v

################## Marmousi test #########################################

ELF = /net/kana/data_3d/elf_north_sea/data3950.H

elf.H: ${ELF} ${BINDIR}/Clip0.x ${BINDIR}/Zero.x
	< ${ELF} ${BINDIR}/Clip0.x clip=100000 > clip.H
	echo "d2=12.5 o2=50 hff=-1" >> clip.H
	< clip.H ${BINDIR}/Zero.x f2=67 f3=663 nt=300 > zero.H
	< zero.H Window n1=800 > $@

${RESDIR}/elf.v: elf.H
	< elf.H Window n2=1 f2=0 > near.H
	< elf.H Window n2=1 f2=119  > far.H
	< near.H ${OGREY} title="Near offset" out=near.v
	< far.H ${OGREY}  title="Far offset" out=far.v
	vp_SideBySideAniso near.v far.v > junk.v
	< junk.v vppen vpstyle=n xsize=12 ysize=6 > $@

${RESDIR}/marm.v: ${MARMDIR}/marmrefl.H
	< ${MARMDIR}/marmrefl.H Window n2=1 f2=95 > near.H
	< ${MARMDIR}/marmrefl.H Window n2=1 f2=0  > far.H
	< near.H ${OGREY} title="Near offset" out=near.v
	< far.H ${OGREY}  title="Far offset" out=far.v
	vp_SideBySideAniso near.v far.v > junk.v
	< junk.v vppen vpstyle=n xsize=12 ysize=6 > $@

marmcmp.H: ${MARMDIR}/marmrefl.H ${BINDIR}/Shot2cmp.x
	< ${MARMDIR}/marmrefl.H Reverse which=2 | ${BINDIR}/Shot2cmp.x > $@

marmpad.H: marmcmp.H
	< marmcmp.H Reverse which=2 | Pad n2out=104 | Reverse which=2 > $@

marmshot.H: marmpad.H ${BINDIR}/Cmp2shot.x
	< marmpad.H ${BINDIR}/Cmp2shot.x > $@

mshot%.H: marmshot.H ${BINDIR}/Extractshot.x
	< marmshot.H ${BINDIR}/Extractshot.x is=$* > $@

marmint.H: marmfild.H ${BINDIR}/Cmp2shot.x
	< marmfild.H ${BINDIR}/Cmp2shot.x > $@

ishot%.H: marmint.H ${BINDIR}/Extractshot.x ${BINDIR}/Mutter.x
	< marmint.H ${BINDIR}/Extractshot.x is=$* | \
	${BINDIR}/Mutter.x tp=0.001 slope0=0.001 > $@

${RESDIR}/mshot.v ${RESDIR}/mshot.v3: mshot224.H ishot224.H
	< mshot224.H ${OGREY} pclip=99 title=Input out=m.v
	< ishot224.H ${OGREY} pclip=99 title=Interpolated out=i.v
	vp_Movie m.v i.v > ${RESDIR}/mshot.v3
	vp_OverUnderAniso m.v i.v > junk.v
	< junk.v vppen vpstyle=n xsize=7.5 ysize=10 > ${RESDIR}/mshot.v

${RESDIR}/mishot.v: ishot225.H
	< ishot224.H ${OGREY} pclip=99 title=Interpolated out=i.v
	< i.v vppen vpstyle=n xsize=7.5 ysize=5 > $@

mcmp%.H: marmpad.H
	< marmpad.H Window n3=1 f3=$* > cmp.H
	< cmp.H Reverse which=2 | Window n2=103 > rcmp.H
	Cat rcmp.H cmp.H axis=2 > $@

icmp%.H: marmfild.H ${BINDIR}/Mutter.x
	< marmfild.H Window n3=1 f3=$* > cmp.H
	< cmp.H Reverse which=2 | Window n2=103 > rcmp.H
	Cat rcmp.H cmp.H axis=2 | ${BINDIR}/Mutter.x tp=0.001 slope0=0.001 > $@

${RESDIR}/mcmp.v ${RESDIR}/mcmp.v3: mcmp224.H icmp224.H
	< mcmp224.H ${OGREY} pclip=99 title=Input out=m.v
	< icmp224.H ${OGREY} pclip=99 title=Interpolated out=i.v
	vp_Movie m.v i.v > ${RESDIR}/mcmp.v3
	vp_OverUnderAniso m.v i.v > junk.v
	< junk.v vppen vpstyle=n xsize=7.5 ysize=10 > ${RESDIR}/mcmp.v

marmmask.H: marmpad.H ${BINDIR}/Holes.x
	<marmpad.H ${BINDIR}/Holes.x perc=0 >junks.H mask=junkmask.H
	<junkmask.H Transp > $@

marmnmo.H: marmpad.H ${BINDIR}/ConstNMO.x marm.p 
	< marmpad.H ${BINDIR}/ConstNMO.x par=marm.p > marmnmo.H	

elfnmo.H: elf.H ${BINDIR}/ConstNMO.x elf.p 
	< elf.H ${BINDIR}/ConstNMO.x par=elf.p > $@

marmsmb.H: marmnmo.H ${BINDIR}/Vscan.x marm.p 
	< marmnmo.H ${BINDIR}/Vscan.x par=marm.p semblance=y > $@

elfsmb.H: elfnmo.H ${BINDIR}/Vscan.x elf.p 
	< elfnmo.H ${BINDIR}/Vscan.x par=elf.p semblance=y > $@

marmvpk.H: marmsmb.H ${BINDIR}/Blindpick.x
	< marmsmb.H Reverse which=2 | \
	${BINDIR}/Blindpick.x eps=0.1 lam=0.1 | \
	${TRANSP} plane=23 > $@

marmnmod.H: marmpad.H ${DATADIR}/marmvel.H ${BINDIR}/NMO.x ${BINDIR}/Mutter.x
	<  marmpad.H ${BINDIR}/Mutter.x slope0=0.0018 > marmmut.H
	< marmmut.H ${BINDIR}/NMO.x vel=${DATADIR}/marmvel.H > $@

beinmod.H: ${DATADIR}/bei.H ${DATADIR}/beivel.H \
	${BINDIR}/NMO.x ${BINDIR}/Mutter.x
	< ${DATADIR}/bei.H Transp plane=23 > bei23.H
	< bei23.H ${BINDIR}/Mutter.x slope0=1.4 > beimut.H
	< beimut.H ${BINDIR}/NMO.x vel=${DATADIR}/beivel.H > bei23.H