makefile

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

include ${SEPINC}/SEP.top

UTUBEFLAGS = see_progress=y
UF90LIBS = ${GEELIB}

RESULTSER = data im1 fm1 fm6 pm1 early1 schwab1 \
	mall sall seabdat seabold seabnew \
	waves autowaves thin42 laplac \
	bin500 bin50 norm500 norm50 \
	elev raindata lapinter lapstat precinter
RESULTSCR = helix helix1


RESDIR = ./Fig
DATDIR = ./Dat

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

RAINDIR = ${DATALIB}/fun/chernobyl

MM = min1=-185.556 max1=193.18275 min2=-127.262 max2=127.25044 

LATOPTS = option=minimal

GRAPH = Graph clip=1 max2=1 min2=-1 crowd=.99 > junk.H
GRAPH2 = Graph >/dev/null dash=1,0
DOTS = Dots >/dev/null dots=2 connect=0 strings=1 labelsz=14
DOTY = Dots dots=0 connect=0 gaineach=0 constsep=1 strings=0 clip=1.2 > /dev/null
LBL = label1="Iterations"
TAVPLOT = Ta2vplot crowd=.85  wantaxis=0 >/dev/null
MPAD = prec=0 padin=2 padout=2
DOTARGS2 = dots=2 connect=0 gaineach=0 overlap=1.30 labelsz=16
DOTARGS = dots=0 strings=0 connect=1 constsep=1  \
	labels="model-space:data-space:de-spike:input" \
	labelsz=6 Xclip=2.5 gaineach=0 overlap=.95
PLOTSIG = Dots ${DOTARGS} > /dev/null
CUBE = Cubeplot flat=n  \
	point1=0.9 point2=0.9  frame1=100 frame2=200 frame3=8 >/dev/null
CUBE1 = Cubeplot flat=n  \
	point1=0.9 point2=0.9  frame1=100 frame2=280 frame3=11 >/dev/null

COLOR=y

GREY = Grey gainpanel=every pclip=100 wantaxis=n >/dev/null \
	crowd=.88 gpow=.7 label1=' ' label2=' '

.PRECIOUS: alias.P data.H head.H

default : view

${RESDIR}/data.v alias.P dat.v: ${BINDIR}/Alias.x ${BINDIR}/Ideal.x
	 ${BINDIR}/Alias.x tmax=80. pairs=1   n1=200 > alias.P
	 ${BINDIR}/Alias.x tmax=80. pairs=0   n1=200 > alias.H
	<alias.H  ${DOTY}  strings=1 out=dat.v
	${BINDIR}/Ideal.x tmax=80. pairs=0 n1=200 > ideal.H
	< ideal.H ${DOTY} connect=1 out=ideal.v
	vp_Overlay dat.v ideal.v > both.v
	vp_SideBySideAniso both.v dat.v > ${RESDIR}/data.v

aliasbin.H: alias.P ${BINDIR}/Bin.x
	< alias.P ${BINDIR}/Bin.x n1=200 > $@

aliasbin.P: aliasbin.H ${BINDIR}/Aliasreg.x
	<  aliasbin.H ${BINDIR}/Aliasreg.x > $@

aliasbin.v:  aliasbin.H  mlines.H ${BINDIR}/Miss.x
	<aliasbin.H ${BINDIR}/Miss.x filt=mlines.H niter=200 prec=1 >int.H
	< int.H ${DOTY}  connect=1 out=$@

${RESDIR}/bm$*.v ${RESDIR}/bm%.v3 bm%.H: aliasbin.H mlines.H \
	${BINDIR}/Miss.x  dat.v
	<aliasbin.H ${BINDIR}/Miss.x mmov=bm$*.H \
	filt=mlines.H niter=301 prec=1 >int.H
	<bm$*.H Transp plane=23     > transp.H
	<transp.H ${GRAPH} out=${RESDIR}/bm$*.v3 fastplot=20
	<int.H ${DOTY}  connect=1 out=int.v
	vp_Overlay dat.v int.v > ${RESDIR}/bm$*.v

${RESDIR}/im%.v ${RESDIR}/im%.v3 im%.H: \
	dat.v alias.P  ${BINDIR}/Invint.x  
	<alias.P ${BINDIR}/Invint.x n1=200 mm=im$*.H > int.H method=1 filter=$*
	<im$*.H Transp plane=23     > transp.H
	<transp.H ${GRAPH} out=${RESDIR}/im$*.v3 fastplot=20
	<int.H ${DOTY}  connect=1 out=int.v
	vp_Overlay dat.v int.v > ${RESDIR}/im$*.v

${RESDIR}/fm%.v ${RESDIR}/fm%.v3 fm%.H: \
	dat.v  alias.P ${BINDIR}/Invint.x  
	<alias.P ${BINDIR}/Invint.x n1=200 mm=fm$*.H > int.H \
	method=2 filter=$* eps=.01
	<fm$*.H Transp plane=23     > transp.H
	<transp.H ${GRAPH} out=${RESDIR}/fm$*.v3 fastplot=20
	<int.H ${DOTY}  connect=1 out=int.v
	vp_Overlay dat.v int.v > ${RESDIR}/fm$*.v

${RESDIR}/pm%.v ${RESDIR}/pm%.v3 pm%.H: \
	dat.v  alias.P ${BINDIR}/Invint.x  
	<alias.P ${BINDIR}/Invint.x n1=200 mm=pm$*.H > int.H \
	method=2 filter=$* pef=1 eps=.01
	<pm$*.H Transp plane=23     > transp.H
	<transp.H ${GRAPH} out=${RESDIR}/pm$*.v3 fastplot=20
	<int.H ${DOTY}  connect=1 out=int.v
	vp_Overlay dat.v int.v > ${RESDIR}/pm$*.v

NITER = 250
ITERATION = 5
${RESDIR}/early%.v : im%.H fm%.H    
	< im$*.H Window n2=1 f2=${ITERATION}   > iwind.H
	< fm$*.H Window n2=1 f2=${ITERATION}   > fwind.H
	< fm$*.H Window n2=1 f2=${NITER}       > icorr.H
	Merge <iwind.H fwind.H         icorr.H space=no axis=2 > junk.H
	< junk.H ${DOTY} connect=1 overlap=2.0 title=" " \
	labels=regul.${ITERATION}:precon.${ITERATION}:niter=${NITER} \
	labelsz=10 out=${RESDIR}/early$*.v

${RESDIR}/schwab%.v :  ${BINDIR}/Modnorm.x  im%.H fm%.H 
	<im$*.H ${BINDIR}/Modnorm.x > junk1.H
	<fm$*.H ${BINDIR}/Modnorm.x > junk2.H
	Merge <junk1.H junk2.H         space=n axis=1 > junk.H
	< junk.H Transp plane=12  > transp.H
	<transp.H Window f1=0 n1=100 > wind.H
	< wind.H Graph ${LBL} symbol='md' symbolsz=10 symbolsz=7 \
	title="Model Residual Vector Length" \
	min2=0 min1=-1 out=${RESDIR}/schwab$*.v >/dev/null

${RESDIR}/beam%.v :  ${BINDIR}/Modnorm.x  bm%.H fm%.H 
	<bm$*.H ${BINDIR}/Modnorm.x > junk1.H
	<fm$*.H ${BINDIR}/Modnorm.x > junk2.H
	Merge <junk1.H junk2.H         space=n axis=1 > junk.H
	< junk.H Transp plane=12  > transp.H
	<transp.H Window f1=0 n1=40 > wind.H
	< wind.H Graph ${LBL} symbol='bd' symbolsz=10 symbolsz=7 \
	title="Model Residual Vector Length" \
	min2=0 min1=-1 out=${RESDIR}/beam$*.v >/dev/null

#-------------- 3D ------------------------------------

head.H:  ${DATDIR}/head.H
	Window < ${DATDIR}/head.H n2=2 > head.H

dat.H:  ${DATDIR}/dat.H
	Window < ${DATDIR}/dat.H n2=150 f2=50 > dat.H

offset.H: ${DATDIR}/head.H
	Window < ${DATDIR}/head.H n2=2 f2=2 > $@

aoff.H: ${DATDIR}/head.H
	Window < ${DATDIR}/head.H n2=1 f2=2 > $@

head1.H: head.H aoff.H ${BINDIR}/Select.x
	< head.H ${BINDIR}/Select.x header=aoff.H hmin=500 hmax=600 > $@

dat1.H: dat.H aoff.H ${BINDIR}/Select.x
	< dat.H ${BINDIR}/Select.x header=aoff.H hmin=500 hmax=600 > $@

${RESDIR}/offset.v:  offset.H
	Window <   offset.H n2=1 f2=0 > offx.H
	Window <   offset.H n2=1 f2=1 > offy.H
	Cmplx offx.H offy.H > off.H 
	Graph < off.H  symbol='+' symbolsz=2 out=${RESDIR}/offset.v \
	label1="Offset" label2="Azimuth" title="Offset Geometry"  > /dev/null

${RESDIR}/cmp.v:  head.H
	Window <   head.H n2=1 f2=0 > cmpx.H
	Window <   head.H n2=1 f2=1 > cmpy.H
	Cmplx cmpx.H cmpy.H > cmp.H 
	Graph < cmp.H  symbol='+' symbolsz=2 out=${RESDIR}/cmp.v \
	label1="In-line midpoint" label2="Cross-line midpoint" title="CMP Geometry"  > /dev/null

${RESDIR}/cmp1.v:  head1.H
	Window <   head1.H n2=1 f2=0 > cmpx.H
	Window <   head1.H n2=1 f2=1 > cmpy.H
	Cmplx cmpx.H cmpy.H > cmp.H 
	Graph < cmp.H  symbol='+' symbolsz=2 out=${RESDIR}/cmp1.v \
	label1="In-line midpoint" label2="Cross-line midpoint" title="CMP Geometry"  > /dev/null

bin.H fold.H: grid.p dat.H head.H ${BINDIR}/Bin2.x
	${BINDIR}/Bin2.x \
	par=grid.p hff=head.H fold=fold.H interp=1 <dat.H > junk.H
	<junk.H Transp plane=31 > bin.H

bin-1.H fold-1.H: grid1.p dat1.H head1.H ${BINDIR}/Bin2.x
	${BINDIR}/Bin2.x \
	par=grid1.p hff=head1.H fold=fold-1.H interp=1 <dat1.H > junk.H
	<junk.H Transp plane=31 > bin-1.H

bin%.H: grid.p dat.H head.H ${BINDIR}/Bin2.x
	${BINDIR}/Bin2.x \
	par=grid.p hff=head.H interp=$* inv=1 <dat.H > bin$*.H

lin.H lfold.H: grid.p dat.H head.H ${BINDIR}/Bin2.x
	${BINDIR}/Bin2.x \
	par=grid.p hff=head.H fold=lfold.H interp=2 <dat.H > junk.H
	<junk.H Transp plane=31 > lin.H

elin.H: grid.p dat.H head.H ${BINDIR}/Bin2.x
	${BINDIR}/Bin2.x par=grid.p hff=head.H interp=2 inv=1 <dat.H > elin.H

inv%.H: grid.p dat.H head.H ${BINDIR}/Inv2.x
	${BINDIR}/Inv2.x \
	par=grid.p hff=head.H niter=5 interp=$* inv=2 <dat.H > junk.H
	<junk.H Transp plane=31 > inv$*.H

in1-%.H: grid1.p dat1.H head1.H ${BINDIR}/Inv2.x
	${BINDIR}/Inv2.x \
	par=grid1.p hff=head1.H niter=5 interp=$* inv=2 <dat1.H > junk.H
	<junk.H Transp plane=31 > in1-$*.H

in2-%.H: grid1.p dat1.H head1.H ${BINDIR}/Inv2.x
	${BINDIR}/Inv2.x \
	par=grid1.p hff=head1.H niter=5 interp=$* inv=3 <dat1.H > junk.H
	<junk.H Transp plane=31 > in2-$*.H

vbin.H: grid.p dat.H head.H ${BINDIR}/Inv2.x
	${BINDIR}/Inv2.x \
	par=grid.p hff=head.H niter=5 interp=1 inv=1 <dat.H > junk.H
	<junk.H Transp plane=31 > vbin.H

vbin-1.H: grid1.p dat1.H head1.H ${BINDIR}/Inv2.x
	${BINDIR}/Inv2.x \
	par=grid1.p hff=head1.H niter=5 interp=1 inv=1 <dat1.H > junk.H
	<junk.H Transp plane=31 > vbin-1.H

mres.H: grid.p dat.H head.H ${BINDIR}/Inv2.x
	${BINDIR}/Inv2.x \
	par=grid.p hff=head.H niter=5 interp=1 inv=0 <dat.H > mres.H

${RESDIR}/fold.v: fold.H
	<fold.H Taplot allpos=y pclip=100 > junk.A
	<junk.A Ta2vplot transp=n yreverse=n wantscalebar=1 \
	labelsz=10 title="Fold Map" out=${RESDIR}/fold.v > /dev/null

${RESDIR}/fold1.v: fold-1.H
	<fold-1.H Taplot allpos=y pclip=100 > junk.A
	<junk.A Ta2vplot transp=n yreverse=n wantscalebar=1 \
	labelsz=10 title="Fold Map" out=${RESDIR}/fold1.v > /dev/null

${RESDIR}/lfold.v: lfold.H
	<lfold.H Taplot allpos=y pclip=100 > junk.A
	<junk.A Ta2vplot transp=n yreverse=n wantscalebar=1 \
	labelsz=10 title="Linear Fold Map" out=${RESDIR}/lfold.v > /dev/null

${RESDIR}/hist.v: fold.H ${BINDIR}/Histogram.x
	<fold.H ${BINDIR}/Histogram.x > hist.H
	<hist.H Graph title="Fold Distribution" >/dev/null out=${RESDIR}/hist.v

${RESDIR}/hist1.v: fold-1.H ${BINDIR}/Histogram.x
	<fold-1.H ${BINDIR}/Histogram.x > hist.H
	<hist.H Graph title="Fold Distribution" >/dev/null \
	out=${RESDIR}/hist1.v

${RESDIR}/lhist.v: lfold.H ${BINDIR}/Histogram.x
	<lfold.H ${BINDIR}/Histogram.x > lhist$*.H
	<lhist.H Graph title="Linear Fold Distribution" >/dev/null out=${RESDIR}/lhist.v

${RESDIR}/bin.v: bin.H
	<bin.H Byte > bin.A
	<bin.A ${CUBE} title="Nearest Neighbor" out=${RESDIR}/bin.v

${RESDIR}/bin1.v: bin-1.H
	<bin-1.H Byte > bin.A
	<bin.A ${CUBE1} title=" " out=${RESDIR}/bin1.v

${RESDIR}/regbin.v: inv1.H
	<inv1.H Byte > inv1.A
	<inv1.A ${CUBE} title="Nearest Neighbor: Regularized Optimization" out=${RESDIR}/regbin.v

${RESDIR}/regbin1.v: in1-1.H
	<in1-1.H Byte > inv1.A
	<inv1.A ${CUBE1} title="Nearest Neighbor: Regularized Optimization" out=${RESDIR}/regbin1.v

${RESDIR}/invbin.v: vbin.H
	<vbin.H Byte > vbin.A
	<vbin.A ${CUBE} title="Nearest Neighbor: Optimization" out=${RESDIR}/invbin.v

${RESDIR}/invbin1.v: vbin-1.H
	<vbin-1.H Byte > vbin.A
	<vbin.A ${CUBE1} title=" " out=${RESDIR}/invbin1.v

${RESDIR}/lin.v: lin.H
	<lin.H Byte > lin.A
	<lin.A ${CUBE} title="Linear" out=${RESDIR}/lin$*.v

${RESDIR}/elin.v: elin.H
	<elin.H Transp plane=31 > junk.H
	<junk.H Byte > elin.A
	<elin.A ${CUBE} title="Linear: Data normalized" out=${RESDIR}/elin.v

${RESDIR}/reglin.v: inv2.H
	<inv2.H Byte > inv2.A
	<inv2.A ${CUBE} title="Linear: Regularized Optimization" out=${RESDIR}/reglin.v

${RESDIR}/reglin1.v: in1-2.H
	<in1-2.H Byte > inv2.A
	<inv2.A ${CUBE} title="Linear: Regularized Optimization" out=${RESDIR}/reglin1.v

${RESDIR}/regpef1.v: in2-2.H
	<in2-2.H Byte > inv3.A
	<inv3.A ${CUBE} title=" " out=${RESDIR}/regpef1.v

${RESDIR}/regmin.v: inv-4.H
	<inv-4.H Byte > inv-4.A
	<inv-4.A ${CUBE} title="Muir: Regularized Optimization" out=${RESDIR}/regmin.v

${RESDIR}/regtin.v: inv4.H
	<inv4.H Byte > inv4.A
	<inv4.A ${CUBE} title="Lagrange: Regularized Optimization" out=${RESDIR}/regtin.v