makefile2

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

	Interp type=1 d2out=.005 d3out=.005   >$@ maxsize=100


%.A: %.H
	< $*.H Byte gainpanel=all > $@

%-hol3.H %-msk3.H: %.H ${BINDIR}/Holes2.x %.p
	< $*.H ${BINDIR}/Holes2.x par=$*.p mask=$*-msk3.H >$*-hol3.H

%-dip3.H: %.H ${BINDIR}/Smooth3.x 
	< $*.H ${BINDIR}/Smooth3.x nw=3 > $@

qdome-dip.H: qdome-hol3.H qdome-msk3.H ${BINDIR}/Smooth3.x 
	< qdome-hol3.H ${BINDIR}/Smooth3.x nw=3 \
	mask=qdome-msk3.H nliter=10 niter=100 > $@

%-pmis.H: %-hol3.H %-msk3.H %-dip.H ${BINDIR}/Miss3.x %-band.p 
	< $*-hol3.H ${BINDIR}/Miss3.x nw=3 niter=200 \
	mask=$*-msk3.H dip=$*-dip.H | \
	Bandpass par=$*-band.p > $@

%-pdif.H: %-pmis.H %.H %-band.p
	< $*.H Bandpass par=$*-band.p > junk1.H
	Add scale=1,-1 $*-pmis.H junk1.H > junk2.H
	Cat $*-pmis.H junk.H > $@

%-pch3.H: %-dip.H ${BINDIR}/Patch3.x
	< $*-dip.H Window n4=1 f4=1 > q3.H
	< $*-dip.H Window n4=1 | \
	${BINDIR}/Patch3.x np=61 nq=61 pmax=4 qmax=4 q=q3.H > $@

bank5.H nh.H lag.H: ${BINDIR}/Filterbank5.x lage1.p
	${BINDIR}/Filterbank5.x np=61 nq=61 pmax=4 qmax=4 nh=nh.H lags=lag.H \
	nt=40 nx=40 niter=10 eps=0.0009 par=lage1.p > bank5.H

%-mis3.H: %-hol3.H %-msk3.H %-pch3.H bank5.H lag.H nh.H ${BINDIR}/NMiss.x
	< $*-hol3.H ${BINDIR}/NMiss.x niter=100 mask=$*-msk3.H \
	lag=lag.H nh=nh.H filt=bank5.H pch=$*-pch3.H > $@

%-rnd.H: %.H ${BINDIR}/Random3.x %-band.p
	< $*.H ${BINDIR}/Random3.x | Bandpass par=$*-band.p > $@

%-txt3.H: %-pch3.H %-rnd.H bank5.H lag.H nh.H ${BINDIR}/NHelicon.x
	< $*-rnd.H ${BINDIR}/NHelicon.x \
	lag=lag.H nh=nh.H adj=0 div=1 filt=bank5.H pch=$*-pch3.H > $@

pdome%.A: qdome-dip3.H
	< qdome-dip3.H Window n4=1 f4=$* | Byte gainpanel=each pclip=100 > $@

pdome0.v: pdome0.A
	< pdome0.A ${CUBE} color=j title="X slope" out=$@

pdome1.v: pdome1.A
	< pdome1.A ${CUBE} color=j title="Y slope" out=$@

udome%.A: qdome-dip.H
	< qdome-dip.H Window n4=1 f4=$* | Byte gainpanel=each pclip=100 > $@

udome0.v: udome0.A
	< udome0.A ${CUBE} color=j wantscalebar=y title="X slope" out=$@

udome1.v: udome1.A
	< udome1.A ${CUBE} color=j wantscalebar=y title="Y slope" out=$@

${RESDIR}/qslope.v ${RESDIR}/qslope.v3: pdome0.v pdome1.v udome0.v udome1.v
	vp_SideBySideAniso pdome0.v pdome1.v >pdome.v
	vp_SideBySideAniso udome0.v udome1.v >udome.v
	vp_OverUnderAniso pdome.v udome.v > ${RESDIR}/qslope.v
	vp_Movie pdome.v udome.v >  ${RESDIR}/qslope.v3

qdome.v: qdome.A
	< qdome.A ${CUBE} title="Qdome Model" out=$@

qhole.v: qdome-hol3.A
	< qdome-hol3.A ${CUBE} title="60% removed" out=$@

${RESDIR}/qdome.v: qdome.v qhole.v
	vp_SideBySideIso qdome.v qhole.v >$@

${RESDIR}/qmiss.v: qmiss.A qhole.v
	< qmiss.A ${CUBE} title="Interpolated" out=miss.v
	vp_SideBySideIso qhole.v miss.v >$@

qmiss.v: qdome-mis3.A 
	< qdome-mis3.A ${CUBE} title="Interpolated" out=qmiss.v

${RESDIR}/pmiss.v: qdome-pmis.A 
	< qdome-pmis.A ${CUBE} title="Interpolated" out=$@

pmiss.v: qdome-pdif.A
	< qdome-pdif.A Window n3=50 | ${CUBE} title="Interpolated" out=pint.v
	< qdome-pdif.A Window f3=50 | ${CUBE} title="Error" out=pdif.v
	vp_SideBySideAniso pint.v pdif.v > $@

qmiss%.v: qmiss%.A
	< qmiss$*.A ${CUBE} title="Interpolated" out=qmiss$*.v

spike%.H: spike%.p
	Spike par=spike$*.p > $@

mfilt%.H: ${BINDIR}/Plane5.x
	 ${BINDIR}/Plane5.x px=$* py=0.5 rhox=1 rhoy=1 \
	nt=40 nx=20 niter=10 eps=0.001 > $@

mlane%.H: mfilt%.H spike6.H 
	< spike6.H Helicon filt=mfilt$*.H div=1 adj=0 > junk1.H
	< junk1.H Bandpass fhi=0.25 > junk2.H
	< junk2.H  Helicon filt=mfilt$*.H div=1 adj=1 > $@

mlane%.v: mlane%.H
	< mlane$*.H Byte gainpanel=all | ${CUBE2} title="Px=$* Py=0.5" out=$@

${RESDIR}/cube.v ${RESDIR}/cube.v3: \
	mlane-2.v mlane-1.75.v mlane-1.5.v mlane-1.25.v \
	mlane-1.v mlane-0.75.v mlane-0.5.v mlane-0.25.v \
	mlane0.v \
	mlane0.25.v mlane0.5.v mlane0.75.v mlane1.v \
	mlane1.25.v mlane1.5.v mlane1.75.v mlane2.v
	vp_SideBySideIso mlane0.75.v mlane-0.75.v > ${RESDIR}/cube.v
	vp_Movie \
	mlane-2.v mlane-1.75.v mlane-1.5.v mlane-1.25.v \
	mlane-1.v mlane-0.75.v mlane-0.5.v mlane-0.25.v \
	mlane0.v \
	mlane0.25.v mlane0.5.v mlane0.75.v mlane1.v \
	mlane1.25.v mlane1.5.v mlane1.75.v mlane2.v > ${RESDIR}/cube.v3

#################### Constant velocity test #####################

cup.H: ${BINDIR}/Cup.x spike.p
	${BINDIR}/Cup.x par=spike.p | Bandpass flo=10 fhi=50 | \
	Transp | Bandpass fhi=50 | Transp >$@

${RESDIR}/cup.v: cup.H
	< cup.H ${OGREY} title=Model out=$@

data.%.H: %.H ${BINDIR}/PreConstKirch.x
	< $*.H Halfint inv=1 | \
	${BINDIR}/PreConstKirch.x zero=y inv=y \
	h0=0 dh=0.008 nh=61 vel=1.5 > $@
	echo "n4=1 o3=0 d3=0.008 n3=61" >> $@

nmo.%.H: data.%.H ${BINDIR}/ConstNMO.x
	< data.$*.H Transp plane=23 > data23.H
	< data23.H ${BINDIR}/ConstNMO.x v0=1.5  > data32.H
	< data32.H Transp plane=23 > $@

shot.%.H: nmo.%.H ${BINDIR}/Cmp2shot.x
	< nmo.$*.H Transp plane=23 | ${BINDIR}/Cmp2shot.x > $@

synshot.H: synnmod.H ${BINDIR}/Cmp2shot.x
	< synnmod.H ${TRANSP} plane=23 | ${BINDIR}/Cmp2shot.x > $@

beishot.H: beinmod.H ${BINDIR}/Cmp2shot2.x
	< beinmod.H ${TRANSP} plane=23 | ${BINDIR}/Cmp2shot2.x > $@

synshot3.H: synshot.H
	< synshot.H Window n3=3 f3=325 > $@

${RESDIR}/cupdata.v ${RESDIR}/cupdata.v3: nmo.cup.H cuphole.H
	< nmo.cup.H Window n3=1 f3=0  f1=50 n1=100 | \
	${OGREY} title="Zero offset" out=d1.v
	< nmo.cup.H Window n3=1 f3=25 f1=50 n1=100 | \
	${OGREY} title="Middle offset" out=d2.v
	< nmo.cup.H Window n3=1 f3=50 f1=50 n1=100 | \
	${OGREY} title="Far offset" out=d3.v
	vp_OverUnderAniso d1.v d2.v d3.v > d.v
	< cuphole.H Window n3=1 f3=0  f1=50 n1=100 | \
	${OGREY} title="Zero offset" out=n1.v
	< cuphole.H Window n3=1 f3=25 f1=50 n1=100 | \
	${OGREY} title="Middle offset" out=n2.v
	< cuphole.H Window n3=1 f3=50 f1=50 n1=100 | \
	${OGREY} title="Far offset" out=n3.v
	vp_OverUnderAniso n1.v n2.v n3.v > n.v
	vp_Movie d.v n.v > ${RESDIR}/cupdata.v3
	vp_SideBySideAniso d.v n.v > junk.v
	< junk.v vppen \
	vpstyle=n txscale=3 xsize=7.5 ysize=10 > ${RESDIR}/cupdata.v

${RESDIR}/tslice.v ${RESDIR}/tslice.v3: nmo.cup.H cuphole.H
	< nmo.cup.H Window n1=1 f1=74 | \
	${OGREY} transp=n title="t=0.3" label1=midpoint label2=offset out=t1.v
	< nmo.cup.H Window n1=1 f1=99 | \
	${OGREY} transp=n title="t=0.4" label1=midpoint label2=offset out=t2.v
	< nmo.cup.H Window n1=1 f1=124 | \
	${OGREY} transp=n title="t=0.5" label1=midpoint label2=offset out=t3.v
	vp_OverUnderAniso t1.v t2.v t3.v > t.v
	 < cuphole.H Window n1=1 f1=74 | \
	${OGREY} transp=n title="t=0.3" label1=midpoint label2=offset out=n1.v
	< cuphole.H Window n1=1 f1=99 | \
	${OGREY} transp=n title="t=0.4" label1=midpoint label2=offset out=n2.v
	< cuphole.H Window n1=1 f1=124 | \
	${OGREY} transp=n title="t=0.5" label1=midpoint label2=offset out=n3.v
	vp_OverUnderAniso n1.v n2.v n3.v > n.v
	vp_Movie t.v n.v > ${RESDIR}/tslice.v3
	vp_SideBySideAniso t.v n.v > junk.v
	< junk.v vppen \
	vpstyle=n txscale=3 xsize=7.5 ysize=10 > ${RESDIR}/tslice.v


cuphole.H mask.H: nmo.cup.H ${BINDIR}/Shotholes.x 
	< nmo.cup.H ${BINDIR}/Shotholes.x perc=0.9 >cuphole.H mask=mask.H

shothole.H: cuphole.H ${BINDIR}/Cmp2shot.x
	< cuphole.H Transp plane=23 | ${BINDIR}/Cmp2shot.x > $@

shotmask.H: mask.H ${BINDIR}/Cmp2shot.x
	< mask.H Transp plane=23 | Transp | Transp plane=23 | \
	${BINDIR}/Cmp2shot.x | Transp | Transp plane=23 > $@

#cuphole.H mask.H: nmo.cup.H ${BINDIR}/Holes.x 
#	< nmo.cup.H ${BINDIR}/Holes.x perc=0.9 >cuphole.H mask=mask.H

aftd.H: nmo.cup.H ${BINDIR}/Stretch.x ${BINDIR}/FFT1.x ${BINDIR}/Shift.x
	< nmo.cup.H Window f1=25 | \
	${BINDIR}/Stretch.x nout=256 | ${BINDIR}/FFT1.x > fft.H
	< fft.H ${BINDIR}/Shift.x | Transp | Transp plane=23 > $@

shotaftd.H: shot.cup.H ${BINDIR}/Stretch.x ${BINDIR}/FFT1.x ${BINDIR}/Shift.x
	< shot.cup.H Window f1=50 | \
	${BINDIR}/Stretch.x nout=256 | ${BINDIR}/FFT1.x > fft.H
	< fft.H ${BINDIR}/Shift.x | Transp | Transp plane=23 > $@

astck.H: aftd.H ${BINDIR}/Finstack.x
	< aftd.H ${BINDIR}/Finstack.x > $@

admo.H: aftd.H ${BINDIR}/Fincon.x
	< aftd.H Window n2=1 f2=60  | \
	${BINDIR}/Fincon.x h0=0.48 nh=1 dh=-0.008 > $@

astkd.H: astck.H nmo.cup.H ${BINDIR}/FFT2.x ${BINDIR}/Stretch.x
	< astck.H ${TRANSP} | ${BINDIR}/FFT2.x > junk1.H
	< junk1.H ${BINDIR}/Stretch.x inv=1 > junk2.H
	< nmo.cup.H Window n1=25 n3=1 > junk1.H
	Cat axis=1 junk1.H junk2.H > $@

admod.H: admo.H nmo.cup.H ${BINDIR}/FFT2.x ${BINDIR}/Stretch.x
	< admo.H ${TRANSP} plane=23 | ${TRANSP} | ${BINDIR}/FFT2.x > junk1.H
	< junk1.H ${BINDIR}/Stretch.x inv=1 > junk2.H
	< nmo.cup.H Window n1=25 > junk1.H
	Cat axis=1 junk1.H junk2.H > $@

admod1.H: admo.H nmo.cup.H ${BINDIR}/FFT2.x ${BINDIR}/Stretch.x
	< admo.H ${TRANSP} | ${BINDIR}/FFT2.x > junk1.H
	< junk1.H ${BINDIR}/Stretch.x inv=1 > junk2.H
	< nmo.cup.H Window n1=25 n3=1 > junk1.H
	Cat axis=1 junk1.H junk2.H > $@

astup.H: nmo.cup.H
	< nmo.cup.H ${TRANSP} plane=23 | Stack norm=no > $@

acc.H: aftd.H ${BINDIR}/Osmooth.x
	<  aftd.H ${BINDIR}/Osmooth.x eps=0.1 > $@

fftd.H: cuphole.H ${BINDIR}/Stretch.x ${BINDIR}/FFT1.x ${BINDIR}/Shift.x
	< cuphole.H Window f1=50 | \
	${BINDIR}/Stretch.x nout=256 | ${BINDIR}/FFT1.x > fft.H
	< fft.H ${BINDIR}/Shift.x | Transp | Transp plane=23 > $@

shotfftd.H: shothole.H ${BINDIR}/Stretch.x ${BINDIR}/FFT1.x ${BINDIR}/Shift.x
	< shothole.H Window f1=50 | \
	${BINDIR}/Stretch.x nout=256 | ${BINDIR}/FFT1.x > fft.H
	< fft.H ${BINDIR}/Shift.x | Transp | Transp plane=23 > $@

fill.H: fftd.H mask.H ${BINDIR}/Ofilp.x 
	< fftd.H ${BINDIR}/Ofilp.x known=mask.H niter=500 > $@

film.H: fftd.H mask.H ${BINDIR}/Ofilm.x 
	< fftd.H ${BINDIR}/Ofilm.x known=mask.H niter=500 ns=8 > $@

shotfill.H: shotfftd.H shotmask.H ${BINDIR}/Sfil.x 
	< shotfftd.H ${BINDIR}/Sfil.x known=shotmask.H niter=1000 > $@

hfill.H: fftd.H mask.H ${BINDIR}/Hfilp.x 
	< fftd.H ${BINDIR}/Hfilp.x known=mask.H niter=500 > $@

cfill.H: fftd.H mask.H ${BINDIR}/Ofilc.x acc.H 
	< fftd.H ${BINDIR}/Ofilc.x cc=acc.H known=mask.H niter=500 > $@

fftd%.H: fftd.H 
	< fftd.H Window n3=1 f3=$* > $@

ffil%.H: fftd%.H mask.H ${BINDIR}/Ofilp.x 
	< fftd$*.H ${BINDIR}/Ofilp.x known=mask.H niter=500 > $@

shotfftd%.H: shotfftd.H 
	< shotfftd.H Window n3=1 f3=$* > $@

shotaftd%.H: shotaftd.H 
	< shotaftd.H Window n3=1 f3=$* > $@

shot3.H: shotaftd.H
	< shotaftd.H Window n2=3 f2=78 > $@

${RESDIR}/shot3.v ${RESDIR}/shot3.v3: shot.cup.H simp.H
	< shot.cup.H Window n3=1 f3=78 | \
	${OGREY} label2=offset crowd1=0.85 title="shot 1" out=shot1.v
	< shot.cup.H Window n3=1 f3=79 | \
	${OGREY} label2=offset crowd1=0.85 title="shot 2" out=shot2.v
	< shot.cup.H Window n3=1 f3=80 | \
	${OGREY} label2=offset crowd1=0.85 title="shot 3" out=shot3.v
	< simp.H Byte pclip=99.5 > diff.A
	< diff.A Window n3=1 f3=0 | Ta2vplot >/dev/null \
	label1=time label2=offset crowd1=0.85 \
	title="(shot1+shot3)/2" out=diff1.v
	< diff.A Window n3=1 f3=1 | Ta2vplot >/dev/null \
	label1=time label2=offset crowd1=0.85 \
	title="error: (shot1+shot3)/2 - shot2" out=diff2.v
	vp_Movie shot1.v shot2.v shot3.v diff1.v diff2.v > ${RESDIR}/shot3.v3
	vppen vpstyle=n big=n gridnum=2,2 \
	< shot1.v shot3.v diff1.v diff2.v > ${RESDIR}/shot3.v

shotin%.H: shot3.H ${BINDIR}/Infill.x
	< shot3.H ${BINDIR}/Infill.x eps=$* > $@

shotex%.H: shot3.H ${BINDIR}/Inextr.x
	< shot3.H ${BINDIR}/Inextr.x eps=$* > $@

shotexmask.H:
	Spike n1=61 n2=1 > junk1.H
	Spike n1=61 n2=2 k1=-1 > junk2.H
	Cat axis=2 junk1.H junk2.H > $@

shotls.H: shot3.H shotexmask.H ${BINDIR}/Sfil.x 
	< shot3.H Window n2=1 | Transp plane=23 > shot1.H
	Spike n1=61 n2=2 n3=256 k1=-1 | Rtoc > junk2.H
	Cat axis=2 shot1.H junk2.H | \
	${BINDIR}/Sfil.x known=shotexmask.H niter=1000 > $@

shotfill%.H: shotfftd%.H shotmask.H ${BINDIR}/Sfil.x 
	< shotfftd$*.H ${BINDIR}/Sfil.x known=shotmask.H niter=1000 > $@

${RESDIR}/fslice.v ${RESDIR}/fslice.v3: fill.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