makefile

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

sall.% : 
	@-cd miss; $(MAKE) sall.$*;

mall.% : 
	@-cd miss; $(MAKE) mall.$*;

seabdat.% : 
	@-cd seab; $(MAKE) seabdat.$*;

seabold.% : 
	@-cd seab; $(MAKE) seabold.$*;

seabnew.% : 
	@-cd seab; $(MAKE) seabnew.$*;

${RESDIR}/auto%.v: \
	${BINDIR}/Helicon.x ${BINDIR}/Autocorr.x ${BINDIR}/Wilson.x %.H %.P
	Spike n1=10 n2=10 k1=5 k2=5 d1=1 d2=1 > spik.H
	<$*.H ${BINDIR}/Autocorr.x > auto$*.H
	<spik.H ${BINDIR}/Helicon.x adj=0 div=0 filt=$*.H > filt.H
	<filt.H ${BINDIR}/Helicon.x adj=1 div=0 filt=$*.H > auto.H
	<auto$*.H ${BINDIR}/Wilson.x par=$*.P >fact$*.H
	<spik.H ${BINDIR}/Helicon.x adj=0 div=0 filt=fact$*.H > fact.H
	<filt.H ${BINDIR}/Helicon.x adj=0 div=1 filt=fact$*.H > spik2.H
	<filt.H  ${GREY} out=filt.v  title=filter
	<auto.H  ${GREY} out=auto.v  title=autocorrelation
	<fact.H  ${GREY} out=fact.v  title="Wilson factor"
	<spik2.H ${GREY} out=spik.v  title="Wilson inversion"
	vp_SideBySideAniso filt.v auto.v > top.v 
	vp_SideBySideAniso fact.v spik.v > bot.v
	vp_OverUnderAniso top.v bot.v > junk.v
	<junk.v  vppen vpstyle=n txscale=2.4 > ${RESDIR}/auto$*.v

${RESDIR}/laplac.v : \
	${BINDIR}/Helicon.x ${BINDIR}/Autocorr.x ${BINDIR}/Wilson.x \
	thinfac42.H thinplate.H laplac.H
	Spike n1=40 n2=40 k1=20 k2=20 d1=1 d2=1 > spik.H
	<spik.H ${BINDIR}/Helicon.x adj=0 div=0 filt=laplac.H > lap1.H
	<spik.H ${BINDIR}/Helicon.x adj=1 div=0 filt=laplac.H > lap2.H
	Add lap1.H lap2.H > lap.H
	<spik.H ${BINDIR}/Helicon.x adj=0 div=0 filt=thinplate.H > bih1.H
	<spik.H ${BINDIR}/Helicon.x adj=1 div=0 filt=thinplate.H > bih2.H
	Add bih1.H bih2.H > bih.H
	<spik.H ${BINDIR}/Helicon.x adj=0 div=0 filt=thinfac42.H > fact.H
	<bih.H  ${BINDIR}/Helicon.x adj=0 div=1 filt=thinfac42.H > spk1.H
	<spk1.H ${BINDIR}/Helicon.x adj=1 div=1 filt=thinfac42.H > spk2.H
	<lap.H   Window n1=10 n2=10 f1=15 f2=15 | ${GREY} out=filt.v  title=laplacian
	<bih.H   Window n1=10 n2=10 f1=15 f2=15 | ${GREY} out=auto.v  title=autocorrelation
	<fact.H  Window n1=10 n2=10 f1=15 f2=15 | ${GREY} out=fact.v  title="Wilson factor"
	<spk2.H  Window n1=10 n2=10 f1=15 f2=15 | ${GREY} out=spik.v  title="Wilson inversion"
	vp_SideBySideAniso filt.v auto.v > top.v 
	vp_SideBySideAniso fact.v spik.v > bot.v
	vp_OverUnderAniso top.v bot.v > junk.v
	<junk.v  vppen vpstyle=n txscale=2.4 > $@

%.H: %.HH
	<$*.HH Dd esize=4 | Window >$*.H
%.P: %.H
	grep lag $*.H > $*.P
	grep n1 $*.H >> $*.P

# Mathematica creates a wrong Bounding Box, and we need to fix it with sed
${RESDIR}/helix1.ps: ${RESDIR}/helix.ps
	cat ${RESDIR}/helix.warning > ${RESDIR}/helix1.warning	
	cat helix.s  > $@	
	sed -e 1,2d ${RESDIR}/helix.ps >>$@

border.V: ${RAINDIR}/border.HH
	< ${RAINDIR}/border.HH Window n1=1 f1=0 n2=338 > x.H
	< ${RAINDIR}/border.HH Window n1=1 f1=1 n2=338 > y.H
	Cmplx x.H y.H > xy.H
	< xy.H Graph ${MM} title=" " plotcol=7 plotfat=2 > border1.V
	< ${RAINDIR}/border.HH Window n1=1 f1=0 f2=338 n2=234 > x.H
	< ${RAINDIR}/border.HH Window n1=1 f1=1 f2=338 n2=234 > y.H
	Cmplx x.H y.H > xy.H
	< xy.H Graph ${MM} title=" " plotcol=7 plotfat=2 > border2.V
	< ${RAINDIR}/border.HH Window n1=1 f1=0 f2=572  > x.H
	< ${RAINDIR}/border.HH Window n1=1 f1=1 f2=572 > y.H
	Cmplx x.H y.H > xy.H
	< xy.H Graph ${MM} title=" " plotcol=7 plotfat=2 > border3.V
	Vppen < border1.V border2.V border3.V \
	erase=once vpstyle=n big=n >$@

${RESDIR}/elev.v: ${RAINDIR}/elevation.HH ${BINDIR}/Igrad.x border.V
	< ${RAINDIR}/elevation.HH ${BINDIR}/Igrad.x > elev.H
	< elev.H Grey title=Elevation transp=n yreverse=n \
	wantaxis=n wantlabel=n wheretitle=t wherexlabel=b > elev.V
	Vppen < elev.V border.V  \
	erase=once vpstyle=n big=n >/dev/null out=$@

${RESDIR}/alldata.v: ${RAINDIR}/alldata.HH
	< ${RAINDIR}/alldata.HH Window n1=1 f1=0 > x.H
	< ${RAINDIR}/alldata.HH Window n1=1 f1=1  > y.H
	Cmplx x.H y.H > xy.H
	< xy.H Graph ${MM} symbol=x title="All data locations" \
	>/dev/null out=$@

obsdata.v: ${RAINDIR}/obsdata.HH border.V
	< ${RAINDIR}/obsdata.HH Window n1=1 f1=0  > x.H
	< ${RAINDIR}/obsdata.HH Window n1=1 f1=1  > y.H
	Cmplx x.H y.H > xy.H
	< xy.H Graph ${MM} symbol=o title="Observed data locations" \
	plotcol=7 >obs.V
	Vppen < obs.V border.V \
	erase=once vpstyle=n big=n >/dev/null out=$@

data.v: ${RAINDIR}/obsdata.HH ${RAINDIR}/coord.HH border.V
	< ${RAINDIR}/obsdata.HH Window n1=1 f1=0 > x.H
	< ${RAINDIR}/obsdata.HH Window n1=1 f1=1 > y.H
	Cmplx x.H y.H > xy.H
	< xy.H Graph ${MM} symbol=o title="Data locations" \
	plotcol=7 >obs.V
	< ${RAINDIR}/coord.HH Window n1=1 f1=0 > x.H
	< ${RAINDIR}/coord.HH Window n1=1 f1=1 > y.H
	Cmplx x.H y.H > xy.H
	< xy.H Graph ${MM} symbol=x title=" " >coord.V
	Vppen < obs.V coord.V border.V  \
	erase=once vpstyle=n big=n >/dev/null out=$@

${RESDIR}/raindata.v: obsdata.v data.v
	vp_SideBySideIso data.v obsdata.v > $@

${RESDIR}/wilson.txt: ${BINDIR}/Wilsontest.x
	${BINDIR}/Wilsontest.x > $@

lapinter%.H: ${RAINDIR}/elevation.HH ${RAINDIR}/obsdata.HH ${BINDIR}/Lapinter.x
	< ${RAINDIR}/obsdata.HH Window n1=2 | Transp > xy.H
	< ${RAINDIR}/obsdata.HH Window f1=2 > raindat.H
	time < ${RAINDIR}/elevation.HH \
	${BINDIR}/Lapinter.x coord=xy.H dat=raindat.H niter=$* > $@

lapinter%.v: lapinter%.H
	< lapinter$*.H Grey wantscalebar=y yreverse=n transp=n allpos=y \
	color=I pclip=100 title="Laplace $*" >/dev/null out=grey.v
	< lapinter$*.H Contour nc=4 c0=50 dc=100 yreverse=n transp=n \
	wanttitle=n plotcol=7 wantframe=y wantaxis=n \
	screenwd=12.17 >/dev/null out=cont.v
	vp_Overlay grey.v cont.v > $@

#lapinter%.v: lapinter%.H
#	< lapinter$*.H Grey wantscalebar=y yreverse=n transp=n \
#	color=j pclip=100 wheretitle=t wherexlabel=b \
#	title="Model-space: $* iterations" >/dev/null out=$@

${RESDIR}/lapinter.v: lapinter10.v lapinter100.v lapinter1000.v lapinter10000.v
	vppen < lapinter10.v lapinter100.v lapinter1000.v lapinter10000.v \
	vpstyle=n size=r gridnum=2,2 >$@

lapdat%.H: lapinter%.H ${RAINDIR}/coord.HH ${BINDIR}/Extract.x
	< ${RAINDIR}/coord.HH Window | Transp > xy.H
	< lapinter$*.H ${BINDIR}/Extract.x coord=xy.H > lapdat$*.H

${RESDIR}/lapstat.v: lapdat10000.H ${RAINDIR}/predict.HH 
	Cmplx ${RAINDIR}/predict.HH lapdat10000.H > xy.H
	Graph < xy.H symbol=* title="Laplace Correlation" >/dev/null out=$@

tridat.H: tridata.H ${RAINDIR}/coord.HH ${BINDIR}/Extract.x
	< ${RAINDIR}/coord.HH Window | Transp > xy.H
	< tridata.H ${BINDIR}/Extract.x coord=xy.H > $@

${RESDIR}/tristat.v: tridat.H ${RAINDIR}/predict.HH 
	Cmplx ${RAINDIR}/predict.HH tridat.H > xy.H
	Graph < xy.H symbol=* title="Triang Correlation" >/dev/null out=$@

lapstat.v: lapdat10000.H ${RAINDIR}/predict.HH 
	Cmplx ${RAINDIR}/predict.HH lapdat10000.H > xy.H
	Graph < xy.H symbol=* title="Correlation" >/dev/null out=$@

rand500.H: ${BINDIR}/Random.x
	${BINDIR}/Random.x tmax=80. n1=500 nd=50 > $@

rand50.H: ${BINDIR}/Random.x
	${BINDIR}/Random.x tmax=80. n1=50 nd=500 > $@

splbin%.H: rand%.H ${BINDIR}/Bin.x
	< rand$*.H ${BINDIR}/Bin.x n1=$* > $@

${RESDIR}/bin%.v: splbin%.H
	< splbin$*.H ${DOTS} out=$@

envint%.H emm%.H: rand%.H ${BINDIR}/Invinter.x ${BINDIR}/Spline1.x
	${BINDIR}/Spline1.x tension=0.01 > flt.H
	time < rand$*.H ${BINDIR}/Invinter.x niter=20 \
	n1=$* filt=flt.H mm=emm$*.H > envint$*.H

snvint%.H smm%.H: rand%.H ${BINDIR}/Invinter.x \
	${BINDIR}/Spline3.x ${BINDIR}/Splinefilter.x
	${BINDIR}/Spline3.x tension=0.01 > flt.H
	time < rand$*.H ${BINDIR}/Invinter.x niter=20 \
	n1=$* filt=flt.H mm=mm.H spline=y nw=4 > invint.H
	< mm.H ${BINDIR}/Splinefilter.x post=y nw=4 > smm$*.H
	< invint.H ${BINDIR}/Splinefilter.x post=y nw=4 > snvint$*.H

alias%.H: ${BINDIR}/Alias2.x
	${BINDIR}/Alias2.x tmax=80. n1=$* > $@

enorm%.H: emm%.H alias%.H ${BINDIR}/Modnorm1.x
	 < emm$*.H ${BINDIR}/Modnorm1.x ideal=alias$*.H > $@

snorm%.H: smm%.H alias%.H ${BINDIR}/Modnorm1.x
	 < smm$*.H ${BINDIR}/Modnorm1.x ideal=alias$*.H > $@

${RESDIR}/norm%.v: enorm%.H snorm%.H
	Cat axis=2 enorm$*.H snorm$*.H | \
	${GRAPH2} title="Model convergence $*" \
	legend=1 legendloc="tr" legendtitle=" " legendsz=10 \
	label1=Iterations label2="Model misfit" \
	curvelabel="Linear:B-Spline" out=$@

${RESDIR}/waves.v : ${BINDIR}/Helicon.x waves.H 
	Spike n1=20 n2=20 nsp=2 k1=3,7 k2=7,3 mag=1,-1  > spike.H
	Spike n1=20 n2=20 nsp=2 k1=19,15 k2=13,6 mag=1,-1 > inp2.H
	<spike.H ${BINDIR}/Helicon.x filt=waves.H > inp1.H
	Add inp1.H inp2.H > inp.H
	<inp.H ${BINDIR}/Helicon.x   filt=waves.H div=1 > div.H
	<inp.H ${GREY} out=inp.v  title=input
	<div.H ${GREY} out=div.v  title=input/filter
	vp_SideBySideAniso inp.v div.v > junk.v
	<junk.v  vppen vpstyle=n txscale=1.4  > $@

thinfac%.H: thinplate.H ${BINDIR}/Wilson.x fact%.p
	< thinplate.H ${BINDIR}/Wilson.x par=fact$*.p > $@

${RESDIR}/thin%.v : ${BINDIR}/Helicon.x thinfac%.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 ${BINDIR}/Helicon.x filt=thinfac$*.H > inp1.H
	Add inp1.H inp2.H > inp.H
	<inp.H ${BINDIR}/Helicon.x       div=1 filt=thinfac$*.H > div.H
	<div.H ${BINDIR}/Helicon.x adj=1 div=1 filt=thinfac$*.H > div2.H
	<inp.H  ${GREY} out=inp.v  clip=1 title=input
	<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 > junk.v
	<junk.v  vppen vpstyle=n txscale=2.4 > $@

precinter%.H: ${RAINDIR}/elevation.HH thinfac42.H \
	      ${RAINDIR}/obsdata.HH ${BINDIR}/Pinter1.x
	< ${RAINDIR}/obsdata.HH Window n1=2 | Transp  > xy.H
	< ${RAINDIR}/obsdata.HH Window f1=2 > dat.H
	time < ${RAINDIR}/elevation.HH ${BINDIR}/Pinter1.x \
	niter=$* coord=xy.H dat=dat.H filt=thinfac42.H > $@

TRIFILES = list_struct delaunay heap predicates
TRIINCLUDS = ${addprefix ${SRCDIR}/, ${addsuffix .h, $(TRIFILES)}}
TRISOURCES = ${addprefix ${SRCDIR}/, ${addsuffix .c, $(TRIFILES)}}
TRIOBJECTS = ${addprefix ${OBJDIR}/, ${addsuffix .o, $(TRIFILES)}}

SEPLIB = -L$(SEPROOT)/lib -lsep
SEPROOT = /usr/local/SEP

# without graphics
${BINDIR}/Sept2a.x: ${SRCDIR}/sept2a.c $(TRISOURCES)
	$(MAKE) CFLAGS="-DWITH_TIME" ${OBJDIR}/sept2a.o $(TRIOBJECTS) 
	$(CC) -o $@ ${OBJDIR}/sept2a.o $(TRIOBJECTS) $(SEPLIB) -lm

# seplib i/o
${OBJDIR}/sept2a.o: ${SRCDIR}/sept2a.c 
	$(CC) -c $(CFLAGS) -I$(SRCDIR) -I$(SEPROOT)/include \
	-DSOURCE='"'$(SRCDIR)/sept2a.c'"' $(SRCDIR)/sept2a.c
	${MV} sept2a.o ${OBJDIR}

tridata.H edge.H: ${BINDIR}/Sept2a.x ${RAINDIR}/obsdata.HH
	< ${RAINDIR}/obsdata.HH Window | Transp > xyz.H
	< xyz.H ${BINDIR}/Sept2a.x n1=376 n2=253 \
	d1=1.009975 d2=1.009975 o1=-185.556375 o2=-127.261523 \
	> tridata.H edgeout=edge.H

${RESDIR}/trinter.v: tridata.H edge.H
	<edge.H Graph plotcol=1 title=" " wantaxis=n wantlabel=n \
	${MM} > edge.V
	< tridata.H Grey title=Triangulation transp=n yreverse=n \
	wheretitle=t wherexlabel=b color=I > data.V
	Vppen < data.V edge.V  \
	erase=once vpstyle=n big=n >/dev/null out=$@

precinter%.V: precinter%.H
	< precinter$*.H Grey wantscalebar=y yreverse=n transp=n \
	color=I pclip=100 title="Preconditioning $*" >$@

precinter%.v: precinter%.H
	< precinter$*.H Grey wantscalebar=y yreverse=n transp=n allpos=y \
	color=I pclip=100 title="Preconditioning $*" >/dev/null out=grey.v
	< precinter$*.H Contour nc=4 c0=50 dc=100 yreverse=n transp=n \
	wanttitle=n plotcol=7 wantframe=y wantaxis=n \
	screenwd=12.17 >/dev/null out=cont.v
	vp_Overlay grey.v cont.v > $@

#
#precinter%.v: precinter%.H
#	< precinter$*.H Grey wantscalebar=y yreverse=n transp=n \
#	color=j pclip=100 wheretitle=t wherexlabel=b \
#	title="Preconditioning: $* iterations" >/dev/null out=$@

${RESDIR}/precinter.v: precinter10.v precinter100.v
	vppen < precinter10.v precinter100.v \
	vpstyle=n size=r gridnum=2,1 >$@

clean: jclean
	@-${TOUCH} junk.P
	@-${RM}       *.P
	@-cd miss; $(MAKE) clean;
	@-cd seab; $(MAKE) clean;
	@-${RM} ${BINDIR}/*.x
	@-${RM} ${OBJDIR}/*.o

include ${SEPINC}/SEP.bottom