RESULTSER = prc02000 prc12000 prc20000 mest02000 mest12000 mest20000 pest02000 pest12000 pest20000 or00deg pest00deg orstk peststk migstk zsbty1 zsbty2 zdt2h1 zdt2h2 z00dg1 z00dg2 zpt00dg1 zpt00dg2 zstk zzo
RESULTSCR = cig02000 cig12000 cig20000 or02000 or12000 or20000
RESULTSNR =

RESDIR=./Fig
BINDIR=./Bin
DATADIR=./Data
PARDIR=./Par
B=${BINDIR}
R=${RESDIR}
D=${DATADIR}
P=${PARDIR}

UF90LIBS=-lsep2df90 -lsepparf90 -lseppar -lsupersetf90 -lsuperset -lsepgeef90
UF90FLAGS=-O1
UPSTEXFLAGS=invras=y
LATOPTS  = style=html

####################################################
#
# Binaries
#
avo:AVO.x
	mv -f AVO.x ${B}/.

ult:Ultimatch-new-all.x
	mv -f Ultimatch-new-all.x ${B}/.

####################################################
#
#    H files  
#
prim_simulated.H: ${B}/AVO.x ${P}/mutei.P
	${B}/AVO.x par=${P}/mutei.P solv=1 <${D}/prim_data_h.H >${D}/$@

primest01.H: ${B}/Ultimatch-new-all.x ${D}/prim_data_h.H ${D}/prim_simulated.H
	${B}/Ultimatch-new-all.x < ${D}/prim_data_h.H \
                multiples_in=${D}/prim_simulated.H \
                mu=1. mode=2 nsmode=2 wmode=1 nfilt=11,11 npatch=50,1 \
                op=1 niter=50 2D=1 eps=30.0 > ${D}/$@

multest01.H: ${D}/prim_data_h.H ${D}/primest01.H
	Math file1=${D}/prim_data_h.H file2=${D}/primest01.H exp="file1-file2" > ${D}/$@

multest02.H: ${B}/Ultimatch-new-all.x ${D}/prim_data_h.H ${D}/primest01.H ${D}/multest01.H
	${B}/Ultimatch-new-all.x < ${D}/prim_data_h.H \
                primaries_in=${D}/primest01.H multiples_in=${D}/multest01.H mu=1. \
                primaries_out=${D}/primest02.H \
                mode=2 nsmode=2 wmode=1 nfilt=21,11 npatch=50,11 \
                op=2 niter=60 2D=1 eps=30.0 > ${D}/$@

primest02_subt.H: ${D}/prim_data_h.H ${D}/multest02.H
	Math file1=${D}/prim_data_h.H file2=${D}/multest02.H exp="file1-file2" > ${D}/$@

stk_mig2ang.H: ${D}/mig2ang.H
	<${D}/mig2ang.H Stack > ${D}/$@

stk_prim_data_h.H: ${D}/prim_data_h.H
	<${D}/prim_data_h.H Stack > ${D}/$@

stk_primest02_subt.H: ${D}/primest02_subt.H
	<${D}/primest02_subt.H Stack > ${D}/$@

####################################################
####################################################
#############  END H FILES  ########################
####################################################

####################################################
####################################################
#############  FIGURES  ############################
####################################################
#
# CIGs
#
or02000.v: 
	<${D}/prim_data_h.H Window f3=75 n3=1 | Grey newclip=1 bclip=-0.8 eclip=0.8 title=' ' label1='Depth(m)' label2='angle(deg)' > /dev/null out=${R}/$@

or12000.v: 
	<${D}/prim_data_h.H Window f3=450 n3=1 | Grey newclip=1 bclip=-0.8 eclip=0.8 title=' ' label1='Depth(m)' label2='angle(deg)' > /dev/null out=${R}/$@

or20000.v:
	<${D}/prim_data_h.H Window f3=750 n3=1 | Grey newclip=1 bclip=-0.8 eclip=0.8 title=' ' label1='Depth(m)' label2='angle(deg)' > /dev/null out=${R}/$@

prc02000.v: 
	<${D}/prim2_prc.H Window f3=75 n3=1 | Grey newclip=1 bclip=-0.8 eclip=0.8 title=' ' label1='Depth(m)' label2='angle(deg)' > /dev/null out=${R}/$@

prc12000.v: 
	<${D}/prim_simulated.H Window f3=450 n3=1 | Grey newclip=1 bclip=-0.8 eclip=0.8 title=' ' label1='Depth(m)' label2='angle(deg)' > /dev/null out=${R}/$@

prc20000.v: 
	<${D}/prim_simulated.H Window f3=750 n3=1 | Grey newclip=1 bclip=-0.8 eclip=0.8 title=' ' label1='Depth(m)' label2='angle(deg)' > /dev/null out=${R}/$@

pest02000.v: 
	<${D}/primest02_subt.H Window f3=75 n3=1 | Grey newclip=1 bclip=-0.8 eclip=0.8 title=' ' label1='Depth(m)' label2='angle(deg)' > /dev/null out=${R}/$@

pest12000.v:
	<${D}/primest02_subt.H Window f3=450 n3=1 | Grey newclip=1 bclip=-0.8 eclip=0.8 title=' ' label1='Depth(m)' label2='angle(deg)' > /dev/null out=${R}/$@

pest20000.v: 
	<${D}/primest02_subt.H Window f3=750 n3=1 | Grey newclip=1 bclip=-0.8 eclip=0.8 title=' ' label1='Depth(m)' label2='angle(deg)' > /dev/null out=${R}/$@

mest02000.v:
	<${D}/multest02.H Window f3=75 n3=1 | Grey newclip=1 bclip=-0.8 eclip=0.8 title=' ' label1='Depth(m)' label2='angle(deg)' > /dev/null out=${R}/$@

mest12000.v:
	<${D}/multest02.H Window f3=450 n3=1 | Grey newclip=1 bclip=-0.8 eclip=0.8 title=' ' label1='Depth(m)' label2='angle(deg)' > /dev/null out=${R}/$@

mest20000.v:
	<${D}/multest02.H Window f3=750 n3=1 | Grey newclip=1 bclip=-0.8 eclip=0.8 title=' ' label1='Depth(m)' label2='angle(deg)' > /dev/null out=${R}/$@

####################################################
#
# CAGs
#
or00deg.v:
	Window f2=86 n2=1 < ${D}/prim_data_h.H | Grey newclip=1 bclip=-1. eclip=1. title=' ' label1='Depth(m)' label2='CMP_x(m)' > /dev/null out=${R}/$@

pest00deg.v:
	Window f2=86 n2=1 < ${D}/primest02_subt.H | Grey newclip=1 bclip=-1. eclip=1. title=' ' label1='Depth(m)' label2='CMP_x(m)' > /dev/null out=${R}/$@

###################################################
#
# Stacks
#
migstk.v: 
	<${D}/stk_mig2ang.H Grey newclip=1 bclip=-120. eclip=120. title=' ' label1='Depth(m)' label2='CMP_x(m)' > /dev/null out=${R}/$@

orstk.v: 
	<${D}/stk_prim_data_h.H Grey newclip=1 bclip=-120. eclip=120. title=' ' label1='Depth(m)' label2='CMP_x(m)' > /dev/null out=${R}/$@

peststk.v:
	<${D}/stk_primest02_subt.H Grey newclip=1 bclip=-120. eclip=120. title=' ' label1='Depth(m)' label2='CMP_x(m)' > /dev/null out=${R}/$@

###################################################
#
# Group CIGs
#
cig02000.v: or02000.v pest02000.v mest02000.v prc02000.v
	vp_SideBySideAniso ${R}/or02000.v ${R}/pest02000.v > 1.v
	vp_SideBySideAniso ${R}/mest02000.v ${R}/prc02000.v > 2.v
	vp_OverUnderAniso 1.v 2.v > ${R}/$@

cig12000.v: or12000.v pest12000.v mest12000.v prc12000.v
	vp_SideBySideAniso ${R}/or12000.v ${R}/pest12000.v > 1.v
	vp_SideBySideAniso ${R}/mest12000.v ${R}/prc12000.v > 2.v
	vp_OverUnderAniso 1.v 2.v > ${R}/$@

cig20000.v: or20000.v pest20000.v mest20000.v prc20000.v
	vp_SideBySideAniso ${R}/or20000.v ${R}/pest20000.v > 1.v
	vp_SideBySideAniso ${R}/mest20000.v ${R}/prc20000.v > 2.v
	vp_OverUnderAniso 1.v 2.v > ${R}/$@

###################################################
#
# Zoom stacks and 0 deg
#
zsbty1.v:
	Window min1=3500. max1=4700. max2=5500. < ${D}/stk_primest02_subt.H | Grey newclip=1 bclip=-40. eclip=40. title=' ' label2='CMP_x(m)' label1='depth(m)' > /dev/null out=${R}/$@

zdt2h1.v:
	Window min1=3500. max1=4700. max2=5500. < ${D}/stk_prim_data_h.H | Grey newclip=1 bclip=-40. eclip=40. title=' ' label2='CMP_x(m)' label1='depth(m)' > /dev/null out=${R}/$@

zsbty2.v:
	Window min1=3700. max1=4700. min2=10000. < ${D}/stk_primest02_subt.H | Grey newclip=1 bclip=-60. eclip=60. title=' ' label2='CMP_x(m)' label1='depth(m)' > /dev/null out=${R}/$@

zdt2h2.v:
	Window min1=3700. max1=4700. min2=10000. < ${D}/stk_prim_data_h.H | Grey newclip=1 bclip=-60. eclip=60. title=' ' label2='CMP_x(m)' label1='depth(m)' > /dev/null out=${R}/$@

z00dg1.v:
	Window f2=86 n2=1 min1=3500. max1=4700. max3=5500. < ${D}/prim_data_h.H | Grey newclip=1 bclip=-.8 eclip=.8 title=' ' label2='CMP_x(m)' label1='depth(m)'> /dev/null out=${R}/$@

zpt00dg1.v:
	Window f2=86 n2=1 min1=3500. max1=4700. max3=5500. < ${D}/primest02_subt.H | Grey newclip=1 bclip=-.8 eclip=.8 title=' ' label2='CMP_x(m)' label1='depth(m)'> /dev/null out=${R}/$@

z00dg2.v: 
	Window f2=86 n2=1 min1=3700. max1=4700. min3=10000. < ${D}/prim_data_h.H | Grey newclip=1 bclip=-.8 eclip=.8 title=' ' label2='CMP_x(m)' label1='depth(m)'> /dev/null out=${R}/$@

zpt00dg2.v:
	Window f2=86 n2=1 min1=3700. max1=4700. min3=10000. < ${D}/primest02_subt.H | Grey newclip=1 bclip=-.8 eclip=.8 title=' ' label2='CMP_x(m)' label1='depth(m)'> /dev/null out=${R}/$@

####################################################
#
# Group ZA sections
#
zstk.v: zdt2h1.v zsbty1.v zdt2h2.v zsbty2.v
	vp_SideBySideAniso ${R}/zdt2h1.v ${R}/zsbty1.v > 1.v
	vp_SideBySideAniso ${R}/zdt2h2.v ${R}/zsbty2.v > 2.v
	vp_OverUnderAniso 1.v 2.v > ${R}/$@

zzo.v: z00dg1.v zpt00dg1.v z00dg2.v zpt00dg2.v
	vp_SideBySideAniso ${R}/z00dg1.v ${R}/zpt00dg1.v > 1.v
	vp_SideBySideAniso ${R}/z00dg2.v ${R}/zpt00dg2.v > 2.v
	vp_OverUnderAniso 1.v 2.v > ${R}/$@

####################################################
####################################################
############ END FIGURES  ##########################
####################################################