# RTMVEL0akefile
# RTM imaging for velocity model 0
# for vel0.H file
# This -akefile is to test the formation of noises. Thus I will compare RTM using true/smooth velocity model from muted/unmuted data

# Figures
RTMVEL0RESULTSER = rtm1.vel0.mult.fullmuted rtm1.vel0.mult.fullmuted.lowcut \
		rtm1.truevel0.mult.fullmuted rtm1.truevel0.mult.fullmuted.lowcut \
		rtm1.vel0.nomult.fullmuted rtm1.vel0.nomult.fullmuted.lowcut \
		rtm1.truevel0.nomult.fullmuted rtm1.truevel0.nomult.fullmuted.lowcut
RTMVEL0RESULTSCR = 
RTMVEL0RESULTSNR = 
RTMVEL0FIG_pdf = rtm1.vel0.mult.fullmuted.pdf rtm1.vel0.mult.fullmuted.lowcut.pdf \
		 rtm1.truevel0.mult.fullmuted.pdf rtm1.truevel0.mult.fullmuted.lowcut.pdf \
		 rtm1.vel0.nomult.fullmuted.pdf rtm1.vel0.nomult.fullmuted.lowcut.pdf \
		 rtm1.truevel0.nomult.fullmuted.pdf rtm1.truevel0.nomult.fullmuted.lowcut.pdf

################################
## compare RTM images here
## - muted/unmuted data
## - true/smooth migration velocity
## - whether to apply lowcut filter
################################
compare.rtm%.vel0.mult: vel0.v vel0.mig.v rtm%.truevel0.mult.v rtm%.truevel0.mult.fullmuted.v rtm%.vel0.mult.v rtm%.vel0.mult.fullmuted.v rtm%.truevel0.mult.lowcut.v rtm%.truevel0.mult.fullmuted.lowcut.v rtm%.vel0.mult.lowcut.v rtm%.vel0.mult.fullmuted.lowcut.v
	tube $R/vel0.v $R/vel0.mig.v $R/rtm$*.truevel0.mult.v $R/rtm$*.truevel0.mult.fullmuted.v $R/rtm$*.vel0.mult.v $R/rtm$*.vel0.mult.fullmuted.v $R/rtm$*.truevel0.mult.lowcut.v $R/rtm$*.truevel0.mult.fullmuted.lowcut.v $R/rtm$*.vel0.mult.lowcut.v $R/rtm$*.vel0.mult.fullmuted.lowcut.v &

compare.rtm%.vel0.nomult: vel0.v vel0.mig.v rtm%.truevel0.nomult.v rtm%.truevel0.nomult.fullmuted.v rtm%.vel0.nomult.v rtm%.vel0.nomult.fullmuted.v
	tube $R/vel0.v $R/vel0.mig.v $R/rtm$*.truevel0.nomult.v $R/rtm$*.truevel0.nomult.fullmuted.v $R/rtm$*.vel0.nomult.v $R/rtm$*.vel0.nomult.fullmuted.v &


##########################
### generate data
## sh model2d.sh '$(m)' 'outputfile' 'modelprogram' 'velfile' 'soufile' 'parfile' '$(abstop)' 'jsx=4' 'sousnapoutfile' 'tmpdir=.'
##########################
## By default I will use $P/shot0.P file for the whole shots along 2D line.
##########
# with SRM
#wvfld%.mult.H:  vel%.H sou.HH ${P}/shot%.P ${B}/model2d_program.x 
#	<$< ${B}/model2d_program.x par=${P}/shot$*.P abstop=0 source=sou.HH $(dn) snap=$@
data%.vel0.mult.H souwvfld%.vel0.mult.H: vel0.H sou.H $P/shot%.vel0.P $B/model2d_program.x $(SH)/model2d.sh
	sh $(SH)/model2d.sh '$*.vel0.mult' 'data$*.vel0.mult.H' '$B/model2d_program.x' 'vel0.H' 'sou.H' '$P/shot$*.vel0.P' '0' '6' 'souwvfld$*.vel0.mult.H' './temp'
data%.vel0.mult: data%.vel0.mult.H
	<$< Transp plane=13 | Transp plane=34 | Grey title='data$* with SRM' $(Axis) label1="Time (s)" | Tube &
souwvfld%.vel0.mult: souwvfld%.vel0.mult.H
	<$< Grey title='souwvfld$* with SRM' $(Axis) | Tube &

##########
# without SRM
#wvfld%.nomult.H:  vel%.H sou.HH ${P}/shot%.P ${B}/model2d_program.x 
#	<$< ${B}/model2d_program.x par=${P}/shot$*.P abstop=1 source=sou.HH $(dn) snap=$@
data%.vel0.nomult.H souwvfld%.vel0.nomult.H:  vel0.H sou.H $P/shot%.vel0.P $B/model2d_program.x $(SH)/model2d.sh
	sh $(SH)/model2d.sh '$*.vel0.nomult' 'data$*.vel0.nomult.H' '$B/model2d_program.x' 'vel0.H' 'sou.H' '$P/shot$*.vel0.P' '1' '6' 'souwvfld$*.vel0.nomult.H' './temp'
data%.vel0.nomult: data%.vel0.nomult.H
	<$< Transp plane=13 | Transp plane=34 | Grey title='data$* without SRM' $(Axis) label1="Time (s)" | Tube &
souwvfld%.vel0.nomult: souwvfld%.vel0.nomult.H
	<$< Grey title='souwvfld$* without SRM' $(Axis) | Tube &

#################
## muting data
## sh mute2d.sh '$(m)' 'inputfile' 'outputfile' 'parfile' 'jsx=8' 'tmpdir=.'
#################
# with SRM
##full muted
data%.vel0.mult.fullmuted.H: data%.vel0.mult.H $P/fullmute.vel0.P $(SH)/mute2d.sh
	<$< Transp plane=13 |Transp plane=34 > fullmute$*.vel0.mult.in.H
	sh $(SH)/mute2d.sh '$*.vel0.mult.fullmuted' 'fullmute$*.vel0.mult.in.H' 'fullmute$*.vel0.mult.out.H' '$P/fullmute.vel0.P' '8' './temp'
	< fullmute$*.vel0.mult.out.H Transp plane=34 | Transp plane=13 > $@
	Rm fullmute$*.vel0.mult.in.H fullmute$*.vel0.mult.out.H
data%.vel0.mult.fullmuted: data%.vel0.mult.fullmuted.H
	<$< Transp plane=13 | Transp plane=34 | Grey title='fullmuted data$* with SRM' $(Axis) label1="Time (s)" | Tube &

# without SRM
##full muted
data%.vel0.nomult.fullmuted.H: data%.vel0.nomult.H $P/fullmute.vel0.P $(SH)/mute2d.sh
	<$< Transp plane=13 |Transp plane=34 > fullmute$*.vel0.nomult.in.H
	sh $(SH)/mute2d.sh '$*.vel0.nomult.fullmuted' 'fullmute$*.vel0.nomult.in.H' 'fullmute$*.vel0.nomult.out.H' '$P/fullmute.vel0.P' '8' './temp'
	< fullmute$*.vel0.nomult.out.H Transp plane=34 | Transp plane=13 > $@
	Rm fullmute$*.vel0.nomult.in.H fullmute$*.vel0.nomult.out.H
data%.vel0.nomult.fullmuted: data%.vel0.nomult.fullmuted.H
	<$< Transp plane=13 | Transp plane=34 | Grey title='fullmuted data$* with SRM' $(Axis) label1="Time (s)" | Tube &

##################################
### RTM images
##       sh rtm2d.sh '$(m)' 'in.H' 'out.H' 'rtmprogram' 'vel.mig.H' 'soufile' 'parfile' 'jsx=4' 'tmpdir=.' 'sousnapfile.out=sousnap${m}.H' 'recsnapfile.out=recsnap${m}.H' 'reflsnapfile.out=reflsnap${m}.H'
##################################
#########
## SRM
#########
######
## unmuted
######
## from unmuted data along with true migration velocity
rtm%.truevel0.mult.H rtm%.truevel0.mult.sousnap.H rtm%.truevel0.mult.recsnap.H rtm%.truevel0.mult.reflsnap.H: data%.vel0.mult.H ${B}/crtm2d_program.x vel0.H sou.H ${P}/rtm.vel0.P $(SH)/rtm2d.sh
	<$< Interp d2out=10 maxsize=1000 > data$*.truevel0.mult.interp.H
	sh ${SH}/rtm2d.sh '$*.truevel0.mult' 'data$*.truevel0.mult.interp.H' '$@' '${B}/crtm2d_program.x' 'vel0.H' 'sou.H' '${P}/rtm.vel0.P' '6' './temp' 'rtm$*.truevel0.mult.sousnap.H' 'rtm$*.truevel0.mult.recsnap.H' 'rtm$*.truevel0.mult.reflsnap.H'
	Rm data$*.truevel0.mult.interp.H
rtm%.truevel0.mult.v: rtm%.truevel0.mult.H
	<$< Grey title="An true-vel RTM image from unmuted data with SRM" $(Axis) screenratio=$(vel0screenratio) | Vppen vpstyle=n scale=$(vel0screenratio) $(dn) out=$R/$@
rtm%.truevel0.mult: rtm%.truevel0.mult.H
	Sep_cube rtm$*.truevel0.mult.sousnap.H rtm$*.truevel0.mult.recsnap.H rtm$*.truevel0.mult.reflsnap.H nviews=3 &

## from unmuted data along with smooth migration velocity
rtm%.vel0.mult.H rtm%.vel0.mult.sousnap.H rtm%.vel0.mult.recsnap.H rtm%.vel0.mult.reflsnap.H: data%.vel0.mult.H ${B}/crtm2d_program.x vel0.mig.H sou.H ${P}/rtm.vel0.P $(SH)/rtm2d.sh
	<$< Interp d2out=10 maxsize=1000 > data$*.vel0.mult.interp.H
	sh ${SH}/rtm2d.sh '$*.vel0.mult' 'data$*.vel0.mult.interp.H' '$@' '${B}/crtm2d_program.x' 'vel0.mig.H' 'sou.H' '${P}/rtm.vel0.P' '6' './temp' 'rtm$*.vel0.mult.sousnap.H' 'rtm$*.vel0.mult.recsnap.H' 'rtm$*.vel0.mult.reflsnap.H'
	Rm data$*.vel0.mult.interp.H
rtm%.vel0.mult.v: rtm%.vel0.mult.H
	<$< Grey title="An RTM image from unmuted data with SRM" $(Axis) screenratio=$(vel0screenratio) | Vppen vpstyle=n scale=$(vel0screenratio) $(dn) out=$R/$@
rtm%.vel0.mult: rtm%.vel0.mult.H
	Sep_cube rtm$*.vel0.mult.sousnap.H rtm$*.vel0.mult.recsnap.H rtm$*.vel0.mult.reflsnap.H nviews=3 &
#####
## full muted
#####
# from full muted data along with true migration velocity
rtm%.truevel0.mult.fullmuted.H rtm%.truevel0.mult.fullmuted.sousnap.H rtm%.truevel0.mult.fullmuted.recsnap.H rtm%.truevel0.mult.fullmuted.reflsnap.H: data%.vel0.mult.fullmuted.H ${B}/crtm2d_program.x vel0.H sou.H ${P}/rtm.vel0.P $(SH)/rtm2d.sh
	<$< Interp d2out=10 maxsize=1000 > data$*.truevel0.mult.fullmuted.interp.H
	sh ${SH}/rtm2d.sh '$*.truevel0.mult.fullmuted' 'data$*.truevel0.mult.fullmuted.interp.H' '$@' '${B}/crtm2d_program.x' 'vel0.H' 'sou.H' '${P}/rtm.vel0.P' '6' './temp' 'rtm$*.truevel0.mult.fullmuted.sousnap.H' 'rtm$*.truevel0.mult.fullmuted.recsnap.H' 'rtm$*.truevel0.mult.fullmuted.reflsnap.H'
	Rm data$*.truevel0.mult.fullmuted.interp.H
rtm%.truevel0.mult.fullmuted.v: rtm%.truevel0.mult.fullmuted.H
	<$< Grey title="An true-vel RTM image from fullmuted data with SRM" $(Axis) screenratio=$(vel0screenratio) | Vppen vpstyle=n scale=$(vel0screenratio) $(dn) out=$R/$@
rtm%.truevel0.mult.fullmuted: rtm%.truevel0.mult.fullmuted.H
	Sep_cube rtm$*.truevel0.mult.fullmuted.sousnap.H rtm$*.truevel0.mult.fullmuted.recsnap.H rtm$*.truevel0.mult.fullmuted.reflsnap.H nviews=3 &

rtm%.truevel0.mult.fullmuted.spectrum.v: rtm%.truevel0.mult.fullmuted.H $B/Fftw2D_program.x
	<$< Rtoc | $B/Fftw2D_program.x flag=1 dim=2 axis1=1 axis2=1 | Cabs | Grey title="The spectrum of an true-vel RTM image from fullmuted data with SRM" $(Axis) label1="Kz (/m)" label2="Kx (/m)" $(dn) out=$R/$@

# from full muted data along with smooth migration velocity
rtm%.vel0.mult.fullmuted.H rtm%.vel0.mult.fullmuted.sousnap.H rtm%.vel0.mult.fullmuted.recsnap.H rtm%.vel0.mult.fullmuted.reflsnap.H: data%.vel0.mult.fullmuted.H ${B}/crtm2d_program.x vel0.mig.H sou.H ${P}/rtm.vel0.P $(SH)/rtm2d.sh
	<$< Interp d2out=10 maxsize=1000 > data$*.vel0.mult.fullmuted.interp.H
	sh ${SH}/rtm2d.sh '$*.vel0.mult.fullmuted' 'data$*.vel0.mult.fullmuted.interp.H' '$@' '${B}/crtm2d_program.x' 'vel0.mig.H' 'sou.H' '${P}/rtm.vel0.P' '6' './temp' 'rtm$*.vel0.mult.fullmuted.sousnap.H' 'rtm$*.vel0.mult.fullmuted.recsnap.H' 'rtm$*.vel0.mult.fullmuted.reflsnap.H'
	Rm data$*.vel0.mult.fullmuted.interp.H
rtm%.vel0.mult.fullmuted.v: rtm%.vel0.mult.fullmuted.H
	<$< Grey title="An RTM image from fullmuted data with SRM" $(Axis) screenratio=$(vel0screenratio) | Vppen vpstyle=n scale=$(vel0screenratio) $(dn) out=$R/$@
rtm%.vel0.mult.fullmuted: rtm%.vel0.mult.fullmuted.H
	Sep_cube rtm$*.vel0.mult.fullmuted.sousnap.H rtm$*.vel0.mult.fullmuted.recsnap.H rtm$*.vel0.mult.fullmuted.reflsnap.H nviews=3 &

rtm%.vel0.mult.fullmuted.spectrum.v: rtm%.vel0.mult.fullmuted.H $B/Fftw2D_program.x
	<$< Rtoc | $B/Fftw2D_program.x flag=1 dim=2 axis1=1 axis2=1 | Cabs | Grey title="The spectrum of an RTM image from fullmuted data with SRM" $(Axis) label1="Kz (/m)" label2="Kx (/m)" $(dn) out=$R/$@

#########
## NO SRM
#########
######
## unmuted
######
## from unmuted data along with true migration velocity
rtm%.truevel0.nomult.H rtm%.truevel0.nomult.sousnap.H rtm%.truevel0.nomult.recsnap.H rtm%.truevel0.nomult.reflsnap.H: data%.vel0.nomult.H ${B}/crtm2d_program.x vel0.H sou.H ${P}/rtm.vel0.P $(SH)/rtm2d.sh
	<$< Interp d2out=10 maxsize=1000 > data$*.truevel0.nomult.interp.H
	sh ${SH}/rtm2d.sh '$*.truevel0.nomult' 'data$*.truevel0.nomult.interp.H' '$@' '${B}/crtm2d_program.x' 'vel0.H' 'sou.H' '${P}/rtm.vel0.P' '6' './temp' 'rtm$*.truevel0.nomult.sousnap.H' 'rtm$*.truevel0.nomult.recsnap.H' 'rtm$*.truevel0.nomult.reflsnap.H'
	Rm data$*.truevel0.nomult.interp.H
rtm%.truevel0.nomult.v: rtm%.truevel0.nomult.H
	<$< Grey title="An true-vel RTM image from unmuted data without SRM" $(Axis) screenratio=$(vel0screenratio) | Vppen vpstyle=n scale=$(vel0screenratio) $(dn) out=$R/$@
rtm%.truevel0.nomult: rtm%.truevel0.nomult.H
	Sep_cube rtm$*.truevel0.nomult.sousnap.H rtm$*.truevel0.nomult.recsnap.H rtm$*.truevel0.nomult.reflsnap.H nviews=3 &

## from unmuted data along with smooth migration velocity
rtm%.vel0.nomult.H rtm%.vel0.nomult.sousnap.H rtm%.vel0.nomult.recsnap.H rtm%.vel0.nomult.reflsnap.H: data%.vel0.nomult.H ${B}/crtm2d_program.x vel0.mig.H sou.H ${P}/rtm.vel0.P $(SH)/rtm2d.sh
	<$< Interp d2out=10 maxsize=1000 > data$*.vel0.nomult.interp.H
	sh ${SH}/rtm2d.sh '$*.vel0.nomult' 'data$*.vel0.nomult.interp.H' '$@' '${B}/crtm2d_program.x' 'vel0.mig.H' 'sou.H' '${P}/rtm.vel0.P' '6' './temp' 'rtm$*.vel0.nomult.sousnap.H' 'rtm$*.vel0.nomult.recsnap.H' 'rtm$*.vel0.nomult.reflsnap.H'
	Rm data$*.vel0.nomult.interp.H
rtm%.vel0.nomult.v: rtm%.vel0.nomult.H
	<$< Grey title="An RTM image from unmuted data without SRM" $(Axis) screenratio=$(vel0screenratio) | Vppen vpstyle=n scale=$(vel0screenratio) $(dn) out=$R/$@
rtm%.vel0.nomult: rtm%.vel0.nomult.H
	Sep_cube rtm$*.vel0.nomult.sousnap.H rtm$*.vel0.nomult.recsnap.H rtm$*.vel0.nomult.reflsnap.H nviews=3 &
#####
## full muted
#####
# from full muted data along with true migration velocity
rtm%.truevel0.nomult.fullmuted.H rtm%.truevel0.nomult.fullmuted.sousnap.H rtm%.truevel0.nomult.fullmuted.recsnap.H rtm%.truevel0.nomult.fullmuted.reflsnap.H: data%.vel0.nomult.fullmuted.H ${B}/crtm2d_program.x vel0.H sou.H ${P}/rtm.vel0.P $(SH)/rtm2d.sh
	<$< Interp d2out=10 maxsize=1000 > data$*.truevel0.nomult.fullmuted.interp.H
	sh ${SH}/rtm2d.sh '$*.truevel0.nomult.fullmuted' 'data$*.truevel0.nomult.fullmuted.interp.H' '$@' '${B}/crtm2d_program.x' 'vel0.H' 'sou.H' '${P}/rtm.vel0.P' '6' './temp' 'rtm$*.truevel0.nomult.fullmuted.sousnap.H' 'rtm$*.truevel0.nomult.fullmuted.recsnap.H' 'rtm$*.truevel0.nomult.fullmuted.reflsnap.H'
	Rm data$*.truevel0.nomult.fullmuted.interp.H
rtm%.truevel0.nomult.fullmuted.v: rtm%.truevel0.nomult.fullmuted.H
	<$< Grey title="An true-vel RTM image from fullmuted data without SRM" $(Axis) screenratio=$(vel0screenratio) | Vppen vpstyle=n scale=$(vel0screenratio) $(dn) out=$R/$@
rtm%.truevel0.nomult.fullmuted: rtm%.truevel0.nomult.fullmuted.H
	Sep_cube rtm$*.truevel0.nomult.fullmuted.sousnap.H rtm$*.truevel0.nomult.fullmuted.recsnap.H rtm$*.truevel0.nomult.fullmuted.reflsnap.H nviews=3 &

rtm%.truevel0.nomult.fullmuted.spectrum.v: rtm%.truevel0.nomult.fullmuted.H $B/Fftw2D_program.x
	<$< Rtoc | $B/Fftw2D_program.x flag=1 dim=2 axis1=1 axis2=1 | Cabs | Grey title="The spectrum of an true-vel RTM image from fullmuted data without SRM" $(Axis) label1="Kz (/m)" label2="Kx (/m)" $(dn) out=$R/$@

# from full muted data along with smooth migration velocity
rtm%.vel0.nomult.fullmuted.H rtm%.vel0.nomult.fullmuted.sousnap.H rtm%.vel0.nomult.fullmuted.recsnap.H rtm%.vel0.nomult.fullmuted.reflsnap.H: data%.vel0.nomult.fullmuted.H ${B}/crtm2d_program.x vel0.mig.H sou.H ${P}/rtm.vel0.P $(SH)/rtm2d.sh
	<$< Interp d2out=10 maxsize=1000 > data$*.vel0.nomult.fullmuted.interp.H
	sh ${SH}/rtm2d.sh '$*.vel0.nomult.fullmuted' 'data$*.vel0.nomult.fullmuted.interp.H' '$@' '${B}/crtm2d_program.x' 'vel0.mig.H' 'sou.H' '${P}/rtm.vel0.P' '6' './temp' 'rtm$*.vel0.nomult.fullmuted.sousnap.H' 'rtm$*.vel0.nomult.fullmuted.recsnap.H' 'rtm$*.vel0.nomult.fullmuted.reflsnap.H'
	Rm data$*.vel0.nomult.fullmuted.interp.H
rtm%.vel0.nomult.fullmuted.v: rtm%.vel0.nomult.fullmuted.H
	<$< Grey title="An RTM image from fullmuted data without SRM" $(Axis) screenratio=$(vel0screenratio) | Vppen vpstyle=n scale=$(vel0screenratio) $(dn) out=$R/$@
rtm%.vel0.nomult.fullmuted: rtm%.vel0.nomult.fullmuted.H
	Sep_cube rtm$*.vel0.nomult.fullmuted.sousnap.H rtm$*.vel0.nomult.fullmuted.recsnap.H rtm$*.vel0.nomult.fullmuted.reflsnap.H nviews=3 &

rtm%.vel0.nomult.fullmuted.spectrum.v: rtm%.vel0.nomult.fullmuted.H $B/Fftw2D_program.x
	<$< Rtoc | $B/Fftw2D_program.x flag=1 dim=2 axis1=1 axis2=1 | Cabs | Grey title="The spectrum of an RTM image from fullmuted data without SRM" $(Axis) label1="Kz (/m)" label2="Kx (/m)" $(dn) out=$R/$@

#################################
### RTM with SRM after lowcut filtered
#################################
## from unmuted data
# true vel
rtm%.truevel0.mult.lowcut.H: rtm%.truevel0.mult.H $B/Light1.x
	<$< $B/Light1.x >$@
rtm%.truevel0.mult.lowcut.v: rtm%.truevel0.mult.lowcut.H
	<$< Window n3=1 f3=0 | Grey title="An lowcut-filtered true-vel RTM image from unmuted data with SRM" $(Axis) screenratio=$(vel0screenratio) | Vppen vpstyle=n scale=$(vel0screenratio) $(dn) out=$R/$@
# smooth vel
rtm%.vel0.mult.lowcut.H: rtm%.vel0.mult.H $B/Light1.x
	<$< $B/Light1.x >$@
rtm%.vel0.mult.lowcut.v: rtm%.vel0.mult.lowcut.H
	<$< Window n3=1 f3=0 | Grey title="An lowcut-filtered RTM image from unmuted data with SRM" $(Axis) screenratio=$(vel0screenratio) | Vppen vpstyle=n scale=$(vel0screenratio) $(dn) out=$R/$@

## from fullmuted data
# true vel
rtm%.truevel0.mult.fullmuted.lowcut.H: rtm%.truevel0.mult.fullmuted.H $B/Light1.x
	<$< $B/Light1.x >$@
rtm%.truevel0.mult.fullmuted.lowcut.v: rtm%.truevel0.mult.fullmuted.lowcut.H
	<$< Window n3=1 f3=0 | Grey title="An lowcut-filtered true-vel RTM image from muted data with SRM" $(Axis) screenratio=$(vel0screenratio) | Vppen vpstyle=n scale=$(vel0screenratio) $(dn) out=$R/$@
rtm%.truevel0.mult.fullmuted.lowcut.spectrum.v: rtm%.truevel0.mult.fullmuted.lowcut.H $B/Fftw2D_program.x
	<$< Rtoc | $B/Fftw2D_program.x flag=1 dim=2 axis1=1 axis2=1 | Cabs | Grey title="The spectrum of a lowcut-filtered true-vel RTM image from fullmuted data with SRM" $(Axis) label1="Kz (/m)" label2="Kx (/m)" $(dn) out=$R/$

# smooth vel
rtm%.vel0.mult.fullmuted.lowcut.H: rtm%.vel0.mult.fullmuted.H $B/Light1.x
	<$< $B/Light1.x >$@
rtm%.vel0.mult.fullmuted.lowcut.v: rtm%.vel0.mult.fullmuted.lowcut.H
	<$< Window n3=1 f3=0 | Grey title="An lowcut-filtered RTM image from muted data with SRM" $(Axis) screenratio=$(vel0screenratio) | Vppen vpstyle=n scale=$(vel0screenratio) $(dn) out=$R/$@

rtm%.vel0.mult.fullmuted.lowcut.spectrum.v: rtm%.vel0.mult.fullmuted.lowcut.H $B/Fftw2D_program.x
	<$< Rtoc | $B/Fftw2D_program.x flag=1 dim=2 axis1=1 axis2=1 | Cabs | Grey title="The spectrum of a lowcut-filtered RTM image from fullmuted data with SRM" $(Axis) label1="Kz (/m)" label2="Kx (/m)" $(dn) out=$R/$


#################################
### RTM without SRM after lowcut filtered
#################################
## from unmuted data
# true vel
rtm%.truevel0.nomult.lowcut.H: rtm%.truevel0.nomult.H $B/Light1.x
	<$< $B/Light1.x >$@
rtm%.truevel0.nomult.lowcut.v: rtm%.truevel0.nomult.lowcut.H
	<$< Window n3=1 f3=0 | Grey title="An lowcut-filtered true-vel RTM image from unmuted data without SRM" $(Axis) screenratio=$(vel0screenratio) | Vppen vpstyle=n scale=$(vel0screenratio) $(dn) out=$R/$@
# smooth vel
rtm%.vel0.nomult.lowcut.H: rtm%.vel0.nomult.H $B/Light1.x
	<$< $B/Light1.x >$@
rtm%.vel0.nomult.lowcut.v: rtm%.vel0.nomult.lowcut.H
	<$< Window n3=1 f3=0 | Grey title="An lowcut-filtered RTM image from unmuted data without SRM" $(Axis) screenratio=$(vel0screenratio) | Vppen vpstyle=n scale=$(vel0screenratio) $(dn) out=$R/$@

## from fullmuted data
# true vel
rtm%.truevel0.nomult.fullmuted.lowcut.H: rtm%.truevel0.nomult.fullmuted.H $B/Light1.x
	<$< $B/Light1.x >$@
rtm%.truevel0.nomult.fullmuted.lowcut.v: rtm%.truevel0.nomult.fullmuted.lowcut.H
	<$< Window n3=1 f3=0 | Grey title="An lowcut-filtered true-vel RTM image from muted data without SRM" $(Axis) screenratio=$(vel0screenratio) | Vppen vpstyle=n scale=$(vel0screenratio) $(dn) out=$R/$@
rtm%.truevel0.nomult.fullmuted.lowcut.spectrum.v: rtm%.truevel0.nomult.fullmuted.lowcut.H $B/Fftw2D_program.x
	<$< Rtoc | $B/Fftw2D_program.x flag=1 dim=2 axis1=1 axis2=1 | Cabs | Grey title="The spectrum of a lowcut-filtered true-vel RTM image from fullmuted data without SRM" $(Axis) label1="Kz (/m)" label2="Kx (/m)" $(dn) out=$R/$

# smooth vel
rtm%.vel0.nomult.fullmuted.lowcut.H: rtm%.vel0.nomult.fullmuted.H $B/Light1.x
	<$< $B/Light1.x >$@
rtm%.vel0.nomult.fullmuted.lowcut.v: rtm%.vel0.nomult.fullmuted.lowcut.H
	<$< Window n3=1 f3=0 | Grey title="An lowcut-filtered RTM image from muted data without SRM" $(Axis) screenratio=$(vel0screenratio) | Vppen vpstyle=n scale=$(vel0screenratio) $(dn) out=$R/$@

rtm%.vel0.nomult.fullmuted.lowcut.spectrum.v: rtm%.vel0.nomult.fullmuted.lowcut.H $B/Fftw2D_program.x
	<$< Rtoc | $B/Fftw2D_program.x flag=1 dim=2 axis1=1 axis2=1 | Cabs | Grey title="The spectrum of a lowcut-filtered RTM image from fullmuted data without SRM" $(Axis) label1="Kz (/m)" label2="Kx (/m)" $(dn) out=$R/$



################################
## get the image showing origins of RTM noises from rtm1.vel0 where sx=3500
################################
rtm1.vel0.mult.fullmuted.imagingsnaps: rtm1.vel0.mult.fullmuted.sousnap1.pdf rtm1.vel0.mult.fullmuted.recsnap1.pdf rtm1.vel0.mult.fullmuted.reflsnap1.pdf \
	rtm1.vel0.mult.fullmuted.sousnap2.pdf rtm1.vel0.mult.fullmuted.recsnap2.pdf rtm1.vel0.mult.fullmuted.reflsnap2.pdf \
	rtm1.vel0.mult.fullmuted.sousnap3.pdf rtm1.vel0.mult.fullmuted.recsnap3.pdf rtm1.vel0.mult.fullmuted.reflsnap3.pdf

## % must be sou, rec or refl
rtm1.vel0.mult.fullmuted.%snap1.v: rtm1.vel0.mult.fullmuted.H ${TXTDIR}/text-imagingsnap1-vel0-mult
	< rtm1.vel0.mult.fullmuted.$*snap.H Window n3=1 f3=12 | Grey $(Axis) screenratio=$(vel0screenratio) $(dn) out=$R/tmp1x.rtm1.vel0.mult.fullmuted.$*snap1.v
	vp_annotate <$R/tmp1x.rtm1.vel0.mult.fullmuted.$*snap1.v batch=y vpstyle=n scale=$(vel0screenratio) text=${TXTDIR}/text-imagingsnap1-vel0-mult >$R/$@
	rm -f $R/tmp1x.rtm1.vel0.mult.fullmuted.$*snap1.v
rtm1.vel0.mult.fullmuted.%snap2.v: rtm1.vel0.mult.fullmuted.H ${TXTDIR}/text-imagingsnap2-vel0-mult
	< rtm1.vel0.mult.fullmuted.$*snap.H Window n3=1 f3=8 | Grey $(Axis) screenratio=$(vel0screenratio) $(dn) out=$R/tmp1x.rtm1.vel0.mult.fullmuted.$*snap2.v
	vp_annotate <$R/tmp1x.rtm1.vel0.mult.fullmuted.$*snap2.v batch=y vpstyle=n scale=$(vel0screenratio) text=${TXTDIR}/text-imagingsnap2-vel0-mult >$R/$@
	rm -f $R/tmp1x.rtm1.vel0.mult.fullmuted.$*snap2.v
rtm1.vel0.mult.fullmuted.%snap3.v: rtm1.vel0.mult.fullmuted.H ${TXTDIR}/text-imagingsnap3-vel0-mult
	< rtm1.vel0.mult.fullmuted.$*snap.H Window n3=1 f3=4 | Grey $(Axis) screenratio=$(vel0screenratio) $(dn) out=$R/tmp1x.rtm1.vel0.mult.fullmuted.$*snap3.v
	vp_annotate <$R/tmp1x.rtm1.vel0.mult.fullmuted.$*snap3.v batch=y vpstyle=n scale=$(vel0screenratio) text=${TXTDIR}/text-imagingsnap3-vel0-mult >$R/$@
	rm -f $R/tmp1x.rtm1.vel0.mult.fullmuted.$*snap3.v

