include ${SEPINC}/SEP.top

BINDIR = ./Bin
RESDIR = ./Fig
SRCDIR = ./Src
FFTWLIBS = -L/opt/FFTW/lib -lfftw3f
UF90LIBS = -lsupersetf90 -lsepgeef90 -lsuperset -lsepauxf90 -lsepmathf90 -lsep2df90 -lsep3df90 -lsep3d -lsepf90 -lsep -openmp -lsvml ${FFTWLIBS}
FCOLEVEL = 3
# for pompei
UF90FLAGS = -O3 -axP -openmp -Bstatic
# for kook
#UF90FLAGS = -O3 -axW -openmp -static-libcxa -Bstatic

TXTDIR = ./Txt

BPGOM3D_DATA = /net/server/tang/Research/bpgom3d
BPGOM3D_IMAG = /net/server/tang/Research/bpgom3d/save

BPGOM_DATA3D = /data/bpgom/bpgom_data.H
BPGOM_VELO3D = /data/bpgom/bpgom_velo.H

LSM3D_IMAG = /net/server/tang/Research/bpgom3d/lsm

COLORBAR = newclip=1 bclip=4.9 eclip=14.9 wantscalebar=y barlabel="Velocity (km/s)" color=j
LABELFONT = labelrot=n labelfat=2 labelsz=8
LABELAXES = wantaxes=n label1="Depth" label2="Inline" label3="Crossline" d1num=0 d2num=0 d3num=0


COLOR = no



bpgom.cmp.3d.H:
	Window3d < ${BPGOM_DATA3D} n4=396 > $@

bpgom-cmp-3d-time-line%.H: bpgom.cmp.3d.H
	Window3d < $< n3=1 f3=$* > junk1-$*.H
	Transp < junk1-$*.H plane=13 maxsize=2000 > junk2-$*.H
	echo n5=1 d5=1. o5=0. >> junk2-$*.H
	${BINDIR}/fftw_transf.x < junk2-$*.H adj=n f1=2.083333 f2=4 f3=33 f4=35 dt=0.004 radian=n sign=1 datapath=${DATAPATH} > $@

bpgom-cmp-3d-time-fix-line%.H: bpgom-cmp-3d-time-line%.H
	${BINDIR}/preprocessing.x < $< > $@  datapath=${DATAPATH}

bpgom-shot-3d-time-line%.H: bpgom-cmp-3d-time-fix-line%.H
	Create3d < $< keyname1=offset keyname2=cmp >junk1-$*.H
	Headermath <junk1-$*.H verb=1 key1=sx type1=scalar_float eqn1=cmp-offset*0.5 key2=offx type2=scalar_float eqn2=offset      >junk2-$*.H
	Headermath <junk1-$*.H verb=1 key1=sx type1=scalar_float eqn1=cmp+offset*0.5 key2=offx type2=scalar_float eqn2=-1.0*offset >junk3-$*.H
	Sort3d <junk2-$*.H verb=1 nkeys=2 key1=offx og1=0.        dg1=164.04 ng1=189 key2=sx og2=32726 dg2=82.02 ng2=801 >junk4-$*.H
	Sort3d <junk3-$*.H verb=1 nkeys=2 key1=offx og1=-30839.52 dg1=164.04 ng1=188 key2=sx og2=32726 dg2=82.02 ng2=801 >junk5-$*.H
	Infill3d < junk4-$*.H verb=1 > junk6-$*.H
	Infill3d < junk5-$*.H verb=1 > junk7-$*.H
	Cat junk7-$*.H junk6-$*.H axis=2 datapath=${DATAPATH} > $@

bpgom-shot-3d-time-interp-line%.H: bpgom-shot-3d-time-line%.H
	Interp < $< d2out=82.02 maxsize=1500 datapath=${DATAPATH} > $@

bpgom-shot-3d-freq-line%.H: bpgom-shot-3d-time-interp-line%.H
	echo n5=1 o5=0 d5=1 >> $<
	Transp < $< reshape=1,2,4,5 plane=34 maxsize=2000 > transp-$*.H
	${BINDIR}/fftw_transf.x < transp-$*.H adj=y ${BPGOM3D_FREQ} > freq-$*.H
	Transp < freq-$*.H reshape=1,3,5 maxsize=2000 datapath=${DATAPATH} > $@
	echo label1="offx" label2="offy" label3="frequency" label4="sx" label5="sy" >> $@
	Rm transp-$*.H freq-$*.H

bpgom-shot-3d-freq-all.H:
	python update_history.py --history=$@ --command=bpgom-shot-3d-freq-line.H --index_beg=1 --index=125

bpgom-shot-3d-freq-all2.H: bpgom-shot-3d-freq-all.H
	Cp $< $@

bpgom-shot-3d-freq-subsample.H: bpgom3d-shot-3d-freq-all2.H
	Window3d < $< j1=2 j4=2 j5=2 squeeze=n > $@

bpgom3d-csou-planes-subsample-line%.H: bpgom-shot-3d-freq-subsample.H
	Window3d < $< n5=1 f5=$* squeeze=n > junk0.$@
	${BINDIR}/hesdct_efnc.x nws=205 ows=18.8496 dws=0.523599 nsx=401 dsx=164.04 osx=32726 domain="encsou" enctyp="planes" efns=junk1.$@ \
        psx_min=-0.000086073 psx_max=0.000086073 npsx=101 sx0_pos=32726 sx0_neg=98342
	${BINDIR}/encode.x csou=bpgom-csou-3d.H crec=junk0.$@ efns=junk1.$@ wfds=$@ wfdr=bpgom3d-crec-planes-subsample-line$*.H offset=y wgtsou=y wgtrec=y datapath=${DATAPATH}
	Rm junk0.$@ junk1.$@


bpgom3d-bvel.H: ${BINDIR}/bpgom3d_vel_edit.x ${BPGOM_VELO3D}
	${BINDIR}/bpgom3d_vel_edit.x < ${BPGOM_VELO3D} > junk1.$@
	Interp < junk1.$@ o2out=-8202 n2out=2001 type=0 maxsize=5000 > junk2.$@
	Transp < junk2.$@ reshape=1,3 maxsize=2000 > junk3.$@
	Interp < junk3.$@ o2out=41010 n2out=851 type=0 maxsize=2000 > $@
	Rm junk1.$@ junk2.$@ junk3.$@

bpgom3d-bvel-subsample.H: bpgom3d-bvel.H
	Window3d < $< j1=2 j2=2 j3=2 squeeze=n > $@


##################################################
# computing both hx and hy with a smaller target
##################################################
bpgom3d-bimg-planes-subsample-line%.H:
	Window3d < ${DATDIR}/bpgom3d-csou-planes-subsample-all2.H n5=1 f5=$* squeeze=n min1=35000 max1=100000 | Transp plane=25 > ${WOKDIR}/bpgom3d-csou-planes-subsample-line$*.H
	Window3d < ${DATDIR}/bpgom3d-crec-planes-subsample-all2.H n5=1 f5=$* squeeze=n min1=35000 max1=100000 | Transp plane=25 > ${WOKDIR}/bpgom3d-crec-planes-subsample-line$*.H
	python bpgom3d_imag2.py --bindir=${BINDIR} --datdir=${WOKDIR} --resdir=${RESDIR} --tmpdir=${TMPDIR}\
         --csou=bpgom3d-csou-planes-subsample-line$*.H --crec=bpgom3d-crec-planes-subsample-line$*.H --bvel=bpgom3d-bvel-subsample.H --bimg=$@ \
         --jsx_win=2 --nsx_tot=101 \
         --offset=n --comsou=y --oper=ssfpi --nvrf=3 --nws=16 --nxpad_beg=0 --nxpad_end=0 --nypad_beg=100 --nypad_end=100 --ntpx=20 --ntpy=50 \
         --nhx=21 --ohx=-1640.4 --dhx=164.04 --nhy=7 --ohy=-1476.36 --dhy=492.12\
         --image_xmin=50000 --image_xmax=90000 --image_ymin=65616 --image_ymax=82000 --image_zmin=25000 --image_zmax=35000 \
         --memchk=n --queue=default --grepinfo=line$* --datapath=/tmp/


##################################################
# computing both hx and hy with a bigger target
################################################## 
bpgom3d-bimg-planes-subsample-big-line%.H: 
	Window3d < ${DATDIR}/bpgom3d-csou-planes-subsample-all2.H n5=1 f5=$* squeeze=n min1=35000 max1=100000 | Transp plane=25 > ${WOKDIR}/bpgom3d-csou-planes-subsample-big-line$*.H
	Window3d < ${DATDIR}/bpgom3d-crec-planes-subsample-all2.H n5=1 f5=$* squeeze=n min1=35000 max1=100000 | Transp plane=25 > ${WOKDIR}/bpgom3d-crec-planes-subsample-big-line$*.H
	python bpgom3d_imag2.py --bindir=${BINDIR} --datdir=${WOKDIR} --resdir=${RESDIR} --tmpdir=${TMPDIR}\
         --csou=bpgom3d-csou-planes-subsample-big-line$*.H --crec=bpgom3d-crec-planes-subsample-big-line$*.H --bvel=bpgom3d-bvel-subsample.H --bimg=$@ \
         --jsx_win=2 --nsx_tot=101 \
         --offset=n --comsou=y --oper=ssfpi --nvrf=3 --nws=16 --nxpad_beg=0 --nxpad_end=0 --nypad_beg=100 --nypad_end=100 --ntpx=20 --ntpy=50 \
         --nhx=21 --ohx=-1640.4 --dhx=164.04 --nhy=7 --ohy=-1476.36 --dhy=492.12\
         --image_xmin=50000 --image_xmax=90000 --image_ymin=65616 --image_ymax=82000 --image_zmin=15000 --image_zmax=24960 \


#################################################
#????????????????????????????????????????????????
# the shallower part
#################################################
bpgom3d-imag-top-planes-subsample-line%.H: 
	Window3d < ${DATDIR}/bpgom3d-csou-planes-subsample-all2.H n5=1 f5=$* squeeze=n min1=35000 max1=100000 | Transp plane=25 > ${WOKDIR}/bpgom3d-csou-planes-subsample-line$*.H
	Window3d < ${DATDIR}/bpgom3d-crec-planes-subsample-all2.H n5=1 f5=$* squeeze=n min1=35000 max1=100000 | Transp plane=25 > ${WOKDIR}/bpgom3d-crec-planes-subsample-line$*.H
	python bpgom3d_imag2.py --bindir=${BINDIR} --datdir=${WOKDIR} --resdir=${RESDIR} --tmpdir=${TMPDIR}\
         --csou=bpgom3d-csou-planes-subsample-line$*.H --crec=bpgom3d-crec-planes-subsample-line$*.H --bvel=bpgom3d-vmod-invt-full.H --bimg=$@ \
         --jsx_win=3 --nsx_tot=101 \
         --offset=n --comsou=y --oper=ssfpi --nvrf=5 --nws=16 --nxpad_beg=0 --nxpad_end=0 --nypad_beg=100 --nypad_end=100 --ntpx=20 --ntpy=50 \
         --nhx=21 --ohx=-1640.4 --dhx=164.04 \
         --image_xmin=50000 --image_xmax=90000 --image_ymin=65616 --image_ymax=82000 --image_zmax=14940 \
         --memchk=n --queue=default --grepinfo=line$* --datapath=/tmp/


#################################################
#################################################
# Updated velocity
#################################################
#################################################
bpgom3d-imag-invt-planes-subsample-line%.H:
	Window3d < ${DATDIR}/bpgom3d-csou-planes-subsample-all2.H n5=1 f5=$* squeeze=n min1=35000 max1=100000 | Transp plane=25 > ${WOKDIR}/bpgom3d-csou-planes-subsample-line$*.H
	Window3d < ${DATDIR}/bpgom3d-crec-planes-subsample-all2.H n5=1 f5=$* squeeze=n min1=35000 max1=100000 | Transp plane=25 > ${WOKDIR}/bpgom3d-crec-planes-subsample-line$*.H
	python bpgom3d_imag2.py --bindir=${BINDIR} --datdir=${WOKDIR} --resdir=${RESDIR} --tmpdir=${TMPDIR}\
         --csou=bpgom3d-csou-planes-subsample-line$*.H --crec=bpgom3d-crec-planes-subsample-line$*.H --bvel=bpgom3d-vmod-invt-full.H --bimg=$@ \
         --jsx_win=3 --nsx_tot=101 \
         --offset=n --comsou=y --oper=ssfpi --nvrf=5 --nws=16 --nxpad_beg=0 --nxpad_end=0 --nypad_beg=100 --nypad_end=100 --ntpx=20 --ntpy=50 \
         --nhx=21 --ohx=-1640.4 --dhx=164.04 \
         --image_xmin=50000 --image_xmax=90000 --image_ymin=65616 --image_ymax=82000 --image_zmin=15000 --image_zmax=35000 \
         --memchk=n --queue=Q35 --grepinfo=line$* --datapath=/tmp/





bar-vmod.rsf:
	Window3d < ${BPGOM3D_DATA}/bpgom3d-bvel-subsample.H n3=584 min1=42000 max1=95000 min2=65000 max2=85000 > junk1.H
	${BINDIR}/bar.x < junk1.H minval=1.4 maxval=4.6 > junk2.H
	sfbyte < junk2.H gainpanel=a pclip=100 allpos=y bar=$@ > j.T
	Rm junk*.H j.T
	

${RESDIR}/bpgom3d-bvel-full.v: bpgom3d-bvel-subsample.H
	Window3d < bpgom3d-bvel-subsample.H n3=584 min1=42000 max1=95000 min2=65000 max2=85000 > junk1.H
	${BINDIR}/ft2m.x header=junk1.H axis1=y axis2=y axis3=y km=y
	${BINDIR}/ft2m_bin.x < junk1.H km=y > junk2.H
	echo o1=0 o2=0 >> junk2.H
	Window3d < junk2.H | Transp reshape=2,3 > junk3.H
	sfbyte < junk3.H gainpanel=a allpos=y pclip=100 > j.T
	sfcubeplot < j.T bar=bar-vmod.rsf color=j title=" " label1="Depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        scalebar=y barlabel="Velocity" point1=0.64 point2=0.72 o1num=0 d1num=2 n1tic=8 n2tic=0\
        frame1=340 frame2=270 frame3=38 framelabel1=n framelabel2=n framelabel3=n barunit="km/s" labelfat=1 labelsz=6 > junk1.v
	Overlay < junk1.v x=3    y=1.28 type=box width=5.1 height=2.7 out=junk2.v fat=3 > /dev/null
	Overlay < junk1.v x=8.69 y=1.28 type=box width=2.2 height=2.7 out=junk3.v fat=3 > /dev/null
	Overlay < junk1.v x=3    y=6.23 type=box width=5.1 height=2.2 out=junk4.v fat=3 > /dev/null
	vp_Overlay junk1.v junk2.v junk3.v junk4.v > $@
	Rm junk1.H junk2.H junk3.H j.T
	rm -f junk*.v

${RESDIR}/bpgom3d-bvel-full-grey.v:
	Window3d < ${BPGOM3D_DATA}/bpgom3d-bvel-subsample.H n3=584 min1=42000 max1=95000 min2=65000 max2=85000 > junk1.H
	${BINDIR}/ft2m.x header=junk1.H axis1=y axis2=y axis3=y km=y
	${BINDIR}/ft2m_bin.x < junk1.H km=y > junk2.H
	echo o1=0 o2=0 >> junk2.H
	Window3d < junk2.H | Transp reshape=2,3 > junk3.H
	sfbyte < junk3.H gainpanel=a allpos=y pclip=100 > j.T
	sfcubeplot < j.T bar=bar-vmod.rsf title=" " label1="Depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        scalebar=y barlabel="Velocity" point1=0.64 point2=0.72 o1num=0 d1num=2 n1tic=8 n2tic=0\
        frame1=340 frame2=270 frame3=38 framelabel1=n framelabel2=n framelabel3=n barunit="km/s" labelfat=1 labelsz=6 > junk1.v
	Overlay < junk1.v x=3    y=1.28 type=box width=5.1 height=2.7 out=junk2.v fat=3 color=0 > /dev/null
	Overlay < junk1.v x=8.69 y=1.28 type=box width=2.2 height=2.7 out=junk3.v fat=3 color=0 > /dev/null
	Overlay < junk1.v x=3    y=6.23 type=box width=5.1 height=2.2 out=junk4.v fat=3 color=0 > /dev/null
	vp_Overlay junk1.v junk2.v junk3.v junk4.v > $@
	Rm junk1.H junk2.H junk3.H j.T
	rm -f junk*.v


${RESDIR}/bpgom3d-vmod-full.v:
	Window3d < ${BPGOM3D_IMAG}/bpgom3d-vmod-invt-full.H n3=584 min1=42000 max1=95000 min2=65000 max2=85000 > junk1.H
	${BINDIR}/ft2m.x header=junk1.H axis1=y axis2=y axis3=y km=y
	${BINDIR}/ft2m_bin.x < junk1.H km=y > junk2.H
	echo o1=0 o2=0 >> junk2.H
	Window3d < junk2.H | Transp reshape=2,3 > junk3.H
	sfbyte < junk3.H gainpanel=a allpos=y pclip=100 > j.T
	sfcubeplot < j.T bar=bar-vmod.rsf color=j title=" " label1="Depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        scalebar=y barlabel="Velocity" point1=0.64 point2=0.72 o1num=0 d1num=2 n1tic=8 n2tic=0\
        frame1=340 frame2=270 frame3=38 framelabel1=n framelabel2=n framelabel3=n barunit="km/s" labelfat=1 labelsz=6 > junk1.v
	Overlay < junk1.v x=3    y=1.28 type=box width=5.1 height=2.7 out=junk2.v fat=3 > /dev/null
	Overlay < junk1.v x=8.69 y=1.28 type=box width=2.2 height=2.7 out=junk3.v fat=3 > /dev/null
	Overlay < junk1.v x=3    y=6.23 type=box width=5.1 height=2.2 out=junk4.v fat=3 > /dev/null
	vp_Overlay junk1.v junk2.v junk3.v junk4.v > $@
	Rm junk1.H junk2.H junk3.H j.T
	rm -f junk*.v

${RESDIR}/bpgom3d-vmod-full-grey.v:
	Window3d < ${BPGOM3D_IMAG}/bpgom3d-vmod-invt-full.H n3=584 min1=42000 max1=95000 min2=65000 max2=85000 > junk1.H
	${BINDIR}/ft2m.x header=junk1.H axis1=y axis2=y axis3=y km=y
	${BINDIR}/ft2m_bin.x < junk1.H km=y > junk2.H
	echo o1=0 o2=0 >> junk2.H
	Window3d < junk2.H | Transp reshape=2,3 > junk3.H
	sfbyte < junk3.H gainpanel=a allpos=y pclip=100 > j.T
	sfcubeplot < j.T bar=bar-vmod.rsf title=" " label1="Depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        scalebar=y barlabel="Velocity" point1=0.64 point2=0.72 o1num=0 d1num=2 n1tic=8 n2tic=0\
        frame1=340 frame2=270 frame3=38 framelabel1=n framelabel2=n framelabel3=n barunit="km/s" labelfat=1 labelsz=6 > junk1.v
	Overlay < junk1.v x=3    y=1.28 type=box width=5.1 height=2.7 out=junk2.v fat=3 color=0 > /dev/null
	Overlay < junk1.v x=8.69 y=1.28 type=box width=2.2 height=2.7 out=junk3.v fat=3 color=0 > /dev/null
	Overlay < junk1.v x=3    y=6.23 type=box width=5.1 height=2.2 out=junk4.v fat=3 color=0 > /dev/null
	vp_Overlay junk1.v junk2.v junk3.v junk4.v > $@
	Rm junk1.H junk2.H junk3.H j.T
	rm -f junk*.v


${RESDIR}/bpgom3d-gmsk-target.v:
	Window3d < ${BPGOM3D_IMAG}/bpgom3d-gmsk-target-show.H min1=48000 max1=95000 min2=65000 max2=83000 > junk1.H
	${BINDIR}/ft2m.x header=junk1.H axis1=y axis2=y axis3=y km=y
	echo o1=1.83 o2=0 o3=0 >> junk1.H
	Window3d < junk1.H | Transp reshape=2,3 > junk2.H
	sfbyte < junk2.H gainpanel=a pclip=100 bar=bar.rsf >  j.T
	sfcubeplot < j.T bar=bar.rsf title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        scalebar=y barlabel=" " point1=0.52 point2=0.72 o1num=2 d1num=1 o2num=0 d2num=1 n1tic=14 n2tic=7 minval=0 maxval=1\
        frame1=43 frame2=235 frame3=38 framelabel1=n framelabel2=n framelabel3=n barunit=" " labelfat=1 labelsz=8 screenratio=0.48 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.6 > $@
	Rm junk1.H junk2.H junk3.H j.T
	rm -f junk*.v bar.rsf

bar-grad.rsf:
	Window3d < ${BPGOM3D_IMAG}/bpgom3d-grdo-raw-iter1.H min1=48000 max1=95000 min2=65000 max2=83000 > junk1.H
	${BINDIR}/bar.x < junk1.H minval=-1 maxval=1 > junk2.H
	sfbyte < junk2.H bar=$@ pclip=100 > i.T
	Rm junk*.H
	rm -f i.T

${RESDIR}/bpgom3d-grad-raw-target-iter%.v:
	Window3d < ${BPGOM3D_IMAG}/grad-raw-iter$*.H min1=48000 max1=95000 min2=65000 max2=83000 | ${BINDIR}/normalize3d.x > junk1.H
	${BINDIR}/ft2m.x header=junk1.H axis1=y axis2=y axis3=y km=y
	echo o1=1.83 o2=0 o3=0 >> junk1.H
	Window3d < junk1.H | Transp reshape=2,3 > junk2.H
	sfbyte < junk2.H gainpanel=a pclip=100 >  j.T
	sfcubeplot < j.T bar=bar-grad.rsf color=e title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        scalebar=y barlabel=" " point1=0.52 point2=0.72 o1num=2 d1num=1 o2num=0 d2num=1 n1tic=14 n2tic=7 minval=-1 maxval=1\
        frame1=43 frame2=235 frame3=38 framelabel1=n framelabel2=n framelabel3=n barunit=" " labelfat=1 labelsz=8 screenratio=0.48 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.6 > $@
	Rm junk1.H junk2.H junk3.H j.T
	rm -f junk*.v

${RESDIR}/bpgom3d-grad-smt-target-iter%.v:
	Window3d < ${BPGOM3D_IMAG}/grad-smt-iter$*.H min1=48000 max1=95000 min2=65000 max2=83000 | ${BINDIR}/normalize3d.x > junk1.H
	${BINDIR}/ft2m.x header=junk1.H axis1=y axis2=y axis3=y km=y
	echo o1=1.83 o2=0 o3=0 >> junk1.H
	Window3d < junk1.H | Transp reshape=2,3 > junk2.H
	sfbyte < junk2.H gainpanel=a pclip=100 >  j.T
	sfcubeplot < j.T bar=bar-grad.rsf color=e title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        scalebar=y barlabel=" " point1=0.52 point2=0.72 o1num=2 d1num=1 o2num=0 d2num=1 n1tic=14 n2tic=7 minval=-1 maxval=1\
        frame1=43 frame2=235 frame3=38 framelabel1=n framelabel2=n framelabel3=n barunit=" " labelfat=1 labelsz=8 screenratio=0.48 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.6 > $@
	Rm junk1.H junk2.H junk3.H j.T
	rm -f junk*.v

${RESDIR}/bpgom3d-grad-raw-target-grey-iter%.v:
	Window3d < ${BPGOM3D_IMAG}/grad-raw-iter$*.H min1=48000 max1=95000 min2=65000 max2=83000 | ${BINDIR}/normalize3d.x > junk1.H
	${BINDIR}/ft2m.x header=junk1.H axis1=y axis2=y axis3=y km=y
	echo o1=1.83 o2=0 o3=0 >> junk1.H
	Window3d < junk1.H | Transp reshape=2,3 > junk2.H
	sfbyte < junk2.H gainpanel=a pclip=100 >  j.T
	sfcubeplot < j.T bar=bar-grad.rsf title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        scalebar=y barlabel=" " point1=0.52 point2=0.72 o1num=2 d1num=1 o2num=0 d2num=1 n1tic=14 n2tic=7 minval=-1 maxval=1\
        frame1=43 frame2=235 frame3=38 framelabel1=n framelabel2=n framelabel3=n barunit=" " labelfat=1 labelsz=8 screenratio=0.48 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.6 > $@
	Rm junk1.H junk2.H junk3.H j.T
	rm -f junk*.v

${RESDIR}/bpgom3d-grad-smt-target-grey-iter%.v:
	Window3d < ${BPGOM3D_IMAG}/grad-smt-iter$*.H min1=48000 max1=95000 min2=65000 max2=83000 | ${BINDIR}/normalize3d.x > junk1.H
	${BINDIR}/ft2m.x header=junk1.H axis1=y axis2=y axis3=y km=y
	echo o1=1.83 o2=0 o3=0 >> junk1.H
	Window3d < junk1.H | Transp reshape=2,3 > junk2.H
	sfbyte < junk2.H gainpanel=a pclip=100 >  j.T
	sfcubeplot < j.T bar=bar-grad.rsf title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        scalebar=y barlabel=" " point1=0.52 point2=0.72 o1num=2 d1num=1 o2num=0 d2num=1 n1tic=14 n2tic=7 minval=-1 maxval=1\
        frame1=43 frame2=235 frame3=38 framelabel1=n framelabel2=n framelabel3=n barunit=" " labelfat=1 labelsz=8 screenratio=0.48 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.6 > $@
	Rm junk1.H junk2.H junk3.H j.T
	rm -f junk*.v


${RESDIR}/bpgom3d-invt-target-stage%.v:
	Window3d < ${BPGOM3D_IMAG}/bpgom3d-vmod-invt-old2-stage$*.H min1=48000 max1=95000 min2=65000 max2=83000 > junk1.H
	${BINDIR}/ft2m.x header=junk1.H axis1=y axis2=y axis3=y km=y
	echo o1=1.83 o2=0 o3=0 >> junk1.H
	Window3d < junk1.H | Transp reshape=2,3,4 > junk2.H
	sfbyte < junk2.H gainpanel=a pclip=100 allpos=y >  j.T
	sfcubeplot < j.T bar=bar-vmod.rsf color=j title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        scalebar=y barlabel="Velocity" point1=0.52 point2=0.72 o1num=2 d1num=1 o2num=0 d2num=1 n1tic=14 n2tic=7 \
        frame1=43 frame2=235 frame3=38 framelabel1=n framelabel2=n framelabel3=n barunit="km/s" labelfat=1 labelsz=8 screenratio=0.48 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.6 > $@
	Rm junk1.H junk2.H junk3.H j.T
	rm -f junk*.v


# x0 = 12.7998322, y0=19.799747
ORIG_LOC1 = frame1=270 frame2=182 frame3=16
ORIG_LOC2 = frame1=250 frame2=102 frame3=70
ORIG_LOC3 = frame1=80  frame2=58  frame3=84 92

ORIG_INT_LOC1 = frame1=540 frame2=364 frame3=32
ORIG_INT_LOC2 = frame1=500 frame2=204 frame3=140
ORIG_INT_LOC3 = frame1=160 frame2=116 frame3=168

CIG2D_LOC1 = f4=16
CIG2D_LOC2 = f4=70
CIG2D_LOC3 = f4=84 92


${RESDIR}/bpgom3d-bimg-orig%.v:
	Window3d < ${BPGOM3D_IMAG}/bpgom3d-bimg-planes-subsample-target.H n1=1 min1=0 n2=1 min2=0 | Transp reshape=2,3,5 > junk1.H
	${BINDIR}/normalize3d.x < junk1.H > junk2.H
	${BINDIR}/ft2m.x header=junk2.H axis1=y axis2=y axis3=y km=y
	echo o1=0 o2=2.45 o3=0.2 >> junk2.H
	Grey < junk2.H eout=1 gainpanel=a pclip=98 > 1.T
	sfcubeplot < 1.T title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        framelabel1=n framelabel2=n framelabel3=n point1=0.55 point2=0.71 ${ORIG_LOC$*} labelfat=1 labelsz=6 screenratio=0.6 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.8 > $@
	Rm junk1.H junk2.H 1.T
	rm -f junk1.v junk2.v

${RESDIR}/bpgom3d-bimg-cpst%.v:
	Window3d < ${BPGOM3D_IMAG}/bpgom3d-bimg-planes-subsample-target-cpst-sroot.H n1=1 min1=0 n2=1 min2=0 | Transp reshape=2,3,5 > junk1.H
	${BINDIR}/normalize3d.x < junk1.H > junk2.H
	${BINDIR}/ft2m.x header=junk2.H axis1=y axis2=y axis3=y km=y
	echo o1=0 o2=2.45 o3=0.2 >> junk2.H
	Grey < junk2.H eout=1 gainpanel=a pclip=98 > 1.T
	sfcubeplot < 1.T title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        framelabel1=n framelabel2=n framelabel3=n point1=0.55 point2=0.71 ${ORIG_LOC$*} labelfat=1 labelsz=6 screenratio=0.6 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.8 > $@
	Rm junk1.H junk2.H 1.T
	rm -f junk1.v junk2.v

${RESDIR}/bpgom3d-bimg-cond%.v:
	Window3d < ${BPGOM3D_IMAG}/bpgom3d-bimg-target-condition-modeling-sroot.H n1=1 min1=0 n2=1 min2=0 | Transp reshape=2,3,5 > junk1.H
	${BINDIR}/normalize3d.x < junk1.H > junk2.H
	${BINDIR}/ft2m.x header=junk2.H axis1=y axis2=y axis3=y km=y
	echo o1=0 o2=2.45 o3=0.2 >> junk2.H
	Grey < junk2.H eout=1 gainpanel=a pclip=98 > 1.T
	sfcubeplot < 1.T title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        framelabel1=n framelabel2=n framelabel3=n point1=0.55 point2=0.71 ${ORIG_LOC$*} labelfat=1 labelsz=6 screenratio=0.6 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.8 > $@
	Rm junk1.H junk2.H 1.T
	rm -f junk1.v junk2.v

${RESDIR}/bpgom3d-bimg-born%.v:
	Window3d < ${BPGOM3D_IMAG}/bpgom3d-bimg-born-old3.H n1=1 min1=0 n2=1 min2=0 | Transp reshape=2,3,5 > junk1.H
	${BINDIR}/normalize3d.x < junk1.H > junk2.H
	${BINDIR}/ft2m.x header=junk2.H axis1=y axis2=y axis3=y km=y
	echo o1=0 o2=2.45 o3=0.2 >> junk2.H
	Grey < junk2.H eout=1 gainpanel=a pclip=98 > 1.T
	sfcubeplot < 1.T title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        framelabel1=n framelabel2=n framelabel3=n point1=0.55 point2=0.71 ${ORIG_LOC$*} labelfat=1 labelsz=6 screenratio=0.6 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.8 > $@
	Rm junk1.H junk2.H 1.T
	rm -f junk1.v junk2.v

${RESDIR}/bpgom3d-imag-born%.v:
	Window3d < ${BPGOM3D_IMAG}/bpgom3d-imag-invt-born.H n1=1 min1=0 n2=1 min2=0 | Transp reshape=2,3,5 > junk1.H
	${BINDIR}/normalize3d.x < junk1.H > junk2.H
	${BINDIR}/ft2m.x header=junk2.H axis1=y axis2=y axis3=y km=y
	echo o1=0 o2=2.45 o3=0.2 >> junk2.H
	Grey < junk2.H eout=1 gainpanel=a pclip=98 > 1.T
	sfcubeplot < 1.T title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        framelabel1=n framelabel2=n framelabel3=n point1=0.55 point2=0.71 ${ORIG_LOC$*} labelfat=1 labelsz=6 screenratio=0.6 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.8 > $@
	Rm junk1.H junk2.H 1.T
	rm -f junk1.v junk2.v


${RESDIR}/bpgom3d-imag-orig%.v:
	Window3d < ${BPGOM3D_IMAG}/bpgom3d-imag-invt-planes-subsample.H n1=1 min1=0 n2=1 min2=0 | Transp reshape=2,3,5 > junk1.H
	${BINDIR}/normalize3d.x < junk1.H > junk2.H
	${BINDIR}/ft2m.x header=junk2.H axis1=y axis2=y axis3=y km=y
	echo o1=0 o2=2.45 o3=0.2 >> junk2.H
	Grey < junk2.H eout=1 gainpanel=a pclip=98 > 1.T
	sfcubeplot < 1.T title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        framelabel1=n framelabel2=n framelabel3=n point1=0.55 point2=0.71 ${ORIG_LOC$*} labelfat=1 labelsz=6 screenratio=0.6 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.8 > $@
	Rm junk1.H junk2.H 1.T
	rm -f junk1.v junk2.v

${RESDIR}/bpgom3d-imag-cpst%.v:
	Window3d < ${BPGOM3D_IMAG}/bpgom3d-imag-invt-target-cpst.H n1=1 min1=0 n2=1 min2=0 | Transp reshape=2,3,5 > junk1.H
	${BINDIR}/normalize3d.x < junk1.H > junk2.H
	${BINDIR}/ft2m.x header=junk2.H axis1=y axis2=y axis3=y km=y
	echo o1=0 o2=2.45 o3=0.2 >> junk2.H
	Grey < junk2.H eout=1 gainpanel=a pclip=98 > 1.T
	sfcubeplot < 1.T title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        framelabel1=n framelabel2=n framelabel3=n point1=0.55 point2=0.71 ${ORIG_LOC$*} labelfat=1 labelsz=6 screenratio=0.6 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.8 > $@
	Rm junk1.H junk2.H 1.T
	rm -f junk1.v junk2.v

${RESDIR}/bpgom3d-imag2-cpst%.v:
	Window3d < ${BPGOM3D_IMAG}/bpgom3d-imag-invt-target-cpst.H n1=1 min1=0 n2=1 min2=0 | Transp reshape=2,3,5 > junk1.H
	Interp < junk1.H d1out=30 > junk11.H
	Transp < junk11.H          maxsize=2000 | Interp d1out=82.02 | Transp          maxsize=2000 > junk12.H
	Transp < junk12.H plane=13 maxsize=2000 | Interp d1out=82.02 | Transp plane=13 maxsize=2000 > junk13.H
	${BINDIR}/normalize3d.x < junk13.H > junk2.H
	${BINDIR}/ft2m.x header=junk2.H axis1=y axis2=y axis3=y km=y
	echo o1=0 o2=2.45 o3=0.2 >> junk2.H
	Grey < junk2.H eout=1 gainpanel=a pclip=98 > 1.T
	sfcubeplot < 1.T title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        framelabel1=n framelabel2=n framelabel3=n point1=0.55 point2=0.71 ${ORIG_INT_LOC$*} labelfat=1 labelsz=6 screenratio=0.6 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.8 > $@
	Rm junk1.H junk2.H 1.T
	rm -f junk1.v junk2.v

${RESDIR}/bpgom3d-bimg2-cpst%.v:
	Window3d < ${BPGOM3D_IMAG}/bpgom3d-bimg-planes-subsample-target-cpst-sroot.H n1=1 min1=0 n2=1 min2=0 | Transp reshape=2,3,5 > junk1.H
	Interp < junk1.H d1out=30 > junk11.H
	Transp < junk11.H          maxsize=2000 | Interp d1out=82.02 | Transp          maxsize=2000 > junk12.H
	Transp < junk12.H plane=13 maxsize=2000 | Interp d1out=82.02 | Transp plane=13 maxsize=2000 > junk13.H
	${BINDIR}/normalize3d.x < junk13.H > junk2.H
	${BINDIR}/ft2m.x header=junk2.H axis1=y axis2=y axis3=y km=y
	echo o1=0 o2=2.45 o3=0.2 >> junk2.H
	Grey < junk2.H eout=1 gainpanel=a pclip=98 > 1.T
	sfcubeplot < 1.T title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        framelabel1=n framelabel2=n framelabel3=n point1=0.55 point2=0.71 ${ORIG_INT_LOC$*} labelfat=1 labelsz=6 screenratio=0.6 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.8 > $@
	Rm junk1.H junk2.H 1.T
	rm -f junk1.v junk2.v




${RESDIR}/bpgom3d-bhes%.v:
	Window3d < ${BPGOM3D_IMAG}/bpgom3d-bhes-planes-subsample.H n1=1 min1=0 n2=1 min2=0 | Transp reshape=2,3,5 > junk1.H
	${BINDIR}/normalize3d.x < junk1.H > junk2.H
	${BINDIR}/ft2m.x header=junk2.H axis1=y axis2=y axis3=y km=y
	echo o1=0 o2=2.45 o3=0.2 >> junk2.H
	Grey < junk2.H eout=1 gainpanel=a pclip=98 allpos=y > 1.T
	sfcubeplot < 1.T color=j title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        framelabel1=n framelabel2=n framelabel3=n point1=0.55 point2=0.71 ${ORIG_LOC$*} labelfat=1 labelsz=6 screenratio=0.6 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.8 > $@
	Rm junk1.H junk2.H 1.T
	rm -f junk1.v junk2.v

${RESDIR}/bpgom3d-hess%.v:
	Window3d < ${BPGOM3D_IMAG}/bpgom3d-hess-diag-target.H n1=1 min1=0 n2=1 min2=0 | Transp reshape=2,3,5 > junk1.H
	${BINDIR}/normalize3d.x < junk1.H > junk2.H
	${BINDIR}/ft2m.x header=junk2.H axis1=y axis2=y axis3=y km=y
	echo o1=0 o2=2.45 o3=0.2 >> junk2.H
	Grey < junk2.H eout=1 gainpanel=a pclip=98 allpos=y > 1.T
	sfcubeplot < 1.T color=j title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        framelabel1=n framelabel2=n framelabel3=n point1=0.55 point2=0.71 ${ORIG_LOC$*} labelfat=1 labelsz=6 screenratio=0.6 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.8 > $@
	Rm junk1.H junk2.H 1.T
	rm -f junk1.v junk2.v


${RESDIR}/cig2d-bvel-cpst%.v:
	Window3d < ${BPGOM3D_IMAG}/bpgom3d-bimg-target-condition-modeling-sroot.H n2=1 min2=0 n4=1 ${CIG2D_LOC$*} | Transp plane=13 > junk1.H
	${BINDIR}/normalize3d.x < junk1.H > junk2.H
	${BINDIR}/ft2m.x header=junk2.H axis1=y axis2=n axis3=y km=y
	echo o1=0 >> junk2.H
	Window3d < junk2.H n2=1 min2=55000 | Grey title=" " wantaxis1=y labelfat=2 labelsz=8 label1="Relative depth (km)" label2="hx (km)" \
        out=junk1.v screenratio=3.3 labelrot=n > /dev/null
	Window3d < junk2.H n2=1 min2=60000 | Grey title=" " wantaxis1=n labelfat=2 labelsz=8 label2="hx (km)" out=junk2.v screenratio=3.3 > /dev/null
	Window3d < junk2.H n2=1 min2=65000 | Grey title=" " wantaxis1=n labelfat=2 labelsz=8 label2="hx (km)" out=junk3.v screenratio=3.3 > /dev/null
	Window3d < junk2.H n2=1 min2=70000 | Grey title=" " wantaxis1=n labelfat=2 labelsz=8 label2="hx (km)" out=junk4.v screenratio=3.3 > /dev/null
	Window3d < junk2.H n2=1 min2=75000 | Grey title=" " wantaxis1=n labelfat=2 labelsz=8 label2="hx (km)" out=junk5.v screenratio=3.3 > /dev/null
	Window3d < junk2.H n2=1 min2=80000 | Grey title=" " wantaxis1=n labelfat=2 labelsz=8 label2="hx (km)" out=junk6.v screenratio=3.3 > /dev/null
	Window3d < junk2.H n2=1 min2=85000 | Grey title=" " wantaxis1=n labelfat=2 labelsz=8 label2="hx (km)" out=junk7.v screenratio=3.3 > /dev/null
	vppen < junk1.v vpstyle=n align=lb | vppen xshift=0.65 yshift=1 > 1.v
	vppen < junk2.v vpstyle=n align=lb | vppen xshift=4.3  yshift=1. > 2.v
	vppen < junk3.v vpstyle=n align=lb | vppen xshift=7.28  yshift=1 > 3.v
	vppen < junk4.v vpstyle=n align=lb | vppen xshift=10.26 yshift=1 > 4.v
	vppen < junk5.v vpstyle=n align=lb | vppen xshift=13.24 yshift=1 > 5.v
	vppen < junk6.v vpstyle=n align=lb | vppen xshift=16.22 yshift=1 > 6.v
	vppen < junk7.v vpstyle=n align=lb | vppen xshift=19.2 yshift=1 > 7.v
	vppen 1.v 2.v 3.v 4.v 5.v 6.v 7.v erase=n scale=0.2 > $@
	Rm junk1.H junk2.H
	rm -f junk[0-7].v [1-7].v

${RESDIR}/cig2d-invt-cpst%.v:
	Window3d < ${BPGOM3D_IMAG}/bpgom3d-imag-invt-target-cpst.H n2=1 min2=0 n4=1 ${CIG2D_LOC$*} | Transp plane=13 > junk1.H
	${BINDIR}/normalize3d.x < junk1.H > junk2.H
	${BINDIR}/ft2m.x header=junk2.H axis1=y axis2=n axis3=y km=y
	echo o1=0 >> junk2.H
	Window3d < junk2.H n2=1 min2=55000 | Grey title=" " wantaxis1=y labelfat=2 labelsz=8 label1="Relative depth (km)" label2="hx (km)" \
        out=junk1.v screenratio=3.3 labelrot=n > /dev/null
	Window3d < junk2.H n2=1 min2=60000 | Grey title=" " wantaxis1=n labelfat=2 labelsz=8 label2="hx (km)" out=junk2.v screenratio=3.3 > /dev/null
	Window3d < junk2.H n2=1 min2=65000 | Grey title=" " wantaxis1=n labelfat=2 labelsz=8 label2="hx (km)" out=junk3.v screenratio=3.3 > /dev/null
	Window3d < junk2.H n2=1 min2=70000 | Grey title=" " wantaxis1=n labelfat=2 labelsz=8 label2="hx (km)" out=junk4.v screenratio=3.3 > /dev/null
	Window3d < junk2.H n2=1 min2=75000 | Grey title=" " wantaxis1=n labelfat=2 labelsz=8 label2="hx (km)" out=junk5.v screenratio=3.3 > /dev/null
	Window3d < junk2.H n2=1 min2=80000 | Grey title=" " wantaxis1=n labelfat=2 labelsz=8 label2="hx (km)" out=junk6.v screenratio=3.3 > /dev/null
	Window3d < junk2.H n2=1 min2=85000 | Grey title=" " wantaxis1=n labelfat=2 labelsz=8 label2="hx (km)" out=junk7.v screenratio=3.3 > /dev/null
	vppen < junk1.v vpstyle=n align=lb | vppen xshift=0.65 yshift=1 > 1.v
	vppen < junk2.v vpstyle=n align=lb | vppen xshift=4.3  yshift=1. > 2.v
	vppen < junk3.v vpstyle=n align=lb | vppen xshift=7.28  yshift=1 > 3.v
	vppen < junk4.v vpstyle=n align=lb | vppen xshift=10.26 yshift=1 > 4.v
	vppen < junk5.v vpstyle=n align=lb | vppen xshift=13.24 yshift=1 > 5.v
	vppen < junk6.v vpstyle=n align=lb | vppen xshift=16.22 yshift=1 > 6.v
	vppen < junk7.v vpstyle=n align=lb | vppen xshift=19.2 yshift=1 > 7.v
	vppen 1.v 2.v 3.v 4.v 5.v 6.v 7.v erase=n scale=0.2 > $@
	Rm junk1.H junk2.H
	rm -f junk[0-7].v [1-7].v


${RESDIR}/cig2d-bvel-born%.v:
	Window3d < ${BPGOM3D_IMAG}/bpgom3d-bimg-born-old3.H n2=1 min2=0 n4=1 ${CIG2D_LOC$*} | Transp plane=13 > junk1.H
	${BINDIR}/normalize3d.x < junk1.H > junk2.H
	${BINDIR}/ft2m.x header=junk2.H axis1=y axis2=n axis3=y km=y
	echo o1=0 >> junk2.H
	Window3d < junk2.H n2=1 min2=55000 | Grey title=" " wantaxis1=y labelfat=2 labelsz=8 label1="Relative depth (km)" label2="hx (km)" \
        out=junk1.v screenratio=3.3 labelrot=n > /dev/null
	Window3d < junk2.H n2=1 min2=60000 | Grey title=" " wantaxis1=n labelfat=2 labelsz=8 label2="hx (km)" out=junk2.v screenratio=3.3 > /dev/null
	Window3d < junk2.H n2=1 min2=65000 | Grey title=" " wantaxis1=n labelfat=2 labelsz=8 label2="hx (km)" out=junk3.v screenratio=3.3 > /dev/null
	Window3d < junk2.H n2=1 min2=70000 | Grey title=" " wantaxis1=n labelfat=2 labelsz=8 label2="hx (km)" out=junk4.v screenratio=3.3 > /dev/null
	Window3d < junk2.H n2=1 min2=75000 | Grey title=" " wantaxis1=n labelfat=2 labelsz=8 label2="hx (km)" out=junk5.v screenratio=3.3 > /dev/null
	Window3d < junk2.H n2=1 min2=80000 | Grey title=" " wantaxis1=n labelfat=2 labelsz=8 label2="hx (km)" out=junk6.v screenratio=3.3 > /dev/null
	Window3d < junk2.H n2=1 min2=85000 | Grey title=" " wantaxis1=n labelfat=2 labelsz=8 label2="hx (km)" out=junk7.v screenratio=3.3 > /dev/null
	vppen < junk1.v vpstyle=n align=lb | vppen xshift=0.65 yshift=1 > 1.v
	vppen < junk2.v vpstyle=n align=lb | vppen xshift=4.3  yshift=1. > 2.v
	vppen < junk3.v vpstyle=n align=lb | vppen xshift=7.28  yshift=1 > 3.v
	vppen < junk4.v vpstyle=n align=lb | vppen xshift=10.26 yshift=1 > 4.v
	vppen < junk5.v vpstyle=n align=lb | vppen xshift=13.24 yshift=1 > 5.v
	vppen < junk6.v vpstyle=n align=lb | vppen xshift=16.22 yshift=1 > 6.v
	vppen < junk7.v vpstyle=n align=lb | vppen xshift=19.2 yshift=1 > 7.v
	vppen 1.v 2.v 3.v 4.v 5.v 6.v 7.v erase=n scale=0.2 > $@
	Rm junk1.H junk2.H
	rm -f junk[0-7].v [1-7].v

${RESDIR}/cig2d-invt-born%.v:
	Window3d < ${BPGOM3D_IMAG}/bpgom3d-imag-invt-born.H n2=1 min2=0 n4=1 ${CIG2D_LOC$*} | Transp plane=13 > junk1.H
	${BINDIR}/normalize3d.x < junk1.H > junk2.H
	${BINDIR}/ft2m.x header=junk2.H axis1=y axis2=n axis3=y km=y
	echo o1=0 >> junk2.H
	Window3d < junk2.H n2=1 min2=55000 | Grey title=" " wantaxis1=y labelfat=2 labelsz=8 label1="Relative depth (km)" label2="hx (km)" \
        out=junk1.v screenratio=3.3 labelrot=n > /dev/null
	Window3d < junk2.H n2=1 min2=60000 | Grey title=" " wantaxis1=n labelfat=2 labelsz=8 label2="hx (km)" out=junk2.v screenratio=3.3 > /dev/null
	Window3d < junk2.H n2=1 min2=65000 | Grey title=" " wantaxis1=n labelfat=2 labelsz=8 label2="hx (km)" out=junk3.v screenratio=3.3 > /dev/null
	Window3d < junk2.H n2=1 min2=70000 | Grey title=" " wantaxis1=n labelfat=2 labelsz=8 label2="hx (km)" out=junk4.v screenratio=3.3 > /dev/null
	Window3d < junk2.H n2=1 min2=75000 | Grey title=" " wantaxis1=n labelfat=2 labelsz=8 label2="hx (km)" out=junk5.v screenratio=3.3 > /dev/null
	Window3d < junk2.H n2=1 min2=80000 | Grey title=" " wantaxis1=n labelfat=2 labelsz=8 label2="hx (km)" out=junk6.v screenratio=3.3 > /dev/null
	Window3d < junk2.H n2=1 min2=85000 | Grey title=" " wantaxis1=n labelfat=2 labelsz=8 label2="hx (km)" out=junk7.v screenratio=3.3 > /dev/null
	vppen < junk1.v vpstyle=n align=lb | vppen xshift=0.65 yshift=1 > 1.v
	vppen < junk2.v vpstyle=n align=lb | vppen xshift=4.3  yshift=1. > 2.v
	vppen < junk3.v vpstyle=n align=lb | vppen xshift=7.28  yshift=1 > 3.v
	vppen < junk4.v vpstyle=n align=lb | vppen xshift=10.26 yshift=1 > 4.v
	vppen < junk5.v vpstyle=n align=lb | vppen xshift=13.24 yshift=1 > 5.v
	vppen < junk6.v vpstyle=n align=lb | vppen xshift=16.22 yshift=1 > 6.v
	vppen < junk7.v vpstyle=n align=lb | vppen xshift=19.2 yshift=1 > 7.v
	vppen 1.v 2.v 3.v 4.v 5.v 6.v 7.v erase=n scale=0.2 > $@
	Rm junk1.H junk2.H
	rm -f junk[0-7].v [1-7].v


${RESDIR}/adcig2d-bvel-cpst%.v:
	Window3d < ${BPGOM3D_IMAG}/bpgom3d-bimg-adcig-cpst.H n2=1 min2=0 n4=1 ${CIG2D_LOC$*} | Transp plane=13 > junk1.H
	${BINDIR}/normalize3d.x < junk1.H > junk2.H
	${BINDIR}/ft2m.x header=junk2.H axis1=y axis2=n axis3=n km=y
	echo o1=0 >> junk2.H
	Window3d < junk2.H n2=1 min2=55000 | Grey title=" " wantaxis1=y labelfat=2 labelsz=8 label1="Relative depth (km)" label2="\F10 g \F3 (deg.)" \
        out=junk1.v screenratio=3.3 o2num=-30 d2num=30 labelrot=n > /dev/null
	Window3d < junk2.H n2=1 min2=60000 | Grey title=" " wantaxis1=n labelfat=2 labelsz=8 label2="\F10 g \F3 (deg.)" out=junk2.v screenratio=3.3 o2num=-30 d2num=30 > /dev/null
	Window3d < junk2.H n2=1 min2=65000 | Grey title=" " wantaxis1=n labelfat=2 labelsz=8 label2="\F10 g \F3 (deg.)" out=junk3.v screenratio=3.3 o2num=-30 d2num=30 > /dev/null
	Window3d < junk2.H n2=1 min2=70000 | Grey title=" " wantaxis1=n labelfat=2 labelsz=8 label2="\F10 g \F3 (deg.)" out=junk4.v screenratio=3.3 o2num=-30 d2num=30 > /dev/null
	Window3d < junk2.H n2=1 min2=75000 | Grey title=" " wantaxis1=n labelfat=2 labelsz=8 label2="\F10 g \F3 (deg.)" out=junk5.v screenratio=3.3 o2num=-30 d2num=30 > /dev/null
	Window3d < junk2.H n2=1 min2=80000 | Grey title=" " wantaxis1=n labelfat=2 labelsz=8 label2="\F10 g \F3 (deg.)" out=junk6.v screenratio=3.3 o2num=-30 d2num=30 > /dev/null
	Window3d < junk2.H n2=1 min2=85000 | Grey title=" " wantaxis1=n labelfat=2 labelsz=8 label2="\F10 g \F3 (deg.)" out=junk7.v screenratio=3.3 o2num=-30 d2num=30 > /dev/null
	vppen < junk1.v vpstyle=n align=lb | vppen xshift=0.65 yshift=1 > 1.v
	vppen < junk2.v vpstyle=n align=lb | vppen xshift=4.3  yshift=1. > 2.v
	vppen < junk3.v vpstyle=n align=lb | vppen xshift=7.28  yshift=1 > 3.v
	vppen < junk4.v vpstyle=n align=lb | vppen xshift=10.26 yshift=1 > 4.v
	vppen < junk5.v vpstyle=n align=lb | vppen xshift=13.24 yshift=1 > 5.v
	vppen < junk6.v vpstyle=n align=lb | vppen xshift=16.22 yshift=1 > 6.v
	vppen < junk7.v vpstyle=n align=lb | vppen xshift=19.2 yshift=1 > 7.v
	vppen 1.v 2.v 3.v 4.v 5.v 6.v 7.v erase=n scale=0.2 > $@
	Rm junk1.H junk2.H
	rm -f junk[0-7].v [1-7].v

${RESDIR}/adcig2d-invt-cpst%.v:
	Window3d < ${BPGOM3D_IMAG}/bpgom3d-imag-adcig-cpst.H n2=1 min2=0 n4=1 ${CIG2D_LOC$*} | Transp plane=13 > junk1.H
	${BINDIR}/normalize3d.x < junk1.H > junk2.H
	${BINDIR}/ft2m.x header=junk2.H axis1=y axis2=n axis3=n km=y
	echo o1=0 >> junk2.H
	Window3d < junk2.H n2=1 min2=55000 | Grey title=" " wantaxis1=y labelfat=2 labelsz=8 label1="Relative depth (km)" label2="\F10 g \F3 (deg.)" \
        out=junk1.v screenratio=3.3 o2num=-30 d2num=30 labelrot=n > /dev/null
	Window3d < junk2.H n2=1 min2=60000 | Grey title=" " wantaxis1=n labelfat=2 labelsz=8 label2="\F10 g \F3 (deg.)" out=junk2.v screenratio=3.3 o2num=-30 d2num=30 > /dev/null
	Window3d < junk2.H n2=1 min2=65000 | Grey title=" " wantaxis1=n labelfat=2 labelsz=8 label2="\F10 g \F3 (deg.)" out=junk3.v screenratio=3.3 o2num=-30 d2num=30 > /dev/null
	Window3d < junk2.H n2=1 min2=70000 | Grey title=" " wantaxis1=n labelfat=2 labelsz=8 label2="\F10 g \F3 (deg.)" out=junk4.v screenratio=3.3 o2num=-30 d2num=30 > /dev/null
	Window3d < junk2.H n2=1 min2=75000 | Grey title=" " wantaxis1=n labelfat=2 labelsz=8 label2="\F10 g \F3 (deg.)" out=junk5.v screenratio=3.3 o2num=-30 d2num=30 > /dev/null
	Window3d < junk2.H n2=1 min2=80000 | Grey title=" " wantaxis1=n labelfat=2 labelsz=8 label2="\F10 g \F3 (deg.)" out=junk6.v screenratio=3.3 o2num=-30 d2num=30 > /dev/null
	Window3d < junk2.H n2=1 min2=85000 | Grey title=" " wantaxis1=n labelfat=2 labelsz=8 label2="\F10 g \F3 (deg.)" out=junk7.v screenratio=3.3 o2num=-30 d2num=30 > /dev/null
	vppen < junk1.v vpstyle=n align=lb | vppen xshift=0.65 yshift=1 > 1.v
	vppen < junk2.v vpstyle=n align=lb | vppen xshift=4.3  yshift=1. > 2.v
	vppen < junk3.v vpstyle=n align=lb | vppen xshift=7.28  yshift=1 > 3.v
	vppen < junk4.v vpstyle=n align=lb | vppen xshift=10.26 yshift=1 > 4.v
	vppen < junk5.v vpstyle=n align=lb | vppen xshift=13.24 yshift=1 > 5.v
	vppen < junk6.v vpstyle=n align=lb | vppen xshift=16.22 yshift=1 > 6.v
	vppen < junk7.v vpstyle=n align=lb | vppen xshift=19.2 yshift=1 > 7.v
	vppen 1.v 2.v 3.v 4.v 5.v 6.v 7.v erase=n scale=0.2 > $@
	Rm junk1.H junk2.H
	rm -f junk[0-7].v [1-7].v




${RESDIR}/bpgom3d-cig3d-orig.v:
	Window3d < ${BPGOM3D_IMAG}/bpgom3d-bimg-planes-subsample-target-interp.H n3=1 min3=60000 n4=1 min4=70000 | Transp reshape=2,3,5 > junk11.H
	Window3d < ${BPGOM3D_IMAG}/bpgom3d-bimg-planes-subsample-target-interp.H n3=1 min3=70000 n4=1 min4=75000 | Transp reshape=2,3,5 > junk21.H
	Window3d < ${BPGOM3D_IMAG}/bpgom3d-bimg-planes-subsample-target-interp.H n3=1 min3=80000 n4=1 min4=80000 | Transp reshape=2,3,5 > junk31.H
	${BINDIR}/normalize3d.x < junk11.H > junk12.H
	${BINDIR}/normalize3d.x < junk21.H > junk22.H
	${BINDIR}/normalize3d.x < junk31.H > junk32.H
	${BINDIR}/ft2m.x header=junk12.H axis1=y axis2=y axis3=y km=y
	${BINDIR}/ft2m.x header=junk22.H axis1=y axis2=y axis3=y km=y
	${BINDIR}/ft2m.x header=junk32.H axis1=y axis2=y axis3=y km=y
	echo o1=0 >> junk12.H
	echo o1=0 >> junk22.H
	echo o1=0 >> junk32.H
	Grey < junk12.H eout=1 gainpanel=a pclip=98 > 1.T
	Grey < junk22.H eout=1 gainpanel=a pclip=98 > 2.T
	Grey < junk32.H eout=1 gainpanel=a pclip=98 > 3.T
	sfcubeplot < 1.T point1=0.80 point2=0.52 screenratio=1.6 frame1=40 frame2=10 frame3=9 title=" " \
        label1="Relative depth" unit1="km" label2="hx" unit2="km" label3="hy" unit3="km" framelabel1=n framelabel2=n framelabel3=n \
        labelfat=1 labelsz=8 o1num=-0.3 d1num=0.3 n1tic=3 o3num=-0.3 d3num=0.3 n3tic=3 > junk1.v
	sfcubeplot < 2.T point1=0.80 point2=0.52 screenratio=1.6 frame1=45 frame2=10 frame3=9 title=" " \
        wantaxis1=n label2="hx" unit2="km" label3="hy" unit3="km" framelabel1=n framelabel2=n framelabel3=n \
        labelfat=1 labelsz=8 o1num=-0.3 d1num=0.3 n1tic=3 o3num=-0.3 d3num=0.3 n3tic=3 > junk2.v
	sfcubeplot < 3.T point1=0.80 point2=0.52 screenratio=1.6 frame1=39 frame2=10 frame3=9 title=" " \
        wantaxis1=n label2="hx" unit2="km" label3="hy" unit3="km" framelabel1=n framelabel2=n framelabel3=n \
        labelfat=1 labelsz=8 o1num=-0.3 d1num=0.3 n1tic=3 o3num=-0.3 d3num=0.3 n3tic=3 > junk3.v
	vppen < junk1.v vpstyle=n align=lu | vppen xshift=1.5   yshift=0 > 1.v
	vppen < junk2.v vpstyle=n align=lu | vppen xshift=7.49  yshift=0 > 2.v
	vppen < junk3.v vpstyle=n align=lu | vppen xshift=12.5  yshift=0 > 3.v
	vppen 1.v 2.v 3.v erase=n scale=0.3 > $@
	Rm junk*.H 1.T 2.T 3.T
	rm -f junk1.v junk2.v junk3.v 1.v 2.v 3.v

${RESDIR}/bpgom3d-cig3d-orig-cpst.v:
	Window3d < ${BPGOM3D_IMAG}/bpgom3d-bimg-planes-subsample-target-cpst-interp-sroot.H n3=1 min3=60000 n4=1 min4=70000 | Transp reshape=2,3,5 > junk11.H
	Window3d < ${BPGOM3D_IMAG}/bpgom3d-bimg-planes-subsample-target-cpst-interp-sroot.H n3=1 min3=70000 n4=1 min4=75000 | Transp reshape=2,3,5 > junk21.H
	Window3d < ${BPGOM3D_IMAG}/bpgom3d-bimg-planes-subsample-target-cpst-interp-sroot.H n3=1 min3=80000 n4=1 min4=80000 | Transp reshape=2,3,5 > junk31.H
	${BINDIR}/normalize3d.x < junk11.H > junk12.H
	${BINDIR}/normalize3d.x < junk21.H > junk22.H
	${BINDIR}/normalize3d.x < junk31.H > junk32.H
	${BINDIR}/ft2m.x header=junk12.H axis1=y axis2=y axis3=y km=y
	${BINDIR}/ft2m.x header=junk22.H axis1=y axis2=y axis3=y km=y
	${BINDIR}/ft2m.x header=junk32.H axis1=y axis2=y axis3=y km=y
	echo o1=0 >> junk12.H
	echo o1=0 >> junk22.H
	echo o1=0 >> junk32.H
	Grey < junk12.H eout=1 gainpanel=a pclip=98 > 1.T
	Grey < junk22.H eout=1 gainpanel=a pclip=98 > 2.T
	Grey < junk32.H eout=1 gainpanel=a pclip=98 > 3.T
	sfcubeplot < 1.T point1=0.80 point2=0.52 screenratio=1.6 frame1=40 frame2=10 frame3=9 title=" " \
        label1="Relative depth" unit1="km" label2="hx" unit2="km" label3="hy" unit3="km" framelabel1=n framelabel2=n framelabel3=n \
        labelfat=1 labelsz=8 o1num=-0.3 d1num=0.3 n1tic=3 o3num=-0.3 d3num=0.3 n3tic=3 > junk1.v
	sfcubeplot < 2.T point1=0.80 point2=0.52 screenratio=1.6 frame1=45 frame2=10 frame3=9 title=" " \
        wantaxis1=n label2="hx" unit2="km" label3="hy" unit3="km" framelabel1=n framelabel2=n framelabel3=n \
        labelfat=1 labelsz=8 o1num=-0.3 d1num=0.3 n1tic=3 o3num=-0.3 d3num=0.3 n3tic=3 > junk2.v
	sfcubeplot < 3.T point1=0.80 point2=0.52 screenratio=1.6 frame1=39 frame2=10 frame3=9 title=" " \
        wantaxis1=n label2="hx" unit2="km" label3="hy" unit3="km" framelabel1=n framelabel2=n framelabel3=n \
        labelfat=1 labelsz=8 o1num=-0.3 d1num=0.3 n1tic=3 o3num=-0.3 d3num=0.3 n3tic=3 > junk3.v
	vppen < junk1.v vpstyle=n align=lu | vppen xshift=1.5   yshift=0 > 1.v
	vppen < junk2.v vpstyle=n align=lu | vppen xshift=7.49  yshift=0 > 2.v
	vppen < junk3.v vpstyle=n align=lu | vppen xshift=12.5  yshift=0 > 3.v
	vppen 1.v 2.v 3.v erase=n scale=0.3 > $@
	Rm junk*.H 1.T 2.T 3.T
	rm -f junk1.v junk2.v junk3.v 1.v 2.v 3.v

${RESDIR}/bpgom3d-cig3d-born.v:
	Window3d < ${BPGOM3D_IMAG}/bpgom3d-bimg-born-hxhy-line1.H n3=1 min3=60000 | Transp reshape=2,3,5 | Interp d1out=60 n1out=334 > junk11.H
	Window3d < ${BPGOM3D_IMAG}/bpgom3d-bimg-born-hxhy-line2.H n3=1 min3=70000 | Transp reshape=2,3,5 | Interp d1out=60 n1out=334 > junk21.H
	Window3d < ${BPGOM3D_IMAG}/bpgom3d-bimg-born-hxhy-line3.H n3=1 min3=80000 | Transp reshape=2,3,5 | Interp d1out=60 n1out=334 > junk31.H
	${BINDIR}/normalize3d.x < junk11.H > junk12.H
	${BINDIR}/normalize3d.x < junk21.H > junk22.H
	${BINDIR}/normalize3d.x < junk31.H > junk32.H
	${BINDIR}/ft2m.x header=junk12.H axis1=y axis2=y axis3=y km=y
	${BINDIR}/ft2m.x header=junk22.H axis1=y axis2=y axis3=y km=y
	${BINDIR}/ft2m.x header=junk32.H axis1=y axis2=y axis3=y km=y
	echo o1=0 >> junk12.H
	echo o1=0 >> junk22.H
	echo o1=0 >> junk32.H
	Grey < junk12.H eout=1 gainpanel=a pclip=98 > 1.T
	Grey < junk22.H eout=1 gainpanel=a pclip=98 > 2.T
	Grey < junk32.H eout=1 gainpanel=a pclip=98 > 3.T
	sfcubeplot < 1.T point1=0.80 point2=0.52 screenratio=1.6 frame1=40 frame2=10 frame3=9 title=" " \
        label1="Relative depth" unit1="km" label2="hx" unit2="km" label3="hy" unit3="km" framelabel1=n framelabel2=n framelabel3=n \
        labelfat=1 labelsz=8 o1num=-0.3 d1num=0.3 n1tic=3 o3num=-0.3 d3num=0.3 n3tic=3 > junk1.v
	sfcubeplot < 2.T point1=0.80 point2=0.52 screenratio=1.6 frame1=45 frame2=10 frame3=9 title=" " \
        wantaxis1=n label2="hx" unit2="km" label3="hy" unit3="km" framelabel1=n framelabel2=n framelabel3=n \
        labelfat=1 labelsz=8 o1num=-0.3 d1num=0.3 n1tic=3 o3num=-0.3 d3num=0.3 n3tic=3 > junk2.v
	sfcubeplot < 3.T point1=0.80 point2=0.52 screenratio=1.6 frame1=39 frame2=10 frame3=9 title=" " \
        wantaxis1=n label2="hx" unit2="km" label3="hy" unit3="km" framelabel1=n framelabel2=n framelabel3=n \
        labelfat=1 labelsz=8 o1num=-0.3 d1num=0.3 n1tic=3 o3num=-0.3 d3num=0.3 n3tic=3 > junk3.v
	vppen < junk1.v vpstyle=n align=lu | vppen xshift=1.5   yshift=0 > 1.v
	vppen < junk2.v vpstyle=n align=lu | vppen xshift=7.49  yshift=0 > 2.v
	vppen < junk3.v vpstyle=n align=lu | vppen xshift=12.5  yshift=0 > 3.v
	vppen 1.v 2.v 3.v erase=n scale=0.3 > $@
	Rm junk*.H 1.T 2.T 3.T
	rm -f junk1.v junk2.v junk3.v 1.v 2.v 3.v

${RESDIR}/bpgom3d-cig3d-bhes.v:
	Window3d < ${BPGOM3D_IMAG}/bpgom3d-bhes-planes-subsample-interp.H n3=1 min3=60000 n4=1 min4=70000 | Transp reshape=2,3,5 > junk11.H
	Window3d < ${BPGOM3D_IMAG}/bpgom3d-bhes-planes-subsample-interp.H n3=1 min3=70000 n4=1 min4=75000 | Transp reshape=2,3,5 > junk21.H
	Window3d < ${BPGOM3D_IMAG}/bpgom3d-bhes-planes-subsample-interp.H n3=1 min3=80000 n4=1 min4=80000 | Transp reshape=2,3,5 > junk31.H
	${BINDIR}/normalize3d.x < junk11.H > junk12.H
	${BINDIR}/normalize3d.x < junk21.H > junk22.H
	${BINDIR}/normalize3d.x < junk31.H > junk32.H
	${BINDIR}/ft2m.x header=junk12.H axis1=y axis2=y axis3=y km=y
	${BINDIR}/ft2m.x header=junk22.H axis1=y axis2=y axis3=y km=y
	${BINDIR}/ft2m.x header=junk32.H axis1=y axis2=y axis3=y km=y
	echo o1=0 >> junk12.H
	echo o1=0 >> junk22.H
	echo o1=0 >> junk32.H
	Grey < junk12.H eout=1 gainpanel=a pclip=98 allpos=y > 1.T
	Grey < junk22.H eout=1 gainpanel=a pclip=98 allpos=y > 2.T
	Grey < junk32.H eout=1 gainpanel=a pclip=98 allpos=y > 3.T
	sfcubeplot < 1.T point1=0.80 point2=0.52 screenratio=1.6 frame1=40 frame2=10 frame3=9 title=" " color=j \
        label1="Relative depth" unit1="km" label2="hx" unit2="km" label3="hy" unit3="km" framelabel1=n framelabel2=n framelabel3=n \
        labelfat=1 labelsz=8 o1num=-0.3 d1num=0.3 n1tic=3 o3num=-0.3 d3num=0.3 n3tic=3 > junk1.v
	sfcubeplot < 2.T point1=0.80 point2=0.52 screenratio=1.6 frame1=45 frame2=10 frame3=9 title=" " color=j \
        wantaxis1=n label2="hx" unit2="km" label3="hy" unit3="km" framelabel1=n framelabel2=n framelabel3=n \
        labelfat=1 labelsz=8 o1num=-0.3 d1num=0.3 n1tic=3 o3num=-0.3 d3num=0.3 n3tic=3 > junk2.v
	sfcubeplot < 3.T point1=0.80 point2=0.52 screenratio=1.6 frame1=39 frame2=10 frame3=9 title=" " color=j \
        wantaxis1=n label2="hx" unit2="km" label3="hy" unit3="km" framelabel1=n framelabel2=n framelabel3=n \
        labelfat=1 labelsz=8 o1num=-0.3 d1num=0.3 n1tic=3 o3num=-0.3 d3num=0.3 n3tic=3 > junk3.v
	vppen < junk1.v vpstyle=n align=lu | vppen xshift=1.5   yshift=0 > 1.v
	vppen < junk2.v vpstyle=n align=lu | vppen xshift=7.49  yshift=0 > 2.v
	vppen < junk3.v vpstyle=n align=lu | vppen xshift=12.5  yshift=0 > 3.v
	vppen 1.v 2.v 3.v erase=n scale=0.3 > $@
	Rm junk*.H 1.T 2.T 3.T
	rm -f junk1.v junk2.v junk3.v 1.v 2.v 3.v




${RESDIR}/bpgom3d-born-visi.v:
	Window3d < ${BPGOM3D_IMAG}/bpgom3d-born-visibility.H n1=1 min1=66 | Transp > junk1.H
	${BINDIR}/normalize2d.x < junk1.H > junk2.H
	echo o1=-328 d1=32.8 o2=-328 d2=32.8 >> junk2.H
	Grey < junk2.H pclip=100 color= allpos=y title=" " label1="psy (\F10 m\F3 s/m)" label2="psx (\F10 m\F3 s/m)" labelrot=n \
        wantscalebar=y barlabel=" " bclip=0. eclip=1.001 newclip=1 screenratio=0.85 labelfat=2 labelsz=8 out=$@ \
        grid=1 o1num=-300 o2num=-300 > /dev/null
	#Overlay < junk1.v x=3.4    y=3.1 type=box width=4.95 height=4.05 fat=3 out=junk2.v > /dev/null
	#vp_Overlay junk1.v junk2.v > $@
	#Rm junk1.H junk2.H
	#rm -f junk1.v junk2.v 

${RESDIR}/bpgom3d-born-tsou.v:
	Window3d < ${BPGOM3D_IMAG}/bpgom3d-tsou-born-planes-10-10.H | Grey eout=1 gainpanel=a pclip=100 > 1.T
	${BINDIR}/ft2m.x header=1.T axis1=n axis2=y axis3=y km=y
	echo o2=1.95 o3=-0.3 >> 1.T
	sfcubeplot < 1.T point1=0.63 point2=0.69 title=" " label1="Time" unit1="s" label2="Inline" unit2="km" label3="Crossline" unit3="km"\
        frame1=1000 frame2=132 frame3=60 \
        framelabel1=n framelabel2=n framelabel3=n labelfat=2 labelsz=6 > $@
	Rm 1.T
	rm -f junk1.v

${RESDIR}/bpgom3d-born-trec.v:
	Window3d < ${BPGOM3D_IMAG}/bpgom3d-trec-born-planes-10-10.H | Grey eout=1 gainpanel=a pclip=98 > 1.T
	${BINDIR}/ft2m.x header=1.T axis1=n axis2=y axis3=y km=y
	echo o2=1.95 o3=-0.3 >> 1.T
	sfcubeplot < 1.T point1=0.63 point2=0.69 title=" " label1="Time" unit1="s" label2="Inline" unit2="km" label3="Crossline" unit3="km"\
        frame1=1000 frame2=132 frame3=60 \
        framelabel1=n framelabel2=n framelabel3=n labelfat=2 labelsz=6 > $@
	Rm 1.T
	rm -f junk1.v


# Reflectivity imaging
${RESDIR}/lsm3d-target.v:
	Window3d < ${BPGOM3D_IMAG}/bpgom3d-vmod-invt-full.H n3=584 min1=42000 max1=95000 min2=65000 max2=85000 > junk1.H
	${BINDIR}/ft2m.x header=junk1.H axis1=y axis2=y axis3=y km=y
	${BINDIR}/ft2m_bin.x < junk1.H km=y > junk2.H
	echo o1=0 o2=0 >> junk2.H
	Window3d < junk2.H | Transp reshape=2,3 > junk3.H
	sfbyte < junk3.H gainpanel=a allpos=y pclip=100 > j.T
	sfcubeplot < j.T bar=bar-vmod.rsf color=j title=" " label1="Depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        scalebar=y barlabel="Velocity" point1=0.64 point2=0.72 o1num=0 d1num=2 n1tic=8 n2tic=0\
        frame1=290 frame2=220 frame3=30 framelabel1=n framelabel2=n framelabel3=n barunit="km/s" labelfat=1 labelsz=6 > junk1.v
	Overlay < junk1.v x=4.95 y=3.6  type=box width=2.45 height=1.1  out=junk2.v fat=3 > /dev/null
	Overlay < junk1.v x=8.69 y=3.6  type=box width=1.22 height=1.1  out=junk3.v fat=3 > /dev/null
	Overlay < junk1.v x=4.95 y=6.23 type=box width=2.45 height=1.22 out=junk4.v fat=3 > /dev/null
	vp_Overlay junk1.v junk2.v junk3.v junk4.v > $@
	Rm junk1.H junk2.H junk3.H j.T
	rm -f junk*.v

${RESDIR}/lsm3d-target-grey.v:
	Window3d < ${BPGOM3D_IMAG}/bpgom3d-vmod-invt-full.H n3=584 min1=42000 max1=95000 min2=65000 max2=85000 > junk1.H
	${BINDIR}/ft2m.x header=junk1.H axis1=y axis2=y axis3=y km=y
	${BINDIR}/ft2m_bin.x < junk1.H km=y > junk2.H
	echo o1=0 o2=0 >> junk2.H
	Window3d < junk2.H | Transp reshape=2,3 > junk3.H
	sfbyte < junk3.H gainpanel=a allpos=y pclip=100 > j.T
	sfcubeplot < j.T bar=bar-vmod.rsf title=" " label1="Depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        scalebar=y barlabel="Velocity" point1=0.64 point2=0.72 o1num=0 d1num=2 n1tic=8 n2tic=0\
        frame1=290 frame2=220 frame3=30 framelabel1=n framelabel2=n framelabel3=n barunit="km/s" labelfat=1 labelsz=6 > junk1.v
	Overlay < junk1.v x=4.95 y=3.6  type=box width=2.45 height=1.1  out=junk2.v fat=3 color=0 > /dev/null
	Overlay < junk1.v x=8.69 y=3.6  type=box width=1.22 height=1.1  out=junk3.v fat=3 color=0 > /dev/null
	Overlay < junk1.v x=4.95 y=6.23 type=box width=2.45 height=1.22 out=junk4.v fat=3 color=0 > /dev/null
	vp_Overlay junk1.v junk2.v junk3.v junk4.v > $@
	Rm junk1.H junk2.H junk3.H j.T
	rm -f junk*.v


SPIKE1 = frame1=45 frame2=26 frame3=21
SPIKE2 = frame1=67 frame2=78 frame3=44

${RESDIR}/lsm3d-imag-spike%.v:
	Window3d < ${LSM3D_IMAG}/lsm3d-imag-spike-half.H | Transp reshape=2,3,5 | ${BINDIR}/normalize3d.x > junk1.H
	${BINDIR}/ft2m.x header=junk1.H axis1=y axis2=y axis3=y km=y
	echo o1=0 o2=7.0 o3=0.2 >> junk1.H
	Grey < junk1.H eout=1 gainpanel=a > 1.T
	sfcubeplot < 1.T point1=0.53 point2=0.67 title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        ${SPIKE$*} framelabel1=n framelabel2=n framelabel3=n labelfat=1 labelsz=6 screenratio=0.6 \
        o1num=7 o2num=0 o3num=0.5 d1num=0.5 d2num=0.5 d3num=0.5 n1tic=13 n2tic=5 n3tic=6 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.8 > $@
	Rm junk1.H 1.T
	rm -f junk*.v

${RESDIR}/lsm3d-imag-flat%.v:
	Window3d < ${LSM3D_IMAG}/lsm3d-imag-flat.H | Transp reshape=2,3,5 | ${BINDIR}/normalize3d.x > junk1.H
	${BINDIR}/ft2m.x header=junk1.H axis1=y axis2=y axis3=y km=y
	echo o1=0 o2=7.0 o3=0.2 >> junk1.H
	Grey < junk1.H eout=1 gainpanel=a > 1.T
	sfcubeplot < 1.T point1=0.53 point2=0.67 title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        ${SPIKE$*} framelabel1=n framelabel2=n framelabel3=n labelfat=1 labelsz=6 screenratio=0.6 \
        o1num=7 o2num=0 o3num=0.5 d1num=0.5 d2num=0.5 d3num=0.5 n1tic=13 n2tic=5 n3tic=6 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.8 > $@
	Rm junk1.H 1.T
	rm -f junk*.v

${RESDIR}/lsm3d-invt-noreg-flat%.v:
	Window3d < ${LSM3D_IMAG}/mmov-flat-noreg-100.H n4=1 f4=99 | Transp reshape=2,3,5 | ${BINDIR}/normalize3d.x > junk1.H
	${BINDIR}/ft2m.x header=junk1.H axis1=y axis2=y axis3=y km=y
	echo o1=0 o2=7.0 o3=0.2 >> junk1.H
	Grey < junk1.H eout=1 gainpanel=a > 1.T
	sfcubeplot < 1.T point1=0.53 point2=0.67 title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
	${SPIKE$*} framelabel1=n framelabel2=n framelabel3=n labelfat=1 labelsz=6 screenratio=0.6 \
        o1num=7 o2num=0 o3num=0.5 d1num=0.5 d2num=0.5 d3num=0.5 n1tic=13 n2tic=5 n3tic=6 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.8 > $@
	Rm junk1.H 1.T
	rm -f junk*.v

${RESDIR}/lsm3d-invt-reg-flat%.v:
	Window3d < ${LSM3D_IMAG}/mmov-flat.H n4=1 f4=99 | Transp reshape=2,3,5 | ${BINDIR}/normalize3d.x > junk1.H
	${BINDIR}/ft2m.x header=junk1.H axis1=y axis2=y axis3=y km=y
	echo o1=0 o2=7.0 o3=0.2 >> junk1.H
	Grey < junk1.H eout=1 gainpanel=a > 1.T
	sfcubeplot < 1.T point1=0.53 point2=0.67 title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        ${SPIKE$*} framelabel1=n framelabel2=n framelabel3=n labelfat=1 labelsz=6 screenratio=0.6 \
        o1num=7 o2num=0 o3num=0.5 d1num=0.5 d2num=0.5 d3num=0.5 n1tic=13 n2tic=5 n3tic=6 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.8 > $@
	Rm junk1.H 1.T
	rm -f junk*.v

${RESDIR}/lsm3d-invt-reg30-flat%.v:
	Window3d < ${LSM3D_IMAG}/mmov-dip.H n4=1 f4=99 | Transp reshape=2,3,5 | ${BINDIR}/normalize3d.x > junk1.H
	${BINDIR}/ft2m.x header=junk1.H axis1=y axis2=y axis3=y km=y
	echo o1=0 o2=7.0 o3=0.2 >> junk1.H
	Grey < junk1.H eout=1 gainpanel=a > 1.T
	sfcubeplot < 1.T point1=0.53 point2=0.67 title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        ${SPIKE$*} framelabel1=n framelabel2=n framelabel3=n labelfat=1 labelsz=6 screenratio=0.6 \
        o1num=7 o2num=0 o3num=0.5 d1num=0.5 d2num=0.5 d3num=0.5 n1tic=13 n2tic=5 n3tic=6 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.8 > $@
	Rm junk1.H 1.T
	rm -f junk*.v




HESS_FILTER1 = n1=1 min1=66436 n2=1 min2=73161 n3=1 min3=11940
HESS_FILTER2 = n1=1 min1=70373 n2=1 min2=70209 n3=1 min3=13620
HESS_DIAG_FILTER1 = frame2=9  frame3=46 frame1=16
HESS_DIAG_FILTER2 = frame2=33 frame3=28 frame1=44

bar-hess-diag.rsf:
	Window3d < ${LSM3D_IMAG}/lsm3d-hess-offd-target.H n4=1 min4=0 n5=1 min5=0 n6=1 min6=0 | Transp reshape=2,3,8 > junk1.H
	${BINDIR}/bar.x < junk1.H minval=0 maxval=1. > junk2.H
	sfbyte < junk2.H gainpanel=a pclip=100 allpos=y bar=$@ > j.T
	Rm junk*.H j.T


${RESDIR}/lsm3d-hess-filter%.v: bar-hess-diag.rsf
	Window3d < ${LSM3D_IMAG}/lsm3d-hess-offd-target.H ${HESS_FILTER$*} | Transp reshape=2,3,8 | ${BINDIR}/normalize3d.x > junk1.H
	Interp < junk1.H d1out=100 d2out=100 d3out=100 > junk2.H
	Window3d < ${LSM3D_IMAG}/lsm3d-hess-offd-target.H n4=1 min4=0 n5=1 min5=0 n6=1 min6=0 | Transp reshape=2,3,8 | ${BINDIR}/normalize3d.x > junk3.H
	${BINDIR}/ft2m.x header=junk1.H axis1=y axis2=y axis3=y km=y
	${BINDIR}/ft2m.x header=junk2.H axis1=y axis2=y axis3=y km=y
	${BINDIR}/ft2m.x header=junk3.H axis1=y axis2=y axis3=y km=y
	echo o1=0 o2=7.0 o3=0.2 >> junk3.H
	sfbyte < junk1.H gainpanel=a pclip=100 polarity=y > 1.T
	sfcubeplot < 1.T frame1=15 frame2=5 frame3=7 title=" " label1="Lag z" unit1="km" label2="Lag x" unit2="km" label3="Lag y" unit3="km" \
        framelabel1=n framelabel2=n framelabel3=n o1num=-0.2 d1num=0.2 o2num=-0.2 d2num=0.2 o3num=-0.2 d3num=0.2 n1tic=3 n2tic=3 n3tic=3\
        point1=0.43 point2=0.44 screenratio=0.9 labelfat=1 labelsz=8 cubelinecol=7 > 1.v
	Pad < junk2.H n1out=128 n2out=128 n3out=128 | Rtoc | Ft3d sign1=1 sign2=1 sign3=1 center1=1 center2=1 center3=1 | Cabs > junk0.H
	sfbyte < junk0.H gainpanel=a pclip=100 polarity=y allpos=y > 2.T
	sfcubeplot < 2.T frame1=30 frame2=55 frame3=62 title=" " label1="Kz" unit1="1/km" label2="Kx" unit2="1/km" label3="Ky" unit3="1/km" \
        framelabel1=n framelabel2=n framelabel3=n point1=0.5 0.73 point2=0.5 screenratio=0.9 labelfat=1 labelsz=8 cubelinecol=7 > 2.v
	sfbyte < junk3.H gainpanel=a allpos=y bar=bar.rsf pclip=100 > 3.T
	sfcubeplot < 3.T color=j point1=0.53 0.45 point2=0.67 bar=bar-hess-diag.rsf scalebar=y cubelinecol=7\
        title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        ${HESS_DIAG_FILTER$*} framelabel1=n framelabel2=n framelabel3=n labelfat=1 labelsz=5.33 screenratio=0.6 \
        o1num=7 o2num=0 d1num=0.5 d2num=0.5 o3num=0.5 d3num=0.5 n1tic=13 n2tic=5 n3tic=6 > 3.v
	vppen < 3.v vpstyle=n align=lu | vppen xshift=0.5  yshift=9 scale=1.51 > s1.v
	vppen < 1.v vpstyle=n align=lu | vppen xshift=0.5  yshift=0 > s2.v
	vppen < 2.v vpstyle=n align=lu | vppen xshift=11   yshift=0 > s3.v
	vppen s1.v s2.v s3.v erase=n scale=0.18 > $@


${RESDIR}/lsm3d-hess-grey-filter%.v: bar-hess-diag.rsf
	Window3d < ${LSM3D_IMAG}/lsm3d-hess-offd-target.H ${HESS_FILTER$*} | Transp reshape=2,3,8 | ${BINDIR}/normalize3d.x > junk1.H
	Interp < junk1.H d1out=100 d2out=100 d3out=100 > junk2.H
	Window3d < ${LSM3D_IMAG}/lsm3d-hess-offd-target.H n4=1 min4=0 n5=1 min5=0 n6=1 min6=0 | Transp reshape=2,3,8 | ${BINDIR}/normalize3d.x > junk3.H
	${BINDIR}/ft2m.x header=junk1.H axis1=y axis2=y axis3=y km=y
	${BINDIR}/ft2m.x header=junk2.H axis1=y axis2=y axis3=y km=y
	${BINDIR}/ft2m.x header=junk3.H axis1=y axis2=y axis3=y km=y
	echo o1=0 o2=7.0 o3=0.2 >> junk3.H
	sfbyte < junk1.H gainpanel=a pclip=100 polarity=y > 1.T
	sfcubeplot < 1.T frame1=15 frame2=5 frame3=7 title=" " label1="Lag z" unit1="km" label2="Lag x" unit2="km" label3="Lag y" unit3="km" \
        framelabel1=n framelabel2=n framelabel3=n o1num=-0.2 d1num=0.2 o2num=-0.2 d2num=0.2 o3num=-0.2 d3num=0.2 n1tic=3 n2tic=3 n3tic=3\
        point1=0.43 point2=0.44 screenratio=0.9 labelfat=1 labelsz=8 cubelinecol=7 > 1.v
	Pad < junk2.H n1out=128 n2out=128 n3out=128 | Rtoc | Ft3d sign1=1 sign2=1 sign3=1 center1=1 center2=1 center3=1 | Cabs > junk0.H
	sfbyte < junk0.H gainpanel=a pclip=100 allpos=y > 2.T
	sfcubeplot < 2.T frame1=30 frame2=55 frame3=62 title=" " label1="Kz" unit1="1/km" label2="Kx" unit2="1/km" label3="Ky" unit3="1/km" \
        framelabel1=n framelabel2=n framelabel3=n point1=0.5 0.73 point2=0.5 screenratio=0.9 labelfat=1 labelsz=8 cubelinecol=7 > 2.v
	sfbyte < junk3.H gainpanel=a allpos=y bar=bar.rsf pclip=100 > 3.T
	sfcubeplot < 3.T point1=0.53 0.45 point2=0.67 bar=bar-hess-diag.rsf scalebar=y cubelinecol=7\
        title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        ${HESS_DIAG_FILTER$*} framelabel1=n framelabel2=n framelabel3=n labelfat=1 labelsz=5.33 screenratio=0.6 \
        o1num=7 o2num=0 d1num=0.5 d2num=0.5 o3num=0.5 d3num=0.5 n1tic=13 n2tic=5 n3tic=6 > 3.v
	vppen < 3.v vpstyle=n align=lu | vppen xshift=0.5  yshift=9 scale=1.51 > s1.v
	vppen < 1.v vpstyle=n align=lu | vppen xshift=0.5  yshift=0 > s2.v
	vppen < 2.v vpstyle=n align=lu | vppen xshift=11   yshift=0 > s3.v
	vppen s1.v s2.v s3.v erase=n scale=0.18 > $@


${RESDIR}/lsm3d-refl-pick.v:
	Window3d < ${LSM3D_IMAG}/picks/lsm3d-refl-pick.H | Transp reshape=2,3,5 | ${BINDIR}/normalize3d.x > junk1.H
	${BINDIR}/ft2m.x header=junk1.H axis1=y axis2=y axis3=y km=y
	echo o1=0 o2=7.0 o3=0.2 >> junk1.H
	Grey < junk1.H eout=1 gainpanel=a > 1.T
	sfcubeplot < 1.T point1=0.53 point2=0.67 title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        frame1=55 frame2=72 frame3=19 framelabel1=n framelabel2=n framelabel3=n labelfat=1 labelsz=6 screenratio=0.6 \
        o1num=7 o2num=0 o3num=0.5 d1num=0.5 d2num=0.5 d3num=0.5 n1tic=13 n2tic=5 n3tic=6 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.8 > $@
	Rm junk1.H 1.T
	rm -f junk*.v

${RESDIR}/lsm3d-dipx-pick.v:
	${BINDIR}/a2p.x < ${LSM3D_IMAG}/lsm3d-dipx-pick.H > junk0.H
	Window3d < junk0.H | Transp reshape=2,3,8 > junk1.H
	${BINDIR}/ft2m.x header=junk1.H axis1=y axis2=y axis3=y km=y
	echo o1=0 o2=7.0 o3=0.2 >> junk1.H
	sfbyte < junk1.H gainpanel=a pclip=100 bar=dipx-bar.rsf > 1.T
	sfcubeplot < 1.T scalebar=y color=e bar=dipx-bar.rsf point1=0.53 point2=0.67 barlabel="Slope" title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        frame1=55 frame2=72 frame3=19 framelabel1=n framelabel2=n framelabel3=n labelfat=1 labelsz=6 screenratio=0.6 \
        o1num=7 o2num=0 o3num=0.5 d1num=0.5 d2num=0.5 d3num=0.5 n1tic=13 n2tic=5 n3tic=6 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.8 > $@
	Rm junk1.H 1.T
	rm -f junk*.v dipx-bar.rsf

${RESDIR}/lsm3d-dipx-pick-grey.v:
	${BINDIR}/a2p.x < ${LSM3D_IMAG}/lsm3d-dipx-pick.H > junk0.H
	Window3d < junk0.H | Transp reshape=2,3,8 > junk1.H
	${BINDIR}/ft2m.x header=junk1.H axis1=y axis2=y axis3=y km=y
	echo o1=0 o2=7.0 o3=0.2 >> junk1.H
	sfbyte < junk1.H gainpanel=a pclip=100 bar=dipx-bar.rsf > 1.T
	sfcubeplot < 1.T scalebar=y bar=dipx-bar.rsf point1=0.53 point2=0.67 barlabel="Slope" title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        frame1=55 frame2=72 frame3=19 framelabel1=n framelabel2=n framelabel3=n labelfat=1 labelsz=6 screenratio=0.6 \
        o1num=7 o2num=0 o3num=0.5 d1num=0.5 d2num=0.5 d3num=0.5 n1tic=13 n2tic=5 n3tic=6 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.8 > $@
	Rm junk1.H 1.T
	rm -f junk*.v dipx-bar.rsf


${RESDIR}/lsm3d-dipy-pick.v:
	${BINDIR}/a2p.x < ${LSM3D_IMAG}/lsm3d-dipy-pick.H > junk0.H
	Window3d < junk0.H | Transp reshape=2,3,8 > junk1.H
	${BINDIR}/ft2m.x header=junk1.H axis1=y axis2=y axis3=y km=y
	echo o1=0 o2=7.0 o3=0.2 >> junk1.H
	sfbyte < junk1.H gainpanel=a pclip=100 bar=dipy-bar.rsf > 1.T
	sfcubeplot < 1.T scalebar=y color=e bar=dipy-bar.rsf point1=0.53 point2=0.67 barlabel="Slope" title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        frame1=55 frame2=72 frame3=19 framelabel1=n framelabel2=n framelabel3=n labelfat=1 labelsz=6 screenratio=0.6 \
        o1num=7 o2num=0 o3num=0.5 d1num=0.5 d2num=0.5 d3num=0.5 n1tic=13 n2tic=5 n3tic=6 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.8 > $@
	Rm junk1.H 1.T
	rm -f junk*.v dipy-bar.rsf

${RESDIR}/lsm3d-dipy-pick-grey.v:
	${BINDIR}/a2p.x < ${LSM3D_IMAG}/lsm3d-dipy-pick.H > junk0.H
	Window3d < junk0.H | Transp reshape=2,3,8 > junk1.H
	${BINDIR}/ft2m.x header=junk1.H axis1=y axis2=y axis3=y km=y
	echo o1=0 o2=7.0 o3=0.2 >> junk1.H
	sfbyte < junk1.H gainpanel=a pclip=100 bar=dipy-bar.rsf > 1.T
	sfcubeplot < 1.T scalebar=y bar=dipy-bar.rsf point1=0.53 point2=0.67 barlabel="Slope" title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        frame1=55 frame2=72 frame3=19 framelabel1=n framelabel2=n framelabel3=n labelfat=1 labelsz=6 screenratio=0.6 \
        o1num=7 o2num=0 o3num=0.5 d1num=0.5 d2num=0.5 d3num=0.5 n1tic=13 n2tic=5 n3tic=6 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.8 > $@
	Rm junk1.H 1.T
	rm -f junk*.v dipy-bar.rsf




${RESDIR}/lsm3d-dipflt%.v:
	Window3d < ${LSM3D_IMAG}/lsm3d-dipflt$*.H | Transp reshape=2,3,8 > junk1.H
	${BINDIR}/ft2m.x header=junk1.H axis1=y axis2=y axis3=y km=y
	echo o1=0 o2=7.0 o3=0.2 >> junk1.H
	Grey < junk1.H eout=1 gainpanel=a pclip=100 > 1.T
	sfcubeplot < 1.T point1=0.53 point2=0.67 title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        frame1=54 frame2=72 frame3=19 framelabel1=n framelabel2=n framelabel3=n labelfat=1 labelsz=6 screenratio=0.6 \
        o1num=7 o2num=0 o3num=0.5 d1num=0.5 d2num=0.5 d3num=0.5 n1tic=13 n2tic=5 n3tic=6 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.8 > $@
	Rm junk1.H 1.T
	rm -f junk*.v


${RESDIR}/lsm3d-imag.v:
	Window3d < ${LSM3D_IMAG}/lsm3d-imag-target.H | Transp reshape=2,3,8 | ${BINDIR}/normalize3d.x > junk1.H
	${BINDIR}/ft2m.x header=junk1.H axis1=y axis2=y axis3=y km=y
	echo o1=0 o2=7.0 o3=0.2 >> junk1.H
	Grey < junk1.H eout=1 gainpanel=a pclip=97 > 1.T
	sfcubeplot < 1.T point1=0.53 point2=0.67 title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        frame1=67 frame2=67 frame3=10 framelabel1=n framelabel2=n framelabel3=n labelfat=1 labelsz=6 screenratio=0.6 \
        o1num=7 o2num=0 o3num=0.5 d1num=0.5 d2num=0.5 d3num=0.5 n1tic=13 n2tic=5 n3tic=6 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.8 > $@
	Rm junk1.H 1.T
	rm -f junk*.v

${RESDIR}/lsm3d-mask.v:
	Window3d < ${LSM3D_IMAG}/lsm3d-mask.H | Transp reshape=2,3,8 | ${BINDIR}/normalize3d.x > junk1.H
	${BINDIR}/ft2m.x header=junk1.H axis1=y axis2=y axis3=y km=y
	echo o1=0 o2=7.0 o3=0.2 >> junk1.H
	Grey < junk1.H eout=1 gainpanel=a clip=1 > 1.T
	sfcubeplot < 1.T point1=0.53 point2=0.67 title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        frame1=48 frame2=65 frame3=28 framelabel1=n framelabel2=n framelabel3=n labelfat=1 labelsz=6 screenratio=0.6 \
        o1num=7 o2num=0 o3num=0.5 d1num=0.5 d2num=0.5 d3num=0.5 n1tic=13 n2tic=5 n3tic=6 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.8 > $@
	Rm junk1.H 1.T
	rm -f junk*.v


LSM_LOC1 = frame1=67 frame2=67 frame3=10
LSM_LOC2 = frame1=52 frame2=39 frame3=32
LSM_CLIP = pclip=98
MIG_CLIP = pclip=98
NOREG_CLIP = pclip=98

#LSM_V: ${RESDIR}/lsm3d-a-invt%.v ${RESDIR}/lsm3d-b-invt%.v ${RESDIR}/lsm3d-c-invt%.v ${RESDIR}/lsm3d-noreg-invt%.v ${RESDIR}/lsm3d-imag-mask%.v

${RESDIR}/lsm3d-a-invt%.v:
	Math file1=${LSM3D_IMAG}/lsm3d-invt3.H file2=${LSM3D_IMAG}/lsm3d-mask-display.H exp="file1*file2" > junk0.H
	Window3d < junk0.H | Transp reshape=2,3,8 | ${BINDIR}/normalize3d.x > junk1.H
	${BINDIR}/ft2m.x header=junk1.H axis1=y axis2=y axis3=y km=y
	echo o1=0 o2=7.0 o3=0.2 >> junk1.H
	Grey < junk1.H eout=1 gainpanel=a ${LSM_CLIP} > 1.T
	sfcubeplot < 1.T point1=0.53 point2=0.67 title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        ${LSM_LOC$*} framelabel1=n framelabel2=n framelabel3=n labelfat=1 labelsz=6 screenratio=0.6 \
        o1num=7 o2num=0 o3num=0.5 d1num=0.5 d2num=0.5 d3num=0.5 n1tic=13 n2tic=5 n3tic=6 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.8 > $@
	Rm junk1.H 1.T
	rm -f junk*.v

${RESDIR}/lsm3d-b-invt%.v:
	Math file1=${LSM3D_IMAG}/lsm3d-invt1.H file2=${LSM3D_IMAG}/lsm3d-mask-display.H exp="file1*file2" > junk0.H
	Window3d < junk0.H | Transp reshape=2,3,8 | ${BINDIR}/normalize3d.x > junk1.H
	${BINDIR}/ft2m.x header=junk1.H axis1=y axis2=y axis3=y km=y
	echo o1=0 o2=7.0 o3=0.2 >> junk1.H
	Grey < junk1.H eout=1 gainpanel=a ${LSM_CLIP} > 1.T
	sfcubeplot < 1.T point1=0.53 point2=0.67 title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        ${LSM_LOC$*} framelabel1=n framelabel2=n framelabel3=n labelfat=1 labelsz=6 screenratio=0.6 \
        o1num=7 o2num=0 o3num=0.5 d1num=0.5 d2num=0.5 d3num=0.5 n1tic=13 n2tic=5 n3tic=6 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.8 > $@
	Rm junk1.H 1.T
	rm -f junk*.v

${RESDIR}/lsm3d-c-invt%.v:
	Math file1=${LSM3D_IMAG}/lsm3d-invt6.H file2=${LSM3D_IMAG}/lsm3d-mask-display.H exp="file1*file2" > junk0.H
	Window3d < junk0.H | Transp reshape=2,3,8 | ${BINDIR}/normalize3d.x > junk1.H
	${BINDIR}/ft2m.x header=junk1.H axis1=y axis2=y axis3=y km=y
	echo o1=0 o2=7.0 o3=0.2 >> junk1.H
	Grey < junk1.H eout=1 gainpanel=a ${LSM_CLIP} > 1.T
	sfcubeplot < 1.T point1=0.53 point2=0.67 title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        ${LSM_LOC$*} framelabel1=n framelabel2=n framelabel3=n labelfat=1 labelsz=6 screenratio=0.6 \
        o1num=7 o2num=0 o3num=0.5 d1num=0.5 d2num=0.5 d3num=0.5 n1tic=13 n2tic=5 n3tic=6 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.8 > $@
	Rm junk1.H 1.T
	rm -f junk*.v

${RESDIR}/lsm3d-noreg-invt%.v:
	Math file1=${LSM3D_IMAG}/lsm3d-noreg-invt.H file2=${LSM3D_IMAG}/lsm3d-mask-display.H exp="file1*file2" > junk0.H
	Window3d < junk0.H | Transp reshape=2,3,8 | ${BINDIR}/normalize3d.x > junk1.H
	${BINDIR}/ft2m.x header=junk1.H axis1=y axis2=y axis3=y km=y
	echo o1=0 o2=7.0 o3=0.2 >> junk1.H
	Grey < junk1.H eout=1 gainpanel=a ${NOREG_CLIP}  > 1.T
	sfcubeplot < 1.T point1=0.53 point2=0.67 title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        ${LSM_LOC$*} framelabel1=n framelabel2=n framelabel3=n labelfat=1 labelsz=6 screenratio=0.6 \
        o1num=7 o2num=0 o3num=0.5 d1num=0.5 d2num=0.5 d3num=0.5 n1tic=13 n2tic=5 n3tic=6 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.8 > $@
	Rm junk1.H 1.T
	rm -f junk*.v

${RESDIR}/lsm3d-damp1-invt%.v:
	Math file1=${LSM3D_IMAG}/lsm3d-damp-invt1.H file2=${LSM3D_IMAG}/lsm3d-mask-display.H exp="file1*file2" > junk0.H
	Window3d < junk0.H | Transp reshape=2,3,8 | ${BINDIR}/normalize3d.x > junk1.H
	${BINDIR}/ft2m.x header=junk1.H axis1=y axis2=y axis3=y km=y
	echo o1=0 o2=7.0 o3=0.2 >> junk1.H
	Grey < junk1.H eout=1 gainpanel=a ${NOREG_CLIP}  > 1.T
	sfcubeplot < 1.T point1=0.53 point2=0.67 title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        ${LSM_LOC$*} framelabel1=n framelabel2=n framelabel3=n labelfat=1 labelsz=6 screenratio=0.6 \
        o1num=7 o2num=0 o3num=0.5 d1num=0.5 d2num=0.5 d3num=0.5 n1tic=13 n2tic=5 n3tic=6 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.8 > $@
	Rm junk1.H 1.T
	rm -f junk*.v

${RESDIR}/lsm3d-damp2-invt%.v:
	Math file1=${LSM3D_IMAG}/lsm3d-damp-invt2.H file2=${LSM3D_IMAG}/lsm3d-mask-display.H exp="file1*file2" > junk0.H
	Window3d < junk0.H | Transp reshape=2,3,8 | ${BINDIR}/normalize3d.x > junk1.H
	${BINDIR}/ft2m.x header=junk1.H axis1=y axis2=y axis3=y km=y
	echo o1=0 o2=7.0 o3=0.2 >> junk1.H
	Grey < junk1.H eout=1 gainpanel=a ${NOREG_CLIP}  > 1.T
	sfcubeplot < 1.T point1=0.53 point2=0.67 title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        ${LSM_LOC$*} framelabel1=n framelabel2=n framelabel3=n labelfat=1 labelsz=6 screenratio=0.6 \
        o1num=7 o2num=0 o3num=0.5 d1num=0.5 d2num=0.5 d3num=0.5 n1tic=13 n2tic=5 n3tic=6 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.8 > $@
	Rm junk1.H 1.T
	rm -f junk*.v

${RESDIR}/lsm3d-damp3-invt%.v:
	Math file1=${LSM3D_IMAG}/lsm3d-damp-invt3.H file2=${LSM3D_IMAG}/lsm3d-mask-display.H exp="file1*file2" > junk0.H
	Window3d < junk0.H | Transp reshape=2,3,8 | ${BINDIR}/normalize3d.x > junk1.H
	${BINDIR}/ft2m.x header=junk1.H axis1=y axis2=y axis3=y km=y
	echo o1=0 o2=7.0 o3=0.2 >> junk1.H
	Grey < junk1.H eout=1 gainpanel=a ${NOREG_CLIP}  > 1.T
	sfcubeplot < 1.T point1=0.53 point2=0.67 title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        ${LSM_LOC$*} framelabel1=n framelabel2=n framelabel3=n labelfat=1 labelsz=6 screenratio=0.6 \
        o1num=7 o2num=0 o3num=0.5 d1num=0.5 d2num=0.5 d3num=0.5 n1tic=13 n2tic=5 n3tic=6 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.8 > $@
	Rm junk1.H 1.T
	rm -f junk*.v

${RESDIR}/lsm3d-damp4-invt%.v:
	Math file1=${LSM3D_IMAG}/lsm3d-damp-invt4.H file2=${LSM3D_IMAG}/lsm3d-mask-display.H exp="file1*file2" > junk0.H
	Window3d < junk0.H | Transp reshape=2,3,8 | ${BINDIR}/normalize3d.x > junk1.H
	${BINDIR}/ft2m.x header=junk1.H axis1=y axis2=y axis3=y km=y
	echo o1=0 o2=7.0 o3=0.2 >> junk1.H
	Grey < junk1.H eout=1 gainpanel=a ${NOREG_CLIP}  > 1.T
	sfcubeplot < 1.T point1=0.53 point2=0.67 title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        ${LSM_LOC$*} framelabel1=n framelabel2=n framelabel3=n labelfat=1 labelsz=6 screenratio=0.6 \
        o1num=7 o2num=0 o3num=0.5 d1num=0.5 d2num=0.5 d3num=0.5 n1tic=13 n2tic=5 n3tic=6 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.8 > $@
	Rm junk1.H 1.T
	rm -f junk*.v

${RESDIR}/lsm3d-damp5-invt%.v:
	Math file1=${LSM3D_IMAG}/lsm3d-damp-invt5.H file2=${LSM3D_IMAG}/lsm3d-mask-display.H exp="file1*file2" > junk0.H
	Window3d < junk0.H | Transp reshape=2,3,8 | ${BINDIR}/normalize3d.x > junk1.H
	${BINDIR}/ft2m.x header=junk1.H axis1=y axis2=y axis3=y km=y
	echo o1=0 o2=7.0 o3=0.2 >> junk1.H
	Grey < junk1.H eout=1 gainpanel=a ${NOREG_CLIP}  > 1.T
	sfcubeplot < 1.T point1=0.53 point2=0.67 title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        ${LSM_LOC$*} framelabel1=n framelabel2=n framelabel3=n labelfat=1 labelsz=6 screenratio=0.6 \
        o1num=7 o2num=0 o3num=0.5 d1num=0.5 d2num=0.5 d3num=0.5 n1tic=13 n2tic=5 n3tic=6 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.8 > $@
	Rm junk1.H 1.T
	rm -f junk*.v





${RESDIR}/lsm3d-imag-mask%.v:
	${BINDIR}/normalize3d.x < ${LSM3D_IMAG}/lsm3d-imag-target.H > junk00.H
	Math file1=junk00.H file2=${LSM3D_IMAG}/lsm3d-mask-display.H exp="file1*file2" > junk0.H
	Window3d < junk0.H | Transp reshape=2,3,8 | ${BINDIR}/normalize3d.x > junk1.H
	${BINDIR}/ft2m.x header=junk1.H axis1=y axis2=y axis3=y km=y
	echo o1=0 o2=7.0 o3=0.2 >> junk1.H
	Grey < junk1.H eout=1 gainpanel=a ${MIG_CLIP}  > 1.T
	sfcubeplot < 1.T point1=0.53 point2=0.67 title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        ${LSM_LOC$*} framelabel1=n framelabel2=n framelabel3=n labelfat=1 labelsz=6 screenratio=0.6 \
        o1num=7 o2num=0 o3num=0.5 d1num=0.5 d2num=0.5 d3num=0.5 n1tic=13 n2tic=5 n3tic=6 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.8 > $@
	Rm junk1.H 1.T
	rm -f junk*.v

${RESDIR}/lsm3d-imag-agc%.v:
	${BINDIR}/normalize3d.x < ${LSM3D_IMAG}/lsm3d-imag-target.H > junk00.H
	Math file1=junk00.H file2=${LSM3D_IMAG}/lsm3d-mask-display.H exp="file1*file2" > junk0.H
	Window3d < junk0.H | Transp reshape=2,3,8 | Agc window=80 | ${BINDIR}/normalize3d.x > junk1.H
	${BINDIR}/ft2m.x header=junk1.H axis1=y axis2=y axis3=y km=y
	echo o1=0 o2=7.0 o3=0.2 >> junk1.H
	Grey < junk1.H eout=1 gainpanel=a ${MIG_CLIP}  > 1.T
	sfcubeplot < 1.T point1=0.53 point2=0.67 title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        ${LSM_LOC$*} framelabel1=n framelabel2=n framelabel3=n labelfat=1 labelsz=6 screenratio=0.6 \
        o1num=7 o2num=0 o3num=0.5 d1num=0.5 d2num=0.5 d3num=0.5 n1tic=13 n2tic=5 n3tic=6 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.8 | tube
	Rm junk1.H 1.T
	rm -f junk*.v


${RESDIR}/lsm3d-test-invt%.v:
	Window3d < ${LSM3D_IMAG}/mmov.lsm3d-invt3.H n4=1 f4=99 > junk00.H
	Math file1=junk00.H file2=${LSM3D_IMAG}/lsm3d-mask.H exp="file1*file2" > junk0.H
	Window3d < junk0.H | Transp reshape=2,3,8 | Agc window=40 | ${BINDIR}/normalize3d.x > junk1.H
	${BINDIR}/ft2m.x header=junk1.H axis1=y axis2=y axis3=y km=y
	echo o1=0 o2=7.0 o3=0.2 >> junk1.H
	Grey < junk1.H eout=1 gainpanel=a ${MIG_CLIP}  > 1.T
	sfcubeplot < 1.T point1=0.53 point2=0.67 title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        ${LSM_LOC$*} framelabel1=n framelabel2=n framelabel3=n labelfat=1 labelsz=6 screenratio=0.6 \
        o1num=7 o2num=0 o3num=0.5 d1num=0.5 d2num=0.5 d3num=0.5 n1tic=13 n2tic=5 n3tic=6 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.8 | tube
	Rm junk1.H 1.T
	rm -f junk*.v





RESD_CLIP = clip=0.1

bar-lsm.rsf:
	${BINDIR}/bar.x < ${LSM3D_IMAG}/lsm3d-imag-target.H  minval=-1 maxval=1 > junk1.H
	sfbyte < junk1.H gainpanel=a pclip=100 allpos=y bar=$@ > j.T
	Rm junk*.H j.T

${RESDIR}/lsm3d-a-resd.v:
	Window3d < ${LSM3D_IMAG}/rmov.lsm3d-invt3.H n4=1 f4=99 > junk0.H
	Math file1=junk0.H exp="file1/0.284504" > junk00.H
	Window3d < junk00.H | ${BINDIR}/fix.x | Transp reshape=2,3,8 > junk1.H
	${BINDIR}/ft2m.x header=junk1.H axis1=y axis2=y axis3=y km=y
	echo o1=0 o2=7.0 o3=0.2 >> junk1.H
	sfbyte < junk1.H gainpanel=a ${RESD_CLIP} > 1.T
	sfcubeplot < 1.T point1=0.53 point2=0.67 title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        ${LSM_LOC1} framelabel1=n framelabel2=n framelabel3=n labelfat=1 labelsz=6 screenratio=0.6 \
        o1num=7 o2num=0 o3num=0.5 d1num=0.5 d2num=0.5 d3num=0.5 n1tic=13 n2tic=5 n3tic=6 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.8 > $@
	Rm junk0.H junk1.H 1.T
	rm -f junk*.v

${RESDIR}/lsm3d-b-resd.v: 
	Window3d < ${LSM3D_IMAG}/rmov.lsm3d-invt1.H n4=1 f4=99 > junk0.H
	Math file1=junk0.H exp="file1/0.284504" > junk00.H
	Window3d < junk00.H | ${BINDIR}/fix.x | Transp reshape=2,3,8 > junk1.H
	${BINDIR}/ft2m.x header=junk1.H axis1=y axis2=y axis3=y km=y
	echo o1=0 o2=7.0 o3=0.2 >> junk1.H
	sfbyte < junk1.H gainpanel=a ${RESD_CLIP} > 1.T
	sfcubeplot < 1.T point1=0.53 point2=0.67 title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        ${LSM_LOC1} framelabel1=n framelabel2=n framelabel3=n labelfat=1 labelsz=6 screenratio=0.6 \
        o1num=7 o2num=0 o3num=0.5 d1num=0.5 d2num=0.5 d3num=0.5 n1tic=13 n2tic=5 n3tic=6 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.8 > $@
	Rm junk0.H junk1.H 1.T
	rm -f junk*.v

${RESDIR}/lsm3d-c-resd.v:
	Window3d < ${LSM3D_IMAG}/rmov.lsm3d-invt6.H n4=1 f4=99 > junk0.H
	Math file1=junk0.H exp="file1/0.284504" > junk00.H
	Window3d < junk00.H | ${BINDIR}/fix.x | Transp reshape=2,3,8 > junk1.H
	${BINDIR}/ft2m.x header=junk1.H axis1=y axis2=y axis3=y km=y
	echo o1=0 o2=7.0 o3=0.2 >> junk1.H
	sfbyte < junk1.H gainpanel=a ${RESD_CLIP} > 1.T
	sfcubeplot < 1.T point1=0.53 point2=0.67 title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        ${LSM_LOC1} framelabel1=n framelabel2=n framelabel3=n labelfat=1 labelsz=6 screenratio=0.6 \
        o1num=7 o2num=0 o3num=0.5 d1num=0.5 d2num=0.5 d3num=0.5 n1tic=13 n2tic=5 n3tic=6 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.8 > $@
	Rm junk0.H junk1.H 1.T
	rm -f junk*.v

${RESDIR}/lsm3d-noreg-resd.v: 
	Window3d < ${LSM3D_IMAG}/rmov.lsm3d-noreg-invt.H n4=1 f4=99 > junk0.H
	Math file1=junk0.H exp="file1/0.284504" > junk00.H
	Window3d < junk00.H | ${BINDIR}/fix.x | Transp reshape=2,3,8 > junk1.H
	${BINDIR}/ft2m.x header=junk1.H axis1=y axis2=y axis3=y km=y
	echo o1=0 o2=7.0 o3=0.2 >> junk1.H
	sfbyte < junk1.H gainpanel=a ${RESD_CLIP} > 1.T
	sfcubeplot < 1.T point1=0.53 point2=0.67 title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        ${LSM_LOC1} framelabel1=n framelabel2=n framelabel3=n labelfat=1 labelsz=6 screenratio=0.6 \
        o1num=7 o2num=0 o3num=0.5 d1num=0.5 d2num=0.5 d3num=0.5 n1tic=13 n2tic=5 n3tic=6 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.8 > $@
	Rm junk0.H junk1.H 1.T
	rm -f junk*.v

${RESDIR}/lsm3d-damp-resd%.v:
	Window3d < ${LSM3D_IMAG}/rmov.lsm3d-damp-invt$*.H n4=1 f4=99 > junk0.H
	Math file1=junk0.H exp="file1/0.284504" > junk00.H
	Window3d < junk00.H | ${BINDIR}/fix.x | Transp reshape=2,3,8 > junk1.H
	${BINDIR}/ft2m.x header=junk1.H axis1=y axis2=y axis3=y km=y
	echo o1=0 o2=7.0 o3=0.2 >> junk1.H
	sfbyte < junk1.H gainpanel=a ${RESD_CLIP} > 1.T
	sfcubeplot < 1.T point1=0.53 point2=0.67 title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        ${LSM_LOC1} framelabel1=n framelabel2=n framelabel3=n labelfat=1 labelsz=6 screenratio=0.6 \
        o1num=7 o2num=0 o3num=0.5 d1num=0.5 d2num=0.5 d3num=0.5 n1tic=13 n2tic=5 n3tic=6 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.8 > $@
	Rm junk0.H junk1.H 1.T
	rm -f junk*.v



${RESDIR}/rd0.v:
	cp ${RESDIR}/lsm3d-noreg-resd.v $@

${RESDIR}/rd1.v:
	vp_annotate < ${RESDIR}/lsm3d-damp-resd4.v batch=y text=Txt/text_damp4 > $@

${RESDIR}/rd2.v:
	vp_annotate < ${RESDIR}/lsm3d-damp-resd3.v batch=y text=Txt/text_damp3 > $@

${RESDIR}/rd3.v:
	vp_annotate < Fig/lsm3d-a-resd.v batch=y text=Txt/text_a > $@

${RESDIR}/rd4.v:
	vp_annotate < Fig/lsm3d-b-resd.v batch=y text=Txt/text_b > $@

${RESDIR}/rd5.v:
	vp_annotate < Fig/lsm3d-c-resd.v batch=y text=Txt/text_c > $@

${RESDIR}/lsm3d-resd-iter0.v:
	Window3d < ${LSM3D_IMAG}/lsm3d-resd-iter0.H | Transp reshape=2,3,8 | ${BINDIR}/normalize3d.x > junk1.H
	${BINDIR}/ft2m.x header=junk1.H axis1=y axis2=y axis3=y km=y
	echo o1=0 o2=7.0 o3=0.2 >> junk1.H
	Grey < junk1.H eout=1 gainpanel=a clip=1  > 1.T
	sfcubeplot < 1.T point1=0.53 point2=0.67 title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        ${LSM_LOC1} framelabel1=n framelabel2=n framelabel3=n labelfat=1 labelsz=6 screenratio=0.6 \
        o1num=7 o2num=0 o3num=0.5 d1num=0.5 d2num=0.5 d3num=0.5 n1tic=13 n2tic=5 n3tic=6 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.8 | tube
	Rm junk1.H 1.T
	rm -f junk*.v



${RESDIR}/lsm3d-fobj.v:
	Math file1=${LSM3D_IMAG}/lsm3d-resd-iter0.H exp="file1*file1" > junk00.H
	Stack3d < junk00.H | Stack3d | Transp | Stack3d > junk01.H
	echo o1=0 d1=1 >> junk01.H
	Cat junk01.H ${LSM3D_IMAG}/fobj.lsm3d-noreg-invt.H axis=1 | ${BINDIR}/normalize1d.x > junk1.H
	Cat junk01.H ${LSM3D_IMAG}/fobj.lsm3d-invt3.H      axis=1 | ${BINDIR}/normalize1d.x > junk2.H
	Cat junk01.H ${LSM3D_IMAG}/fobj.lsm3d-invt1.H      axis=1 | ${BINDIR}/normalize1d.x > junk3.H
	Cat junk01.H ${LSM3D_IMAG}/fobj.lsm3d-invt6.H      axis=1 | ${BINDIR}/normalize1d.x > junk4.H
	Cat junk01.H ${LSM3D_IMAG}/fobj.lsm3d-damp-invt4.H axis=1 | ${BINDIR}/normalize1d.x > junk5.H
	Cat junk01.H ${LSM3D_IMAG}/fobj.lsm3d-damp-invt3.H axis=1 | ${BINDIR}/normalize1d.x > junk6.H
	Cat junk1.H junk5.H junk6.H junk2.H junk3.H junk4.H axis=2 > junk0.H
	Graph < junk0.H legend=1 legendloc=tr legendtitle=" " legendsz=10 legendfat=1 plotfat=7,7,7,7,7,7 \
        title=" " label2="Image misfit" label1="Number of iterations" plotcol=7,5,6,3,4,2\
        curvelabel="No regularization":"Damping \F10 e=0.02":"Damping \F10 e=0.05":"Weak dip filtering":"Moderate dip filtering":"Strong dip filtering" labelrot=n \
        screenratio=0.4 labelsz=10 labelfat=1 out=1.v > /dev/null
	vppen < 1.v vpstyle=n scale=0.5 > $@
	Rm junk*.H
	rm -f 1.v





${RESDIR}/lsm3d-imag-cpst.v:
	Window3d < ${LSM3D_IMAG}/lsm3d-imag-cpst.H | Transp reshape=2,3,8 | ${BINDIR}/normalize3d.x > junk1.H
	${BINDIR}/ft2m.x header=junk1.H axis1=y axis2=y axis3=y km=y
	echo o1=0 o2=7.0 o3=0.2 >> junk1.H
	Grey < junk1.H eout=1 gainpanel=a pclip=92 > 1.T
	sfcubeplot < 1.T point1=0.45 point2=0.67 title=" " label1="Relative depth" unit1="km" label2="Inline" unit2="km" label3="Crossline" unit3="km" \
        frame1=38 frame2=65 frame3=41 framelabel1=n framelabel2=n framelabel3=n labelfat=1 labelsz=8 screenratio=0.48 \
        o1num=7 o2num=0 o3num=0.5 d1num=0.5 d2num=0.5 d3num=0.5 n1tic=13 n2tic=5 n3tic=6 > junk1.v
	vppen < junk1.v vpstyle=n scale=0.6 > $@
	Rm junk1.H 1.T
	rm -f junk*.v




${RESDIR}/bpgom3d-gmsk-target.ps:
	${SEP}/bin/pstexpen ${RESDIR}/bpgom3d-gmsk-target.v $@ color=yes fat=1 fatmult=1.5 invras=n force=

${RESDIR}/bpgom3d-grad-raw-target-iter%.ps:
	${SEP}/bin/pstexpen ${RESDIR}/bpgom3d-grad-raw-target-iter$*.v $@ color=yes fat=1 fatmult=1.5 invras=n force= 

${RESDIR}/bpgom3d-grad-smt-target-iter%.ps:
	${SEP}/bin/pstexpen ${RESDIR}/bpgom3d-grad-smt-target-iter$*.v $@ color=yes fat=1 fatmult=1.5 invras=n force= 

${RESDIR}/bpgom3d-invt-target-stage%.ps:
	${SEP}/bin/pstexpen ${RESDIR}/bpgom3d-invt-target-stage$*.v $@ color=yes fat=1 fatmult=1.5 invras=n force=

${RESDIR}/bpgom3d-bhes%.ps:
	${SEP}/bin/pstexpen ${RESDIR}/bpgom3d-bhes$*.v $@ color=yes fat=1 fatmult=1.5 invras=n force=

${RESDIR}/bpgom3d-hess%.ps:
	${SEP}/bin/pstexpen ${RESDIR}/bpgom3d-hess$*.v $@ color=yes fat=1 fatmult=1.5 invras=n force=

${RESDIR}/bpgom3d-cig3d-bhes.ps:
	${SEP}/bin/pstexpen ${RESDIR}/bpgom3d-cig3d-bhes.v $@ color=yes fat=1 fatmult=1.5 invras=n force=

${RESDIR}/lsm3d-target.ps:
	${SEP}/bin/pstexpen ${RESDIR}/lsm3d-target.v $@ color=yes fat=1 fatmult=1.5 invras=n force=

${RESDIR}/lsm3d-hess-filter%.ps:
	${SEP}/bin/pstexpen ${RESDIR}/lsm3d-hess-filter$*.v $@ color=yes fat=1 fatmult=1.5 invras=n force=

${RESDIR}/lsm3d-dipx-pick.ps:
	${SEP}/bin/pstexpen ${RESDIR}/lsm3d-dipx-pick.v $@ color=yes fat=1 fatmult=1.5 invras=n force=

${RESDIR}/lsm3d-dipy-pick.ps:
	${SEP}/bin/pstexpen ${RESDIR}/lsm3d-dipy-pick.v $@ color=yes fat=1 fatmult=1.5 invras=n force=

${RESDIR}/lsm3d-mask.ps:
	${SEP}/bin/pstexpen ${RESDIR}/lsm3d-mask.v $@ color=yes fat=1 fatmult=1.5 invras=n force=

${RESDIR}/lsm3d-fobj.ps:
	${SEP}/bin/pstexpen ${RESDIR}/lsm3d-fobj.v $@ color=yes fat=1 fatmult=1.5 invras=n force=

clean:
	rm -f *.o *.mod .make.dependencies.*

include ${SEPINC}/SEP.bottom
