makefile2

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

include ${SEPINC}/SEP.top
RESDIR = ./Fig
HTMLDIR = .

RESULTSER = otens int splin mesh gal cross \
	    hole seab passfill \
	    sean2 sean2-dip sean2-int sean2-close \
	    qdome cube \
	    cup cupdata tslice fslice all \
	    shot3 shotin elfshot3 elfshotin \
	    sthree bob aliasp2 off-imp fold-win
RESULTSCR = specc specp phase qslope pmiss arg \
	    bin-win smo2-win int2-win int4-win winslice off4
RESULTSNR = 

UF90LIBS = ${GEELIB}
UF77LIBS = ${BEILIB} ${GEMLIB} -lsepfftf -lsepauxf

DOTS = Dots dots=2 connect=0 gaineach=0 overlap=1. >/dev/null 
GREY = Grey gainpanel=every pclip=100 wantaxis=n >/dev/null crowd=.85 gpow=.7 
PLOT = Grey crowd=.85  wantaxis=0 >/dev/null
SEABGREY = Grey crowd=0.8 pclip=100 >/dev/null \
	transp=n yreverse=n label1=longitude label2=latitude
OGREY = Grey >/dev/null label1=time label2=midpoint pclip=99.5
SGREY = Grey pclip=100 wantaxis=n wheretitle=b titlesz=18 >/dev/null crowd=.88
CUBE = Cubeplot frame1=60 frame2=60 frame3=24 \
       flat=y point1=0.66 point2=0.66 >/dev/null
CUBE = Cubeplot frame1=100  frame2=150 frame3=100 \
       flat=y point1=0.6 point2=0.6 >/dev/null
CUBE2 = Cubeplot frame1=20 frame2=20 frame3=20 wantaxes=0 \
       flat=y point1=0.5 point2=0.5 wantcoordlabel=-1 >/dev/null
CUBE1 = Cubeplot flat=y point1=0.75 point2=0.75  \
        frame1=164 frame2=142 frame3=56 >/dev/null
CUBE3 = ${TRANSP} plane=23 | ${CUBE1}

TRANSP = Transp verb=1 max_memory=100

BOOKDIR = ${SEP}/book

DATDIR = ${BOOKDIR}/gee/Data
MARMDIR = ${DATALIB}/oldq2/marmousi
DATADIR = ../Dat

BINDIR = ../Bin/${MTYPE}
OBJDIR = ../Obj/${MTYPE}
SRCDIR = ../Src

LATOPTS = style=amsmath

Makefile : ;

################## 1-D figures #############################################

in.H: in.HH
	< in.HH Dd esize=4 | Window > $@

otens%.H: Spline1.x
	Spline1.x tension=$* > $@

ltens%.H: otens%.H Plus1.x
	< otens$*.H Plus1.x > $@

${RESDIR}/otens.v: ltens0.01.H ltens0.25.H ltens0.5.H ltens0.75.H ltens1.H
	Cat axis=2 ltens1.H ltens0.75.H ltens0.5.H ltens0.25.H ltens0.01.H | \
	${DOTS} out=$@ \
	labels=tension=1:tension=0.75:tension=0.5:tension=0.25:tension=0.01 

int%.H: in.H otens%.H
	< in.H Miss filt=otens$*.H > $@

in.v: in.H
	< in.H ${DOTS} dots=0 labels=given out=$@

int%.v: int%.H 
	<int$*.H ${DOTS} labels=tension=$* out=$@

${RESDIR}/int.v: ${RESDIR}/int.v3 int0.01.H int0.25.H int0.5.H int0.75.H int1.H
	Merge axis=2 space=n  < \
	int1.H int0.75.H int0.5.H int0.25.H int0.01.H in.H > two.H
	< two.H ${DOTS} out=junk.v corners=12 connect=4 radius=0.02 \
	labels=tension=1:tension=0.75:tension=0.5:tension=0.25:tension=0.01:given 
	< junk.v vppen vpstyle=n xsize=6 ysize=8 > $@

${RESDIR}/int.v3: in.v int0.01.v int0.25.v int0.5.v int0.75.v int1.v
	vp_Movie in.v int0.01.v int0.25.v int0.5.v int0.75.v int1.v > $@

##################### 2-D figures #########################################

# 53 -> 35 (43)
stens%.H: ${BINDIR}/Spline.x
	${BINDIR}/Spline.x tension=$* eps=1.e-3 > $@

# 25 -> 16 (20)
gtens%.H: ${BINDIR}/Spline.x
	${BINDIR}/Spline.x tension=$* eps=5.e-3 > $@

splin%.v : gtens%.H 
	Spike n1=40 n2=40 nsp=2 k1=11,16 k2=8,3 mag=1,-1  > spike.H
	Spike n1=40 n2=40 nsp=2 k1=31,28 k2=24,16 mag=1,-1 > inp2.H
	<spike.H Helicon filt=gtens$*.H > inp1.H
	Add inp1.H inp2.H > inp.H
	<inp.H Helicon       div=1 filt=gtens$*.H > div.H
	<div.H Helicon adj=1 div=1 filt=gtens$*.H > div2.H
	<inp.H  ${GREY} out=inp.v  clip=1 title=tension=$*
	<div.H  ${GREY} out=div.v  clip=1 title=input/filter
	<div2.H ${GREY} out=div2.v pclip=100 title="(input/filter)/filter'"
	vp_SideBySideAniso inp.v div.v div2.v > $@

${RESDIR}/splin.v ${RESDIR}/splin.v3: splin0.v splin0.3.v splin0.7.v splin1.v
	vp_OverUnderAniso splin0.3.v splin0.7.v splin1.v > junk.v
	<junk.v  vppen vpstyle=n txscale=4 > ${RESDIR}/splin.v
	vp_Movie splin0.v splin0.3.v splin0.7.v splin1.v > ${RESDIR}/splin.v3

##################### Galilee ##############################################

gal.H: ${DATDIR}/galilee.H ${BINDIR}/Clip.x
	Cp ${DATDIR}/galilee.H copy.H
	< copy.H Transp > transp.H
	<transp.H ${BINDIR}/Clip.x clip=-212. >gal.H

mesh.H: gal.H ${BINDIR}/Bin1.x galilee.p
	<  gal.H ${BINDIR}/Bin1.x par=galilee.p interp=bin2 > mesh.H

fill%.H:   mesh.H ${BINDIR}/Missgal.x gtens%.H
	< mesh.H ${BINDIR}/Missgal.x eps=0.1 filt=gtens$*.H > fill$*.H niter=20

${RESDIR}/mesh.v: mesh.H
	< mesh.H  ${GREY} allpos=y \
	pclip=94 wantaxis=y yreverse=n transp=n \
	out=junk.v title="Binned" label2="South-North (km)" label1="West-East (km)"
	< junk.v vppen vpstyle=n xsize=3 ysize=4 > $@

gal%.v: fill%.H ${BINDIR}/Igrad.x rect.p
	< fill$*.H ${BINDIR}/Igrad.x   par=rect.p > ruf$*.H
	< ruf$*.H  ${GREY} \
	wantaxis=y pclip=94  yreverse=n transp=n \
	out=gal$*.v title="tension=$*"

cross%.H: fill%.H 
	< fill$*.H Window n2=1 f2=300 n1=440 f1=80 > cross$*.H

${RESDIR}/cross.v: cross0.H cross0.3.H cross0.7.H cross1.H
	Cat axis=2 cross1.H cross0.7.H cross0.3.H cross0.H > cross.H
	< cross.H Scale dscale=-1 | ${DOTS} dots=0 connect=4 out=$@ \
	labels=tension=1:tension=0.7:tension=0.3:tension=0 out=junk.v \
	label1="West-East (km)"
	< junk.v vppen vpstyle=n xsize=6 ysize=6 > $@

${RESDIR}/gal.v3: gal0.v gal0.25.v gal0.5.v gal0.75.v gal1.v
	vp_Movie gal0.v gal0.25.v gal0.5.v gal0.75.v gal1.v > $@

${RESDIR}/gal.v: ${RESDIR}/gal.v3 gal0.v gal0.3.v gal0.7.v gal1.v
	vp_SideBySideAniso gal0.v gal0.3.v > top.v
	vp_SideBySideAniso gal0.7.v gal1.v > bot.v
	vp_OverUnderAniso top.v bot.v > junk.v
	< junk.v vppen vpstyle=n txscale=2 xsize=6 ysize=8 > $@

make2.H: ${BINDIR}/Make.x
	${BINDIR}/Make.x  n1=100 n2=20 > make2.H

hole.H: make2.H ${BINDIR}/Hole.x
	<make2.H ${BINDIR}/Hole.x > hole.H

pqhole.H: hole.H ${BINDIR}/Smooth2.x
	< hole.H Window | ${BINDIR}/Smooth2.x mask=hole.H niter=100 > pqhole.H

pqmiss.H: hole.H ${BINDIR}/Missp.x pqhole.H
	< hole.H ${BINDIR}/Missp.x dip=pqhole.H niter=100 > $@

${RESDIR}/hole.v:  make2.H hole.H pqmiss.H pqhole.H
	<make2.H ${PLOT} title=original out=make.v
	<hole.H ${PLOT} title=gapped   out=hole.v
	< pqhole.H Window n3=1 f3=0 | ${PLOT} title="first dip" out=dip1.v
	< pqhole.H Window n3=1 f3=1 | ${PLOT} title="second dip" out=dip2.v
	<pqmiss.H ${PLOT} title=restored out=miss.v
	vp_SideBySideAniso make.v hole.v dip1.v dip2.v miss.v > junk.v
	<junk.v  vppen vpstyle=n txscale=2. >${RESDIR}/hole.v

bin.H: ${DATDIR}/apr18.H interp2.p ${BINDIR}/Binsea.x 
	Cp ${DATDIR}/apr18.H copy.H
	< copy.H  Transp > wind18.H
	<wind18.H ${BINDIR}/Binsea.x par=interp2.p > bin.H

pbin.H: bin.H ${BINDIR}/Smooth.x  
	< bin.H ${BINDIR}/Smooth.x nw=2 mask=bin.H niter=1000 > $@

seab.H: bin.H pbin.H ${BINDIR}/Miss1.x
	< bin.H ${BINDIR}/Miss1.x nw=2 dip=pbin.H niter=1000 > $@

${RESDIR}/seab.v ${RESDIR}/seab.v3: bin.H seab.H
	< bin.H Window \
	min1=-113.11 max1=-112.5736 min2=-16.1474 max2=-15.6086 | \
	${SEABGREY} out=bin.v title=Binned
	< seab.H Bandpass fhi=90 | Window \
	min1=-113.11 max1=-112.5736 min2=-16.1474 max2=-15.6086 | \
	${SEABGREY} out=int.v title=Interpolated
	vp_SideBySideAniso bin.v int.v > ${RESDIR}/seab.v
	vp_Movie bin.v int.v > ${RESDIR}/seab.v3

pass-dat.H: ${DATDIR}/blast.46.H 
	<${DATDIR}/blast.46.H Window f1=150 n1=150  > $@

pass-mask.H: pass-dat.H ${BINDIR}/Mask3.x
	< pass-dat.H ${BINDIR}/Mask3.x eps=0.00001 > $@

pass-dip.H: pass-dat.H pass-mask.H ${BINDIR}/Smooth3.x 
	< pass-dat.H ${BINDIR}/Smooth3.x nw=2 \
	mask=pass-mask.H nliter=10 > $@

pass-mis.H: pass-dat.H pass-mask.H pass-dip.H ${BINDIR}/Miss3.x 
	< pass-dat.H ${BINDIR}/Miss3.x nw=2 niter=200 \
	mask=pass-mask.H dip=pass-dip.H | Bandpass fhi=50 > $@

WIND = Window n3=1
TOVPLOT = | Grey fastplot=20 crowd=.95 wantaxis=n wanttitle=n > /dev/null
${RESDIR}/passfill.v ${RESDIR}/passfill.v3: pass-dat.H pass-mis.H 
	< pass-dat.H Byte clip=.0007 > data.A
	< data.A ${WIND} f3=00  ${TOVPLOT} out=d00.v
	< data.A ${WIND} f3=01  ${TOVPLOT} out=d01.v
	< data.A ${WIND} f3=02  ${TOVPLOT} out=d02.v
	< data.A ${WIND} f3=03  ${TOVPLOT} out=d03.v
	< data.A ${WIND} f3=04  ${TOVPLOT} out=d04.v
	< data.A ${WIND} f3=05  ${TOVPLOT} out=d05.v
	< data.A ${WIND} f3=06  ${TOVPLOT} out=d06.v
	< data.A ${WIND} f3=07  ${TOVPLOT} out=d07.v
	< data.A ${WIND} f3=08  ${TOVPLOT} out=d08.v
	< data.A ${WIND} f3=09  ${TOVPLOT} out=d09.v
	< data.A ${WIND} f3=10  ${TOVPLOT} out=d10.v
	< data.A ${WIND} f3=11  ${TOVPLOT} out=d11.v
	vp_SideBySideAniso d00.v d01.v d02.v d03.v > r1.v
	vp_SideBySideAniso d04.v d05.v d06.v d07.v > r2.v
	vp_SideBySideAniso d08.v d09.v d10.v d11.v > r3.v
	vp_OverUnderAniso r1.v r2.v r3.v  > data.v
	< pass-mis.H Byte clip=.0007 > miss.A
	< miss.A ${WIND} f3=00  ${TOVPLOT} out=d00.v
	< miss.A ${WIND} f3=01  ${TOVPLOT} out=d01.v
	< miss.A ${WIND} f3=02  ${TOVPLOT} out=d02.v
	< miss.A ${WIND} f3=03  ${TOVPLOT} out=d03.v
	< miss.A ${WIND} f3=04  ${TOVPLOT} out=d04.v
	< miss.A ${WIND} f3=05  ${TOVPLOT} out=d05.v
	< miss.A ${WIND} f3=06  ${TOVPLOT} out=d06.v
	< miss.A ${WIND} f3=07  ${TOVPLOT} out=d07.v
	< miss.A ${WIND} f3=08  ${TOVPLOT} out=d08.v
	< miss.A ${WIND} f3=09  ${TOVPLOT} out=d09.v
	< miss.A ${WIND} f3=10  ${TOVPLOT} out=d10.v
	< miss.A ${WIND} f3=11  ${TOVPLOT} out=d11.v
	vp_SideBySideAniso d00.v d01.v d02.v d03.v > r1.v
	vp_SideBySideAniso d04.v d05.v d06.v d07.v > r2.v
	vp_SideBySideAniso d08.v d09.v d10.v d11.v > r3.v
	vp_OverUnderAniso r1.v r2.v r3.v  > miss.v
	vp_SideBySideAniso data.v miss.v > ${RESDIR}/passfill.v
	vp_Movie           data.v miss.v > ${RESDIR}/passfill.v3


sean.H: ${DATALIB}/2d_real/bp/sean.HH ${BINDIR}/Resampletime.x 
	< ${DATALIB}/2d_real/bp/sean.HH Window n3=1 f3=3 n1=500 | \
	Bandpass fhi=50 | Window j1=2 | \
	${BINDIR}/Resampletime.x d1out=.004 | Pad n1=500 > $@
	echo hff=-1 >> $@

sean2.H: sean.H
	< sean.H Window j2=2 > $@

%-pad.H %-mask2.H: %.H ${BINDIR}/LPad.x
	<$*.H ${BINDIR}/LPad.x jump=2 mask=$*-mask2.H > $*-pad.H

${RESDIR}/sean2.v: sean.H sean2.H
	< sean.H Grey out=dat.v \
	title=Original label2=Channel label1="Time (s)" >/dev/null 
	< sean2.H Grey out=mis.v \
	title=Decimated label2=Channel label1="Time (s)" >/dev/null 
	vp_SideBySideAniso dat.v mis.v > $@

sean2-pq.H: ${BINDIR}/Smooth2.x sean2.H 
	< sean2.H  ${BINDIR}/Smooth2.x nw=3 nj=2 \
	dg=0.5 eps=0.08 niter=100 > $@

sean2-p.H: ${BINDIR}/Smoothq.x sean2-pq.H sean2.H
	< sean2.H ${BINDIR}/Smoothq.x q=sean2-pq.H nw=3 nj=2 niter=100 > $@

sean2-pq1.H: ${BINDIR}/Smooth2.x sean2.H sean2-pq.H sean2-p.H
	< sean2.H ${BINDIR}/Smooth2.x nw=3 nj=2 \
	pfile=sean2-p.H qfile=sean2-pq.H \
	eps=0.08 niter=100 | Transp > $@

%-pq2.H: %-pq1.H ${BINDIR}/Double.x
	< $*-pq1.H ${BINDIR}/Double.x nw=4 | Transp > $@

%-mis.H: %-pad.H %-mask2.H %-pq2.H ${BINDIR}/Missp.x 
	< $*-pad.H ${BINDIR}/Missp.x nw=3 \
	dip=$*-pq2.H mask=$*-mask2.H niter=100 > $@

sean2-mis2.H: sean2-mis.H
	< sean2-mis.H Bandpass fhi=70 > $@


${RESDIR}/sean2-dip.v: sean2-pq2.H
	< sean2-pq2.H Window n3=1 f3=1 | Scale dscale=250 | \
	Grey out=dip1.v title="First Dip" pclip=100 wantscalebar=y \
	label2=Channel label1="Time (s)" >/dev/null 
	< sean2-pq2.H Window n3=1 f3=0 | Scale dscale=250 | \
	Grey out=dip2.v title="Second Dip" pclip=100 wantscalebar=y \
	label2=Channel label1="Time (s)" >/dev/null
	vp_SideBySideAniso dip1.v dip2.v > $@

sean2-dip.v: sean2-pq2.H
	< sean2-pq2.H Window n3=1 f3=1 | Scale dscale=250 | \
	Grey color=j out=dip1.v title="First Dip" pclip=100 wantscalebar=y \
	label2=Channel label1="Time (s)" >/dev/null 
	< sean2-pq2.H Window n3=1 f3=0 | Scale dscale=250 | \
	Grey color=j out=dip2.v title="Second Dip" pclip=100 wantscalebar=y \
	label2=Channel label1="Time (s)" >/dev/null
	vp_SideBySideAniso dip1.v dip2.v > $@

${RESDIR}/sean2-int.v: sean2-mis2.H sean.H
	Add scale=1,-1 sean2-mis2.H sean.H > sean-dif.H
	Cat sean2-mis2.H sean-dif.H | Byte > sean2-int.A
	< sean2-int.A Window n3=1 f3=0 | Ta2vplot out=int.v \
	title="Interpolation Result" \
	label2=Channel label1="Time (s)" >/dev/null 
	< sean2-int.A Window n3=1 f3=1 j2=2 f2=1 | Ta2vplot out=dif.v \
	title="Interpolation Error" \
	label2=Channel label1="Time (s)" >/dev/null 
	vp_SideBySideAniso int.v dif.v > $@

${RESDIR}/sean2-close.v ${RESDIR}/sean2-close.v3: sean.H sean2-mis2.H
	< sean.H Window min1=2.5 max1=3.1 max2=120 | \
	Grey title=Original label2=Channel label1="Time (s)" >/dev/null \
	out=orig.v
	< sean2-mis2.H Window min1=2.5 max1=3.1 max2=120 | \
	Grey title=Interpolated label2=Channel label1="Time (s)" >/dev/null \
	out=intp.v
	vp_SideBySideAniso orig.v intp.v > ${RESDIR}/sean2-close.v
	vp_Movie orig.v intp.v > ${RESDIR}/sean2-close.v3

#qdome.H: ${BINDIR}/Qdome.x qdome.p3 
#	${BINDIR}/Qdome.x par=qdome.p3 > qdome.H 

qdome.H: ${BINDIR}/Qdome-bob.x qd.p
	${BINDIR}/Qdome-bob.x >  f.H par=qd.p
	< f.H Smooth rect1=3 diff1=1  | Smooth rect1=4 | \